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INTRODUCTION 


In  recent  years  the  Office  of  Naval  Research  has  sponsored  numerous 
industry  and  university-based  research  projects  aimed  at  furthering  the 
state-of-the-art  of  computer  simulation  methodology.  The  researchers 
performing  that  work  are  among  the  most  distinguished  in  the  field. 

The  numerous  high-quality  journal  and  conference  publications  generated 
by  that  research  attest  to  the  highly  productive  efforts  of  those 
researchers  and  bring  great  credit  to  the  sponsoring  agency. 

But  these  commendable  research  products  have  not  been  sufficiently 
incorporated  into  the  many  computer  simulation  modeling  activities  being 
carried  on  at  various  naval  research  laboratories.  In  general,  naval 
research  simulation  programs  would  receive  high  "marks"  for  their 
programming  structure  and  the  extent  to  which  they  replicate  the 
physical  and  operational  characteristics  of  the  real  systems  that  they 
represent.  But  in  many  cases  they  fail  to  provide  adequate  program 
modules  for  such  statistical  methodology  as  variance  reduction,  statistics 
computation,  the  most  up-to-date  methods  in  random  variate  generation, 
screening  experiments,  experimental  design,  response  surface  methodology, 
optimization  methods,  and  other  proven  statistical  techniques  in 
computer  simulation.  Efforts  from  one  naval  laboratory  to  another  are 
very  "uneven",  and  lack  a  consistent  modeling  structure  and  methodo¬ 
logical  approach. 

This  research  has  addressed  the  development  of  a  simulation  model¬ 
ing  language  that  provides  a  more  consistent  modeling  structure  and 
integrates  the  most  up-to-date,  proven  statistical  methodology  and 
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optimization  approaches.  It  has  sought  to  "marry"  the  most  applicable 
results  from  past  and  ongoing  statistical  methodology  research  with 
the  most  promising  computer  simulation  modeling  approaches.  This 
"marriage"  takes  the  form  of  a  generalized  computer  simulation  program 
that : 

•  Allows  the  simulation  modeler  to  develop  models  of  specific 
operations  with  minimum  programming. 

•  Allows  for  statistical  methodology  features  with  simple  program 
statements . 

•  Provides  a  more  consistent  modeling  approach  from  one  laboratory 
to  the  next,  and  thus  allow  portability  of  models  as  needed. 

This  report  describes  a  generalized  computer  simulation  language 
called  NAVMAP  (Naval  Modeling  and  Analysis  Program) .  This  language  is 
FORTRAN-based,  so  as  to  maximize  its  potential  utility  in  naval 
laboratories.  It  enables  discrete-event,  continuous,  and  state-event 
computer  simulation  modeling  in  an  event-oriented  programming  structure. 
It  is  designed  to  be  compact  and  portable,  with  potential  for  imple¬ 
mentation  on  minicomputers. 
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RESEARCH  OBJECTIVES 

This  research  will  require  three  years  in  total,  and  involve  the 
following  tasks: 

1.  Collection  and  analysis  of  appropriate  ref erences  on  (a)  sta¬ 
tistical  methodology  and  (b)  simulation  modeling  languages 
that  are  needed  for  the  research. 

2.  Design  of  a  simulation  modeling  structure  that  incorporates 
continuous  simulation,  discrete-event  simulation,  statistical 
methodology  and  optimization  techniques. 

3.  Coding  and  testing  individual  program  modules. 

4.  Evaluation  of  the  program  package  with  selected  naval  simu¬ 
lation  models. 

5.  Transmittal  of  the  research  products  to  naval  laboratories. 

The  first  two  phases  of  the  project,  concluding  on  June  30,  1981 

and  covered  in  this  report,  essentially  involves  the  first  three  tasks 
stated  above. 
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RESEARCH  APPROACH 

This  research  project  involved  two  main  elements,  in  terms  of  a 
computer  simulation  modeling  language  for  naval  simulations.  These  two 
elements  were  (1)  statistical  methodology  and  (2)  a  simulation  modeling 
and  programming  structure.  The  following  sections  discuss  the  principal 
issues  embedded  in  each  of  these  elements. 

Statistical  Methodology  in  Computer  Simulation 

A  computer  simulation  model,  particularly  one  involving  stochastic 
elements,  must  provide  for  the  following  functional  capabilities: 

1.  Random  number  generation 

2.  Random  variate  generation 

3.  Statistics  collection  and  reporting 

4.  Variance  reduction 

5.  Input  analysis 

6.  Output  analysis 

7.  Experimental  design 

8.  Optimization 

The  first  four  of  these  functional  capabilities  are  internal  to  the 
computer  simulation,  and  must  be  made  available  within  the  simulation 
modeling  structure.  For  example,  in  a  naval  simulation  model  such  as 
SPEARS  [28],  it  might  be  necessary  to  record  the  range  at  which  "kills" 
of  incoming  enemy  missiles  occur.  If  we  assume  that  an  event  MKILL 
gives  rise  to  such  a  "kill",  and  that  the  kill  distance  XDISK  is  computed 
at  the  instant  in  which  the  "kill"  occurs,  a  single  statement  can  be 
used  to  record  this  value,  such  as  follows: 


CALL  TALLY  (4, XDISK) 
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This  statement  calls  a  subroutine  TALLY  which  updates  the  statistics 
for  variable  4,  which  the  modeler  has  designated  as  XDISK,  the  missile 
kill  distance.  All  statistics  collection  and  summary  operations  would 
then  be  effected  automatically  through  subroutine  TALLY.  Those  statis¬ 
tical  operations  that  are  internal  to  the  operation  of  the  simulation 
model  can  be  incorporated  directly  into  the  simulation  modeling  language, 
allowing  the  modeler  to  call  upon  powerful  statistical  capability  with 
a  mere  handful  of  program  statements.  This  research  has  investigated 
the  most  efficient  techniques  for  accomplishing  these  internal  sta¬ 
tistics  operations  represented  by  functional  capabilities  1  through  4. 

Functional  capabilities  5  through  8  above  take  the  form  of  programs 
that  lie  external  to  the  computer  simulation  model.  For  instance, 
"canned"  programs  would  be  prepared  to  allow  ”goodness-of-f it"  testing 
to  sample  data,  preparatory  to  having  the  modeler  identify  a  distribution 
and  parameter  set  for  a  specified  random  variable  in  the  model.  At 
least  two  statistical  procedures,  the  Chi-Square  test  and  the  Kolmogorov- 
Smirnov  test,  will  be  coded  in  the  next  phase  of  research  for  the 
following  probability  distributions: 


1. 

Exponential 

7. 

Uniform 

2. 

Normal 

8. 

Poisson 

3. 

Lognormal 

9. 

Binomial 

4. 

Gamma 

10. 

Geometric 

5. 

Weibull 

11. 

Negative  Binomial 

6. 

Beta 

12. 

Hypergeometric 

Similar  "canned"  programs  will  be  prepared,  or  identified,  for  such  sta¬ 
tistical  procedures  as  multiple  regression,  analysis  of  variance,  and 
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analysis  of  covariance.  These  statistical  programs  must  be  able  to 
interface  with  the  simulation  language  described  here.  Programs  for 
optimization  of  computer  simulations,  such  as  those  by  Biles  [2]  and 
Smith  [51-53],  will  be  modified  to  interface  naturally  with  the  simulation 
modeling  structure  described  in  this  report. 

Computer  Simulation  Modeling  Structure 

The  development  of  a  generalized  computer  simulation  modeling 
structure  for  naval  simulation  has  involved  the  following  subtasks: 

1.  Identification  of  the  common  features  of  naval  simulation 
models,  through  interviews  with  Navy  laboratory  personnel 
and  analysis  of  existing  naval  simulation  models  (for 
example  [1,  28]). 

2.  Selection  of  a  base  programming  language  that  is  compatible 
with  the  maximum  number  of  naval  laboratory  computers  (for 
example,  FORTRAN  77). 

3.  Formulation  of  modeling  approaches  for 

a.  Continuous  simulation 

b.  Discrete-event  simulation 

c.  State-event  simulation 

d.  Combined  continuous/discrete-event/state-event 
simulation. 

4.  Formulation  of  the  basic  modeling  structure  for  the  overall 
language. 

5.  Coding  the  generalized  simulation  language. 

6.  Evaluation  and  testing  of  the  language. 
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The  crucial  feature  of  simulation  modeling  structure  is  the  adoption 
of  a  "view  of  the  world"  from  which  the  modeler  operates.  Research  to 
date  in  reviewing  naval  simulation  modeling  activity  has  revealed  two 
major  simulation  approaches  in  use  by  naval  personnel:  (a)  digital 
simulation,  in  which  the  entire  simulation  logic  is  represented  via  a 
computer  program;  and  (b)  real-time  simulation,  in  which  there  is  physical 
hardware,  such  as  a  torpedo  or  an  aircraft,  in  the  loop.  This  research 
has  focused  on  digital  computer  simulation,  although  some  of  the  soft¬ 
ware  that  lies  external  to  the  simulation  model  might  very  well  be 
applicable  with  real-time  simulation. 

A  second  aspect  of  the  simulation  modeling  structure  is  discrete 
simulation  versus  continuous  simulation.  Discrete  simulation  occurs 
when  the  dependent  variables  change  discretely  at  specified  points  in 
simulated  time;  for  example,  when  a  mine  detonates,  the  number  of  mines 
in  a  minefield  is  decremented  by  one  unit.  In  continuous  simulation 
the  dependent  variables  may  change  continuously  over  simulated  time; 
for  instance,  the  position  of  a  torpedo  in  a  three-dimensional  space 
representing  the  ocean.  Naval  simulation  modeling  definitely  requires 
both  discrete  and  continuous  simulation.  The  generalized  simulation 
language  developed  during  this  research  project  affords  both  discrete 
and  continuous  simulation  capability. 

A  third  aspect  of  the  simulation  modeling  structure's  "view  of  the 
world"  is  whether  the  discrete  simulation  adopts  (a)  an  event  orienta¬ 
tion,  (b)  an  activity  scanning  orientation,  or  (c)  a  process  orientation. 

The  NAVMAP  language  described  in  this  report  adopts  an  event 
orientation  for  the  discrete  simulation.  In  this  orientation,  events 
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occur  at  instants  in  time.  Events  alter  the  state  of  the  system,  and 
these  state  changes  must  be  so  recorded  and  the  appropriate  statistics 
collected.  Thus,  the  modeler  has  the  task  of  preparing  FORTRAN  sub¬ 
routines  which  execute  the  logical  structure  of  the  events  which  describe 
the  behavior  of  the  system  being  modeled.  NAVMAP  automatically  calendars 
the  events  in  proper  sequence,  and  advances  time  from  one  event  to  the 
next.  A  NAVMAP  model  of  a  given  naval  system  would  typically  consist 
of  a  MAIN  program  and  several  event  subroutines.  The  language  is 
described  in  the  ensuing  sections. 
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NAVMAP 

The  current  status  of  this  project,  after  almost  two  years  of 
effort  under  Contract  N00014-79-C-0757 ,  is  that  a  combined  continuous/ 
discrete-event/state-event  simulation  language  called  NAVMAP  (Naval 
Modeling  and  Analysis  Program)  has  been  coded  in  FORTRAN  and  evaluated. 

The  objective  was  that  this  simulation  language  should  be  highly  flexible, 
compact,  and  portable.  Its  flexibility  is  verified  by  the  fact  that 
it  can  be  employed  in  either  discrete-event  simulation,  continuous 
simulation,  state-event  simulation,  or  a  combination  of  these.  Likewise, 
it  can  be  interfaced  with  FORTRAN-based  statistical  and  optimization 
programs.  The  compactness  of  this  naval  systems  simulation  language 
is  evident  in  the  fact  that  it  contains  only  2600  lines  of  code,  far 
less  that  the  GASP-IV  language  to  which  it  is  comparable  in  flexi¬ 
bility  and  capability.  The  probability  of  the  language  lies  in  the 
fact  that  it  is  FORTRAN-based  and  compact,  which  makes  it  usable  on  any 
computer  with  a  FORTRAN  compiler,  including  most  minicomputers. 

Elements  of  Design 

In  the  design  of  NAVMAP,  several  objectives  have  been  considered. 
Primary  of  those  are:  (1)  flexibility  for  the  interfacing  with 
statistical  and  optimization  procedures,  (2)  portability,  and  (3) 
compactness.  In  this  section,  each  of  these  goals  and  the  steps  for 
achieving  them  are  discussed  separately.  Finally,  in  the  last  part 
of  this  section,  the  data  structure  of  the  language  is  put  into  per¬ 
spective  and  its  features  are  described. 

Interfacing  Flexibility 

Most  of  the  simulation  languages,  like  SLAM  [41]  and  GASP  [40], 
provide  a  summary  report  at  the  end  of  each  simulation  run.  One  could 
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obtain  the  mean,  standard  deviation,  and  other  statistical  measures  on 
each  of  a  number  of  model  variables,  as  well  as  histograms  and  plots  of 
continuous  variables.  However,  more  often  than  not  these  results  are 
useless  with  respect  to  statistical  analysis.  To  make  this  point 
clear,  let  us  proceed  with  a  very  simple  example. 

Consider  an  M/M/1  queueing  system,  and  suppose  it  is  desired  to 
estimate  the  steady  state  mean  waiting  time  of  the  customers.  To 
accomplish  this,  observation  is  made  on  the  waiting  time  of  each 
customer  entering  the  system  and  the  time  he  completes  the  service. 

It  is  well  known  that  this  observation  is  correlated  and  therefore 
the  standard  deviation  provided  in  the  statistical  summary  report  is 
useless.  To  obtain  histograms,  one  must  specify  the  parameters  of 
the  histogram  prior  to  the  run.  Now  suppose  one  has  made  an  expensive 
simulation  run  to  obtain  a  histogram,  and  after  the  run  realizes  that 
the  histogram  is  not  well  shaped  because  of  ill  specification  of  the 
parameters.  To  obtain  a  new  histogram,  one  must  completely  replicate 
the  expensive  simulation  run.  This  could  also  be  true  when  a  plot  of  the 
values  of  a  variable  is  to  be  obtained. 

In  the  event  of  correlated  observations,  one  might  want  to  try 
the  batch  means  procedure  [30-34],  where  one  needs  to  have  access  to  all 
the  observations  produced  during  the  simulation.  Even  in  the  situation 
where  more  advanced  statistical  analysis  is  needed,  such  as  regression 
or  analysis  of  variance,  access  to  individual  observation  might  be 
required.  The  issue  illustrated  here  is  that  one  essentially  needs 
to  have  access  to  individuax  observations  in  order  to  minimize  number 


of  runs  for  analysis.  We  take  the  view  that  the  statistical  analysis 
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and  inferences  about  the  unknown  parameters  of  the  problem  under  study 
must  be  a  separate  module  from  the  simulation  model.  That  is,  the  sole 
purpose  of  the  simulation  must  be  collection  of  observations  on  variables, 
and  the  job  of  analysis  and  preparation  of  data  must  be  assigned  to 
other  modules.  Obviously,  the  existing  simulation  languages  have  short¬ 
comings  in  this  regard.  In  NAVMAP,  we  have  achieved  this  goal  by 
letting  the  user  specify  a  unit  number  for  each  of  the  statistics.  This 
results  in  recording  each  individual  observation  on  that  unit. 

There  are  three  types  of  statistics  collection  devices  in  NAVMAP. 

The  particular  type  employed  in  a  given  instance  depends  on  the  nature 
of  the  model  variable  for  which  statistics  are  being  collected.  The 
first  type  of  statistic  is  a  Tally  statistic,  which  collects  information 
on  an  observation  of  a  non-time  dependent  random  variable.  A  random 
variable  is  considered  non-time  dependent  when  its  value  at  a  precise 
instant  in  simulated  time  is  recorded  without  regard  either  to  the 
length  of  time  that  value  has  persisted  or  the  value  it  has  changed  from. 
For  example,  if  one  is  collecting  observations  on  the  "kill  distance" 
XDISK  at  which  an  incoming  enemy  missile  is  intercepted,  this  value  is 
recorded  at  the  instant  in  time  the  "kill"  occurs,  which  is  marked 
by  the  event  MKILL.  In  the  simulation  of  an  engagement,  observations 
of  XDISK  will  be  recorded  each  time  a  kill  event  occurs.  If  XDISK 
is  the  fourth  variable,  the  observation  would  be  recorded  with  the 
statement  in  MKILL: 

CALL  TALLY  (4, XDISK) 

NAVMAP  automatically  computes  the  minimum,  maximum,  arithmetic  mean, 
standard  deviation  and  number  of  observations  of  each  Tally  variable 


when  one  of  the  functions  TMIN,  TMAX,  TAVG,  TSTD  or  TNUM  is  accessed. 
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This  information  is  automatically  computed  and  printed  at  the  conclu¬ 
sion  of  the  simulation  run.  After  all  simulation  runs  have  been 
executed,  the  user  can  access  the  accumulated  observations  by  reading 
the  data  from  the  selected  storage  unit  identified  on  the  Tally  data 
card  (see  Appendix  A) .  The  symbol  **  is  used  to  separate  data  between 
simulation  runs.  These  observations  are  recorded  unformatted. 

Discrete  time  persistent  statistics  are  employed  for  observing 
the  values  of  either  XX  variables  or  number  of  entities  of  a  file  as 
they  change  and  recording  their  value  and  the  time  of  their  change. 

If  the  unit  number  is  specified,  the  value  and  the  time  of  the  change, 
respectively,  would  be  dumped  on  that  unit.  Again,  these  are  written 
unformatted  and  **  is  used  to  separate  observations  of  each  run. 

DMIN,  DMAX,  DA7G,  DSTD,  and  DPRD  are  functions  that  respectively  give 
minimum,  maximum,  average,  standard  deviation,  and  the  length  of  the 
period  of  the  •ooservations  once  they  are  called. 

Continuous  time  persistent  statistics  do  behave  exactly  as  the 
discrete  time  persistent  except  they  are  for  collecting  observations 
on  SS,  DD  variables  (continuous  variables).  Statistical  collection  can 
be  truncated  by  several  means.  Subroutine  CLEARS  clears  all  the  statis¬ 
tical  arrays  once  it  is  called.  Sometimes  it  is  desirable  to  clear 
some  variables  at  some  specified  time  and  clear  others  at  some  other 
time.  This  can  be  achieved  by  calling  subroutines  TRUND(N) ,  TRUNC(N), 
and  TRUNT(N).  One  can  also  clear  the  statistical  arrays  at  a  pre¬ 
determined  time  by  specifying  this  time  on  the  Initialization  Card. 

Another  means  of  clearing  arrays  is  to  judiciously  discard  some  observation 
according  to  one  of  several  proposed  methods  (for  example,  Schruben[48] ) 
once  the  simulation  runs  are  completed. 
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Portability 

Simulation  languages  which  feature  real  valued  attributes  like 
GASP  IV[40]  and  SLAM  [41],  although  claimed,  are  not  totally  portable. 

The  problem  arises  because  of  the  data  structure  for  maintaining  the 
entities.  Each  entity  consists  of  some  integer  information,  such  as 
pointers,  and  some  real  information,  such  as  attributes.  In  these 
languages  there  are  two  large  arrays  that  are  set  EQUIVALANCE:  one 
is  for  maintaining  integer  information  and  the  other  for  maintaining 
the  real  information  associated  with  each  file  entry. 

On  some  computers  the  number  of  integer  words  per  number  of  real 
words  is  not  1.  That  is,  the  real  word  is  comprised  of  4  bytes  and  the 
integer  word  is  comprised  of  2  bytes.  In  these  machines,  a  problem  arises 
if  two  arrays,  one  integer  and  one  real,  of  the  same  length  are  set 
EQUIVALANCE.  To  illustrate  this  difficulty,  suppose  the  integer  array 
is  ISET(12)  and  the  real  array  is  RSET(12).  Then  as  Figure  1  shows, 
location  2  of  ISET  does  not  correspond  to  location  2  of  RSET,  and 
location  7  of  RSET  does  not  correspond  to  any  location  in  ISET. 

Now  let  NIR  =  number  of  integer  words  per  real  word 

MSET  =  dimension  of  ISET;  a  multiple  of  NIR 

MSET 

Then  clearly  the  dimension  of  RSET  must  be  .  If  L  is  an  integer 

address,  the  corresponding  real  address  should  be  ^“^xR  ~  *  F°r  t*le 
above  example:  NIR  =  2,  MSET  =  10,  and  the  relationship  between  the 
integer  and  real  address  is  clear  from  Figure  2. 

The  problem  can  be  resolved  by  allocating  the  upper  part  of  the 
array  to  integer  values  and  the  lower  part  of  the  array  to  the  real 
values  and  have  a  pointer  in  the  integer  part  to  give  the  address  of 
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Integer  part 
of  entity  1 


Integer  part 
of  entity  2 


Real  part  of  entity  2 


Real  part  of  entity  1 


Figure  2.  Filing  Structure  in  NAVMAP 
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the  corresponding  real  part.  To  illustrate  this  concept,  suppose  every 
entity  is  comprised  of  4  integer  variables  and  1  real  variable,  where 
the  last  information  in  the  integer  part  is  the  pointer  to  the  location 
of  the  real  part.  Then  as  Figure  2  indicates  the  problem  is  alleviated, 
thus  affording  a  simple  basis  for  designing  a  completely  portable  filing 
structure. 

Compactness 

In  a  given  simulation  model  one  might  or  might  not  need  files, 
tallies,  discrete  time  persistant  statistics,  continuous  time  persistent 
statistics,  and  counters.  In  existing  simulation  languages  fixed 
amount  of  storage  are  allocated  for  all  of  the  above  regardless  of 
whether  they  are  used.  We  take  the  view  that  this  procedure  is  extrava¬ 
gant  with  respect  to  storage  space. 

We  observe  that  each  of  the  above  data  elements  integer  part  infor¬ 
mation  and  real  part  information.  We  can  organize  the  integer  information 
in  one  block  and  the  real  information  in  another  block,  and  maintain  a 
pointer  to  relate  the  integer  part  to  real  part.  If  we  keep  a  pointer 
to  the  location  of  the  first  integer  block  of  one  type  (say  files),  we 
can  have  a  dynamic  data  structure  for  that  type  of  data  element.  In 
the  following  paragraphs,  we  define  the  integer  and  real  information 
for  each  type  of  data  element  and  the  manner  that  they  have  been  organ¬ 


ized  in  NAVMAP . 


(a)  Entries 


Integer 

Real 

1. 

Location  of  predecessor 

1. 

Ranking  value 

2. 

Location  of  successor 

2. 

Attribute  1 

3. 

Event  Code 

3. 

Attribute  2 

4. 

Pointer  to  attributes 

• 

N+l 

•  Attribute  N 

Files 

Integer 

1. 

Number  of  entries  in  the  file 

2. 

Location  of  first  entry 

3. 

Location  of  last  entry 

4. 

Ranking  rule 

5. 

Ranking  attribute 

6. 

Location  of  the  corresponding 
statistics  block 

discrete  time  persistent 

Tally 

Integer 

Real 

1. 

Location  of  statistical  block 

1. 

EX 

2. 

Label 

2. 

2 

EXZ 

3. 

Output  device  number 

3. 

Number  of  observation 

4. 

Min. 

5 .  Max . 
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(d)  Discrete  time  persistent  statistics 


Integer 

Real 

1. 

Location  of  statistical  block 

1. 

/Xdt 

2. 

Label 

2. 

/X2dt 

3. 

Variable  Type 

3. 

XLAST 

4. 

Output  device  number 

4. 

MIN 

5. 

MAX 

'* 

6. 

TLAST 

7. 

TCLEAR 

Continuous  time  persistent  statistics 

Integer 

Real 

1. 

Location  of  statistical  block 

1. 

/  Xdt 

2. 

Label 

2. 

Jx2dt 

3. 

Variable  type 

3. 

XLAST 

4. 

Output  device  number 

4. 

MIN 

5. 

MAX 

6. 

TLAST 

7. 

TCLEAR 

(f )  Counters 

_ Integer _ 

1.  Current  count 

2.  Label 

3.  Limit 

This  structural  feature  of  NAVMAP  enables  a  highly  compact  storage  of  a 


simulation  model. 
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Organization  of  NAVMAP 

NAVMAP  is  organized  in  much  the  same  way  as  GASP-IV[40],  but 
with  quite  different  subprograms.  Figure  3  shows  the  organization  of 
NAVMAP.  Table  1  gives  the  subroutines  used  for  executive  control  and 
statistics  collection.  Table  2  shows  the  function  subprograms  employed 
for  statistics  collection.  Table  3  presents  the  random  process  generation 
functions.  Table  4  lists  the  key  NAVMAP  variables. 


Description  of  Routines 
1 )  Main 

The  user  has  the  choice  to  write  his  own  main  program  or  use 
the  default  program  provided  in  NAVMAP.  In  this  routine,  the  length 
of  ISET  and  RSET  are  set  and  four  variables  NCROR,  NPRTR,  LFI ,  and 
LLR  are  initialized.  Subroutine  SOAP  is  then  called.  LFI  must 
always  be  set  equal  to  1.  LLR  must  always  be  set  equal  to  the  length 
of  RSET  plus  1.  NCRDR,  and  NPRTR  are  respectively  number  of  card 
reader  unit  and  printer  unit  at  the  computer  installation.  A 
typical  main  program  is  given  below.  Note  that  the  common  blocks 
are  to  be  written  exactly  as  they  are  in  the  sample  main  program. 


DIMENSION  RSET(5000) 

COMMON /LOC/ LFI , LLR, LFFB , LFTB , LFDSB , LFCSB , LFCTB ,LFAE , LCAL 
COMMON  ISET (5000) 

C0MM0N/GSC1 /NCRDR, NPRTR, SS (99) ,DD(99) ,ATRIB(99) ,SSL(99) ,JJ, 
-DDL (99) ,TN0W,  XX(99) ,DTN0W, ISTOP 
EQUIVALENCE  ( ISET ( 1) , RSET (1) ) 

NCRDR=5 

NPRTR=6 

LFI=1 

LLR=5001 

CALL  SOAP 

STOP 

END 
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FILEM 

CLEARS 

DAVG 

LOCATE 

DRAND 

RMOVE 

TRUND 

DSTD 

ATR 

EXPON 

COPY 

TRUNC 

DPRD 

LPRED 

UNFRM 

SCHD 

TRUNT 

DMIN 

LSUCC 

WEIBL 

LANGUAGE 

TALLY 

TAVG 

DMAX 

LFE 

TRIAG 

LIBRARY 

COUNT 

TSTD 

CAVG 

LLE 

RNORM 

ROUTINES 

TNUM 

CSTD 

NQ 

RLOGN 

TKIN 

CSTD 

ERLING 

TMAX 

CPRD 

CAMA 

CMIN 

BETA 

CMAX 

NPSSN 

IIH3 
CD— CZI 


SUBROUTINE  I  CALLS  J 
SUBROUTINE  I  IS  ACCESSIBLE  TO  J 


Figure  3.  Organization  of  the  NAVMAP  Language 
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Subroutine _ 

FILEM  (IFILE,  A) 

RMOVE  (I,  IFILE,  A) 

COPY  (I,  IFILE,  A) 

LOCATE  (I,  IFILE,  LENT) 

SCHD  (IEVENT,  TIME,  A) 

TALLY  (N,  VAR) 

COUNT  (N,  INC) 

CLEARS 

TRUND(N) 

TRUNC(N) 

TRUNT(N) 

INTLC 

EVENT (I) 

OUTPUT 

STATE 

Table  1. 


Description _ 

Files  the  entity  with  attributes  in 
array  A  in  File  IFILE. 

Removes  Ith  entry  of  the  file  IFILE 
and  loads  its  attributes  into  array  A. 

Copies  the  attributes  of  the  Ith 
entry  of  file  IFILE  into  array  A. 

Gives  the  location  of  the  Ith  entry 
of  file  IFILE  in  LENT. 

Schedules  event  IEVENT  at  time  TIME, 
with  attributes  stored  in  vector  A. 

Collects  statistics  on  Tally  #N 
using  VAR. 

Increments  counter  If N  with  INC. 

Clears  statistical  arrays,  at  TNOW. 

Truncates  discrete  time  persistent 
statistics  #N  at  TNOW. 

Clears  continuous  time  persistent 
statistics  If N  at  TNOW. 

Clears  Tally  #N  statistics  at  TNOW. 

User  written  for  setting  initial 
conditions  of  the  simulation. 

User  written  for  defining  the  events. 

User  written  for  special  output 
preparation  or  any  adhoc  procedure. 

It  is  called  after  the  end  of  each 
simulation  run. 

User  written  subroutine  for  defining 
the  differential  and  difference 
equations . 


Simulation  Executive  Subroutines 
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Functions _ 

TAVG(N) ,  TSTD(N) ,  TNUM(N) , 
TMIN(N) ,  TMAX(N) 


DAVG(N) ,  DSTD(N) ,  DPRD(N) , 
DMIN(N),  DMAX(N) 


CAVG(N) ,  CSTD(N),  CPRD(N) , 

CMIN(N) ,  CMAX(N) 

ATR(LENT,  I) 

LPRED  (LENT),  LSUCC  (ENT) 

NQ  (IFILE) ,  LFE  (IFILE) , 

LLE  (IFILE) 

Table  2.  Statistics  Collection 


Description _ 

Gives  average,  standard  deviation, 
number  of  observations,  minimum, 
and  maximum  of  Tally  //N. 

Gives  average,  standard  deviation, 
length  of  period  of  statistical 
collection,  minimum,  and  maximum  of 
discrete  statistical  time  persistent 
#N. 

Gives  average,  standard  deviation 
length  of  period  of  statistical 
collection,  minimum,  and  maximum  of 
continuous  statistical  time  persis¬ 
tent  #N . 

Gives  the  value  of  the  Ith  attribute 
of  the  entry  located  at  LENT. 

Gives  the  predecessor  and  successor 
of  the  entry  located  at  LENT. 

Gives  the  number  of  entries,  loca¬ 
tion  of  first  entry,  and  location 
of  the  last  entry  of  file  IFILE. 


and  Reporting  Functions 
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Functions _ 

DRAND(IS) 

EXPON (XMEAN , IS) 

UNFRM(ULO , UHI , IS) 

UEIBL(BETA, ALPHA, IS) 

TRI AG ( XLO , XMODE ,XHI , IS) 

RNORM(XMN, STD, IS) 

RLOGN(XMN,STD,IS) 

ERLNG(EMN,XK, IS) 

GAMA ( B ETA , ALPHA , I S ) 

BETA(THETA, PHI , IS) 

NPSSN(XMN.IS) 

Table  3.  Random  Process 


Description _ 

A  pseudo-random  number  oftained  from 
random  number  stream  IS. 

A  sample  from  an  exponential  distri¬ 
bution  with  mea?  XMEAN  using  random 
number  stream  IJ. 

A  sample  from  a  uniform  distribution 
in  the  interval  ULO  to  UHI  using 
random  number  stream  IS. 

A  sample  from  a  Weibull  distribution 
with  scale  parameter  BETA  and  shape 
parameter  ALPHA  using  random  number 
stream  IS. 

A  sample  from  a  triangular  distri¬ 
bution  in  the  interval  XLO  to  XHI 
with  mode  XMODE  using  random  number 
stream  IS. 

A  sample  from  a  normal  distribution 
with  mean  XMN  and  standard  deviation 
STD  using  random  number  stream  IS. 

A  sample  from  a  lognormal  distribu¬ 
tion  with  mean  XMN  and  standard 
deviation  STD  using  random  number 
stream  IS. 

A  sample  from  a  Erlang  distribution 
which  is  the  sum  of  XK  exponential 
samples  each  with  mean  EMN  using 
random  number  stream  IS. 

A  sample  from  a  gamma  distribution 
with  parameters  BETA  and  ALPHA  using 
random  number  stream  IS. 

A  sample  from  a  beta  distribution 
with  parameters  THETA  and  PHI  using 
random  number  stream  IS. 

A  sa^.-ie  from  a  Poisson  distribution 
with  mean  XMN  using  random  number 
stream  IS. 


Generation  Functions 


User  Variables 

NCRDR 

NPRTR  \ 

TNOW 

DTNOW 

ATRIB(I) 

SS(I) 

SSL(I) 

DD(I) 

DDL (I) 

XX(I) 

JJ 

ISTOP 


Description _ 

Installation  card  reader  unit  # 
Installation  line  printer  unit  if 
Current  time 
The  step  size 

Attribute  I  of  current  entity 

Value  of  state  variable  I  at  TNOW 

Value  of  state  variable  I  at  TLAST 

Value  of  the  derivative  of  state 
variable  I  at  TNOW 

Value  of  the  derivative  of  state 
variable  I  at  TLAST 

Value  of  global  variable  I 

An  integer  global  variable 

An  integer  variable  that  terminates 
simulation  when  it  is  equal  to  1 


Table  4.  Key  Simulation  Language  Variables 
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2)  SOAP 

SOAP  is  the  executive  routine  in  NAVMAP. 

3)  DEFALT 

DEFALT  is  used  for  setting  the  default  values  for  NAVMAP 
variables . 

4)  DATAIN 

DATAIN  is  used  for  reading  the  input  cards  and  printing  data 
in  the  way  they  are  read. 

5)  SETFB ,  SETCAL,  SETDS ,  SETCS,  SETCT  and  SETENP 

These  subroutines  are  employed  for  setting  up,  respectively, 
a  file  block,  the  calendar,  the  discrete  time-persistent  block, 
the  continuous  time  persistent  block,  a  counter,  and  the  entry 
pool. 

6)  CLEARF,  CLEARV  AND  INTAL 

These  subroutines  are  used  for  clearing  files  and  variables, 
and  setting  up  some  of  the  initial  conditions  prior  to  each  run. 

7)  EXEC 

EXEC  is  used  for  the  time  advance  mechanism  and  for  solving 
differential  and  difference  equations. 

8)  INTLC 

INTLC  is  a  user  written  subroutine  for  setting  up  the  initial 
conditions  of  each  run. 

9)  DISC  and  EVENT 

These  are  used  for  processing  events.  EVENT  is  a  user 


written  subroutine. 


26 


10)  SSCND  and  KKRSS 

These  subroutines  are  related  to  state-events  and  for  deter¬ 
mining  whether  or  not  a  state-event  has  occurred  within  the  speci¬ 
fied  tolerance. 

11)  STATE 

Subroutine  STATE  is  developed  by  the  user  for  specifying  the 
differential  and  difference  equations  which  describe  the  continuous 
component  of  the  model. 

12)  ECHO  and  SUMRY 

These  subroutines  are  used  respectively  for  printing  the  Echo 
and  Summary  Reports. 

13)  OUTPUT 

This  is  a  user  written  routine  used  for  special  treatment  of 
data  at  the  end  of  each  simulation  run. 

14)  Language  Library  Routines 

Any  of  the  library  function  subprograms  available  on  a  FORTRAN 
compiler  are  accessible  in  NAVMAP. 

Data  input  to  NAVMAP  consists  of  punched  cards.  The  input  data 
is  arranged  in  thirteen  card  types,  depending  on  the  particular  simulation 
function  that  the  data  supports  or  initializes.  Appendix  A  describes 
the  data  formats  for  these  thirteen  types  of  cards. 

As  stated  above,  NAVMAP  possesses  the  capability  for  discrete- 
event,  continuous,  and  state-event  simulation.  Appendix  B  gives  the 


I 

I 
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data  input,  user  program  listing,  and  simulation  output  for  three  rela¬ 
tively  simple  problems  that  demonstrate  the  three  modes  of  simulation 
capability.  Problem  1  is  a  simple  M/M/1  queue;  that  is,  a  single¬ 
channel  queueing  problem  with  Poisson  arrivals  and  exponential  service 
times.  The  first  page  of  the  printout  shows  the  data  input  cards.  The 
next  two  pages  show  the  user-written  FORTRAN  subroutines  that  describe 
the  M/M/1  queueing  system.  The  two  event  subroutines  in  this  discrete- 
event  model  are  ARVL  and  ESRV,  which  correspond  to  the  arrival  of  an 
entity  into  the  system  and  the  completion  of  service,  respectively. 

Next  is  shown  the  printout  of  the  input  data,  which  provides  the  analyst 
a  check  on  the  verity  of  his  input  to  the  model.  Then  two  pages  give 
the  Echo  Report,  which  shows  all  data  and  system  conditions  at  the 
outset  of  the  execution  of  the  model.  Finally,  a  Summary  Report  gives 
the  final  results  of  the  simulation. 

Problem  2  in  Appendix  B  gives  input  data,  user  subprograms  and 
output  results  for  a  sample  continuous  simulation  problem,  the  Pilot 
Ejection  Probl.  .rom  Pritsker  and  Pegden[41].  Key  features  here  are 
the  initialization  of  state  variables  SS(I),  1=1,..., 4  in  Subroutine 
INTLC  and  representation  of  the  state  equations  in  Subroutine  STATE. 
These  equations  model  the  behavior  of  a  pilot  ejection  module  (pilot 
and  seat)  which  is  thrown  from  a  crippled  aircraft.  The  simplicity  of 
the  user  requirements  for  a  continuous  simulation  model  is  illustrated 
here. 

Problem  3  in  Appendix  B  gives  input  data,  user  subprograms  and 


output  results  for  a  sample  continuous/state-event  simulation,  the  Cedar 
Bog  Lake  Problem  from  Pritsker  and  Pegden[41J. 


28 


Appendix  C  gives  a  complete  listing  of  NAVMAP.  The  2607  lines  of 
program  code  represent  a  significant  compression  from  that  of  the  GASP-IV 
language[40]  to  which  it  is  comparable  in  capability. 

Future  efforts  will  concentrate  on  three  areas: 

1.  The  extension  of  NAVMAP  to  include  network  modeling. 

2.  The  development  of  FORTRAN-based  statistical  and  optimization 
programs  that  would  enable  preliminary  data  analysis,  regression, 
correlation,  analysis  of  variance  and  response  surface  analysis. 

3.  Demonstration  of  NAVMAP  with  realistic  naval  systems,  including 
the  following: 

a.  Undersea  systems  vis  a  vis  the  NtJSC  laboratories  at 
Newport,  RI  and  New  London,  CT. 

b.  Minefield  systems  vis  a  vis  the  ARL  at  The  Pennsylvania 
State  University  and  the  Naval  Weapons  Laboratory  at 
Dahlgren,  VA. 

Finally,  NAVMAP  will  be  completely  documented  and  transmitted  to  selected 
naval  laboratories  in  the  form  of  magnetic  tapes. 


REFERENCES 


29 


1.  Bailey,  J.L. ,  and  M.  A.  Weedon,  "The  Minefield  Analysis  with  Hunting 
Evaluation  Model  (MAYHEM)",  NWL  Technical  Report  TR-2742,  U.S.  Naval 
Weapons  Laboratory,  Dahlgren,  Virginia,  May  1972. 

2-  Biles,  W.  E. ,  Optimization  of  Multiple-Response  Simulation  Models",  Final 
Report  to  the  Office  of  Naval  Research,  Contract  N00014-76-C-1021,  NR 
277-237,  University  of  Notre  Dame,  October  1978. 

3.  Biles,  W.  E. ,  "A  Gradient-Regression  Search  Procedure  for  Simulation 
Experimentation",  Proceedings  of  the  1974  Winter  Simulation  Conference, 
Washington,  D.C.,  1974,  p.  491-497. 

4.  Biles,  W.  E. ,  and  J.  J.  Swain,  "Strategies  for  Optimization  of  Multiple  Res¬ 
ponse  Simulation  Models",  Proceedings  of  the  1977  Winter  Simulation  Conference, 
Gaithersburg,  Maryland,  December  1977. 

5.  Biles,  W.  E. ,  and  J.  J.  Swain,  "Mathematical  Programming  and  the  Optimization 
of  Computer  Simulations",  Mathematical  Programming  Study  in  Engineering 
Optimization,  North-Holland  Publishing  Company,  Amsterdam,  Netherlands,  1979. 

6.  Biles,  W.  E.  and  J.  J.  Swain,  Optimization  and  Industrial  Experimentation, 

V'iley- Interscience,  New  York,  1979. 

7.  Box,  G.  E.  P.,  and  K.  B.  Wilson,  "On  the  Experimental  Attainment  of  Optimum 
Conditions",  Journal  of  the  Royal  Statistical  Society,  B,  Vol.  13,  1951,  p. 

1-45. 

8.  Box,  G.  E.  P. ,  and  J.  S.  Hunter,  "Multifactor  Experimental  Designs  for 
Exploring  Response  Surfaces",  Annals  of  Mathemtical  Statistics,  Vol.  28,  1957,  p. 
242. 

9.  Box,  J.  J. ,  "A  New  Method  of  Constrained  Optimization  and  Comparison  with 
other  Methods",  Computer  Journal,  Vol.  8,  1965. 

10.  Brooks,  S.  H. ,  and  M.  R.  Mickey,  "Optimization  Estimation  of  Gradient  Direction 
in  Steepest  Ascent  Experiments",  Biometrics,  Vol.  17,  No.  1,  1961. 

11.  Carter,  G.  and  E.  J.  Ignall,  "A  Variance  Reduction  Technique  for  Simulation", 
Management  Science,  Vol.  21,  1975,  pp.  607-616. 

12.  Conway,  R. ,  "Some  Tactical  Problems  in  Digital  Simulation",  Management  Science, 
Vol.  10,  1963,  pp.  47-61 

13.  Crane,  M.  A.  and  D.  L.  Iglehart,  "Simulating  Stable  Stochastic  Systems  I: 

General  Multiserver  Queues",  J.  ACM,  Vol.  21,  1974,  pp.  103-113. 


14.  Crane,  M.  A.  and  A.  Lemoine,  An  Introduction  to  the  Regenerative  Method  for 
Simulation  Analysis.  Technical  Report  No.  86-23,  California  Analysis  Corp¬ 
oration,  Palo  Alto,  CA,  October  1976. 

15.  Draper,  N.  R. ,  and  H.  Smith,  Applied  Regression  Analysis,  John  Wiley  and  Sons, 


New  York,  1966. 


30 


16.  Dudewicz,  E.  J.,  "New  Procedures  for  Selection  Among  (Simulated)  Alternatives", 
Proceedings .  1977  Winter  Simulation  Conference,  1977,  pp.  58-62. 

17.  Duket,  S.  and  A.  A.  B.  Pritsker,  "Examination  of  Simulation  Output  Using 
Spectral  Methods",  Mathematics  and  Computers  in  Simulation,  Vol.  XX,  1978.  pp. 
53-60 

18.  Farrell,  W. ,  "Literature  Review  and  Bibliography  of  Simulation  Optimization", 
Proceedings,  1977  Winter  Simulation  Conference,  1977,  pp.  116-124. 

19.  Fishman,  G.  S.,  "A  Study  of  Bias  Considerations  in  Simulation  Experiments", 
Operations  Research,  Vol.  20,  1972,  pp.  785-790. 

20.  Fishman,  G.  S. ,  Concepts  and  Methods  in  Discrete  Event  Digital  Simulation,  John 

Wiley,  1973.  .  ~  ~  "  '  ' 

21.  Fishman,  G.  S.,  "Grouping  Observations  in  Digital  Simulation",  Management 
Science,  Vol.  24,  1978,  pp.  510-521. 

22.  Fishman,  G.  S. ,  Principles  of  Discrete  Event  Simulation,  John  Wiley,  1978. 

23.  Fishman,  G.  S. ,  "Statistical  Analysis  for  Queueing  Simulation",  Management 
Science,  Vol.  20,  1973,  pp.  363-369. 

24.  Fishman,  G.  S. ,  "The  Allocation  of  Computer  Time  in  Company  Simulation  Experir- 
ments",  Operations  Research,  Vol.  16,  1968,  pp.  280-295/ 

25.  Hunter,  J.  S.,  and  T.  H.  Naylor,  "Experimental  Designs  for  Computer  Simulation 
Experiments",  Management  Science,  Vol.  16,  No. 5  ,  1970,  pp.  422-434. 

26.  Ignall,  E.  J.,  "On  Experimental  Designs  for  Computer  Simulation  Experiments", 
Management  Science,  Vol.  18,  No.  7,  March  1972,  pp.  384-388. 


27.  Kabak,  I.  W. ,  "Stopping  Rules  for  Queueing  Simulations",  Operations  Research. 
Vol.  16,  1968,  pp.  431-437. 

28.  Kaplan,  D.  J. ,  L.  C.  Davis,  M.  E.  B.  Owens  II,  0.  F.  Forsyth  and  J.  J.  Penick, 
"SPEARS,  an  AAW  Performance  Simulation",  NRL  Report  7958,  Naval  Research 
Laboratory,  Washington,  D.  C. ,  June  1976. 

29.  Kleijnen,  J.  P.  C. ,  Statistical  Techniques  in  Simulation,  Parts  I  and  II,  Marcel 
Dekker,  New  York,  1975. 

Law,  A.  M. ,  Confidence  Intervals  for  Steady-State  Simulations,  1:  A  Survey  of 
Fixed  Sample  Size  Procedures.  Technical  Report  78-5,  University  of  Wisconsin, 
1978. 


30. 


31 


31.  Law,  A.  M. ,  "A  Comparison  of  Two  Techniques  for  Determining  the  Accuracy 
of  Simulation  Output1',  Research  Report  (Department  of  Industrial 
Engineering,  University  of  Wisconsin  -  Madison),  June,  1975. 

32.  Law,  A.  M. ,  "Statistical  Analysis  of  the  Output  Data  from  Terminating 
Simulations",  Research  Report  (  Department  of  Industrial  Engineering, 

University  of  Wisconsin  -  Madison),  April,  1978. 

33.  Law,  A.  M. ,  Carson,  J.  S.,  "A  Sequential  Procedure  for  Determining  the 
Length  of  a  Steady-State  Simulation",  Research  Report  (Dept,  of  Ind.  Engr. , 
University  of  Wisconsin  -  Madison),  April,  1977. 

34.  Law,  A.  M. ,  "Validation  of  Simulation  Models,  I:  An  Overview  and  Survey 

of  Real-World  Practice",  Research  Report  (Department  of  Industrial  Engineering, 
University  of  WIsconsin-Madison) ,  August,  1979. 

35.  Montgomery,  D.  C. ,  Design  and  Analysis  of  Experiments,  John  Wiley  and  Sons, 

New  York,  1976. 

36.  Montgomery,  D.  C. ,  and  V.  M.  Bettencourt,  "Multiple  Response  Surface  Methods 
in  Computer  Simulation",  Simulation,  October,  1977,  pp.  113-121. 

37.  Montgomery,  D.  C. ,  "Methods  for  Factor  Screening  in  Computer  Simulation 
Experiments,"  Research  Report  on  ONR  Contract  N0014-78-C-0312,  Georgia 
Institute  of  Technology,  Atlanta,  GA,  March  1975. 

38.  Pritsker,  A.  A.  B. ,  "GASP",  Enclopedia  of  Computer  Science  and  Technology, 

J.  Belzer,  A.  G.  Holzman,  and  A.  Kent,  Editors,  Vol.  8,  Marcel  Dekker,  Inc., 
1977. 

39.  Pritsker,  A.  A.  B. ,  Modeling  and  Analysis  Using  Q-GERT  Networks,  Halsted 
Press  and  Pritsker  &  Associates,  Inc.,  1977. 

40.  Pritsker,  A.  A.  B. ,  The  GASP  IV  Simulation  Language,  John  Wiley,  1974. 

41.  Pritsker,  A.  A.  B. ,  and  C.  D.  Fegden,  Introduction  to  Simulation  and  SLAM, 
Halsted  Press,  New  York,  1979. 

42.  Russell,  E.  C. ,  Simulating  with  Processes  and  Resources  in  SIMSCRIPT  II.  5, 
CACI,  Inc.,  Arlington,  VA. ,  1974. 

43.  Sabuda,  J.,  F.  H.  Grant,  III  and  A.  A.  B.  Pritsker,  The  GASP  IV/E  User's 
Manual ,  Pritsker  &  Associates,  Inc.,  West  Lafayette,  EN,  1978. 

44.  Schmeiser,  B.  W. ,  Babu,  A.  J.  G. ,  "Beta  Variate  Generation  via  Exponential 
Majorizing  Functions"  Research  Report  (Department  of  Operations  Research 
and  Engineering  Management,  Southern  Methodist  University),  December,  1978. 

45.  Schmeiser,  B.  W. ,  Lai,  R. ,  "Another  Versatile  Family  of  Probability 
Distributions"  Research  Report  (Department  of  Operations  Research  and 
Engineering  Management,  Southern  Methodist  University),  December,  1978. 

46.  Schmeiser,  B.  W. ,  "Methods  for  Modeling  and  Generating  Probabilistic  Com¬ 
ponents  in  Digital  Computer  Simulation  when  the  Standard  Distributions  are 
not  Adequate:  A  Survey"  (Department  of  Operations  Research  and  Engineering 
Management,  Southern  Methodist  University),  December  1977. 


32 


47.  Schmeiser,  B.  W. ,  Lai,  R. ,  "Computer  Generation  of  Bivariate  Gamma 
Random  Vectors"  Research  Report,  (Department  of  Operations  Research, 

Southern  Methodist  University),  1979. 

48.  Schruben,  L.  W. ,  "A  Coverage  Function  for  Interval  Estimates,  "Research 
Report  (School  of  Operations  Research  and  Industrial  Engineering,  Cornell 
University),  February  1979. 

49.  Schruben,  L.  W. ,  "Controlling  Initialization  Bias  in  Simulation  Experi¬ 
ments",  Research  Report  (School  of  Operations  Research  and  Industrial 
Engineering,  Cornell  University),  April,  1979. 

50.  Schruben,  L.  W. ,  Margolin,  B.  H. ,  "Psendo-random  Number  Assignment  in 
Statistically  Designed  Simulation  and  Distribution  Sampling  Experiments", 
Journal  of  the  American  Statistical  Association,  Vol.  73,  No.  363,  1978, 
pp.  504-520. 

51.  Smith,  D.  E.,  "Requirements  of  an  Optimizer  for  Computer  Simulations", 

Naval  Research  Logistics  Quarterly,  Vol.  20,  No.  1,  1973. 

52.  Smith,  D.  E. ,  "An  Empriical  Investigation  of  Optimum-Seeking  in  the 
Computer  Simulation  Situation",  Operations  Research,  Vol.  21,  1973,  pp. 
475-497. 

53.  Smith,  D.  E. ,  "Automated  Response  Surface  Methodology  in  Digital  Computer 
Simulation,"  Volumes  I  and  II,  Research  Report.,  ONR  Contract  N00014- 
74-C-0148,  State  College,  PA,  September  1975. 

54.  Welsh,  J.  and  J.  Elder,  Introduction  to  PASCAL.  Prentice-Hall,  International, 
London,  1979. 


( 

I 

I 


APPENDIX  A 

DESCRIPTION  OF  DATA  INPUT  FORMAT 
FOR  NAVAL  SYSTEMS  SIMULATION  PROGRAM 
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Input  Cards 


1.  Project  Card 

Title  of  the  project,  name  of  the  analyst,  and  the  date  to  appear  on 
the  summary  report  are  specified  on  this  card. 


Cols 

Description 

Format 

1-4 

PROJ 

2A2 

10  -  39 

Title  of  the  project 

15A2 

40  -  59 

Analyst  name 

10A2 

60  -  69 

Date 

5A2 

Discrete 

Card 

Number  of 

the  required  files  in  the  model,  whether  a 

discrete  model  or 

not,  and 
card. 

maximum  cumber  of  attributes  per  entity  are 

specified  on 

this 

Cols 

Description 

Format 

1-4 

DISC 

10  -  11 

20 

Number  of  files 
j  1  we  have  discrete 

0  we  do  not  have  discrete 

12 

11 

30  -  31 

Maximum  number  of  attributes 
per  entity 

12 

Rank  Card 

The  ranking  discipline  of  each  of  the  files  are  specified  on  this 

card 

For  every 

file  there  should  be  one  card. 

Cols 

Description 

Format 

1-4 

RANK 

10  -  11 

File  number 

12 

20 

Ranking  discipline 

11 

1  FIFO 

2  LIFO 

3  HVF 

4  LVF 

30  -  31 

Number  of  ranking  attribute 

12 

Example : 

1.  To  specify  that  file  5  is  ranked  FIFO  we  have: 

1  10  20 

111 
RANK  05  1 
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2.  To  specify  that  file  11  is  ranked  HVF  (High  Value  First)  based  on 
attribute  6  we  have: 


Cols 

] 

10 

20 

30 

1 

4- 

RANK 

11 

3 

06 

4 .  Continuous  Card 

The  information  on  the  continuous  part  of  the  model,  if  any,  is  provided 
through  this  card.  If  the  model  does  not  have  a  continuous  part,  this 
card  must  still  be  used  but  all  the  parameters  are  left  blank. 


Cols 

Description 

Format 

1-4 

CONT 

10  -  11 

No.  of  Differential  Equations 
(NEQD) 

12 

15  -  16 

No.  of  Difference  Equations  (NEQS) 

12 

20  -  21 

No.  of  state  Events  (NSEV) 

12 

30  -  39 

30  -  39 

Minimum  step  size  (DTMIN) 

f:  o .  5 

40  -  49 

"Maximum  step  size  (DTMAX) 

F10.5 

50  -  59 

Absolute  Error  (AERR) 

F10.5 

60  -  69 

Relative  Error  (RERR) 

F10.5 

70 

Indicates  type  of  error  check 

Af 

in  Runge-Kutta  integration  or  in 
state  event  crossing  detection 


when  a  step  size  smaller  than 
DTMIN  is  required.  If  F  is  spec¬ 
ified,  a  fatal  error  occurs.  If 
W  is  specified,  a  warning  message 
is  printed  before  proceeding.  If 
N  is  specified,  execution  proceeds 
with  no  warning  message  given. 

The  default  value,  that  is  if  the 
field  is  left  blank,  is  W. 


The  numerical  integration  accuracy  is  controlled  by  the  specification 
of  AERR  and  RERR.  The  Runge-Kutta-Fehlberg  integration  algorithm  used 
in  the  language  estimates  the  single  step  error  for  each  variable  de¬ 
fined  by  a  differential  equation.  The  Ith  error  estimate  is  compared 
to  TERR  where 


TERR(I)  =  AERR  +  ABS  (SS(I>)  x  RERR. 
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If  the  error  estimate  is  less  than  or  equal  to  TERR(I)  for  each  I,  the 
values  of  SS(I)  are  accepted.  If  not,  the  step  size  is  reduced  and  the 
integration  algorithm  is  reapplied.  There  are  no  default  values  for 
AERR  and  RERR  and  they  must  be  specified  by  the  user.  The  stringent 
values  for  these  could  substantially  increase  the  running  times,  although 
they  result  in  better  accuracy. 

5.  State-event  Card 


The  information  on  state  events  is  provided  on  this  card. 


Cols 

Description 

Format 

1-4 

SEVN 

10  -  11 

Event  number 

12 

20  -  23 

Crossing  variables 

12 

30  -  31 

+i  SS(i) 

-i  DD(i) 

Direction  of  crossing 

12 

40  -  43 

+1  positive  direction 

0  either  direction 
-1  negative  direction 
Crossed  variable 

13 

+1  SS(i) 

-1  DD(i) 

0  a  constant  value 

50  -  59'. 

Constant  value;  if  zero  is 

F10.5 

60  -  69 

specified  in  the  last  field 
Tolerance  of  crossing 

F10.5 

Example 

Below  are  3  examples  of  the  SEVN  cards: 

1.  Define  state-event  1  to  occur  when  SS(3)  crosses  100  in  the  positive 
direction  with  a  tolerance  of  2. 


Cols  1 

10 

20 

30 

40 

50 

60 

4 

4 

4 

4 

4 

4 

4 

SEVN 

01 

+03 

+1 

100. 

2. 

Define  state-event  2 

to  occur 

when 

SS(2)  crosses  DD(1) 

in  the 

tive  direction 

with  a 

tolerance  of 

0.01. 

Cols  1 

10 

20 

30 

40 

50 

60 

4 

4 

4 

4 

4 

4 

4 

SEVN 

02 

+02 

-1 

-01 

.01 

positive  crossing 
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6.  Statistics  Card 


In  this  card  number  of  tallies,  discrete  time  persistent  statistics, 
continuous  time  persistent  statistics,  counters  are  specified.  Also, 
DTSAV,  time  between  saving  continuous  variables  is  declared  on  this 
card. 


Cols 

Description 

Format 

1-4 

STAT 

10  -  11 

No.  of  tallies 

12 

20  -  21 

No.  of  discrete  time 
persistent  statistics 

12 

30  -  31 

No.  continuous  time 
persistent  statistics 

12 

40  -  41 

No.  of  counters 

12 

50  -  59 

The  time  between  observing 
•  the  continuous  time  persistent 
statistics  (DTSAV).  If  it  is 
negative,  the  statistics  are 
collected  every  ABS (DTSAV)  plus 
before  and  after  processing  any 
event.  In  this  situation,  after 

F10.5 

Tally  Card 

processing  the  event  only  the  value 
of  those  continuous  variables  that 
are  changed  within  the  event  are 
collected. 

For  each  tally  one  card  must  be  prepared.  Tally  number  and  name  are 
specified  on  this  card.  For  future  use,  such  as  drawing  histograms, 
the  analyst  might  want  to  have  access  to  every  individual  observation 
By  specifying  an  output  unit  number  in  this  card,  each  individual 
observation  will  be  written  unformatted  on  this  unit.  Care  must  be 


practiced 

to  provide  the  proper  JCL  of 

the  computer  installation. 

Cols 

Description 

Format 

1-4 

TALY 

10  -  11 

Tally  number 

12 

20  -  31 

Tally  name 

6A2 

40  -  41 

Unit  number 

12 

Discrete 

Time  Persistent  Card 

For  each  discrete  time  persistent  statistic  one  card  must  be  prepared. 
Number  and  name  of  the  discrete  time  persistent  are  specified  on  this 
card.  Also,  the  output  unit  number  for  writing  off  each  individual 
observation,  like  the  tally  card,  can  be  specified  on  this  card.  There 
are  two  types  of  discrete  time  persistent  statistics  that  can  be  col¬ 
lected.  One  is  on  the  number  of  the  entries  of  a f ile ,  and  the  other  is 
on  one  of  the  XX  variables. 


39 


Cols 

Description 

Format 

1-4 

DTPS 

10  -  11 

Number 

12 

20  -  31 

Name 

6A2 

40  -  42 

-i  ith  XX 
+i  Number 

variable 
of  entries  of 

13 

file  i 

50  -  51 

Unit  number 

12 

Examp le 

Suppose  that  discrete  time  persistent  number  3  is  for  collecting  statis¬ 
tics  on  the  variation  of  number  of  the  entries  on  file  2  and  it  is  to 


be  named  "NO. 

IN  QUE.” 

.  Plus,  it  is 

desired 

to  "dump" 

out  all  the 

changes  of  the 

number 

of  the  entries 

of  file 

2  and  the 

times  of  changes 

on  output  Unit 

8.  The 

following  card  provides  this  information. 

Cols  1 

10 

20 

40 

50 

4* 

+ 

+ 

DTPS 

03 

NO.  IN  QUE. 

+02 

08 

9.  Continuous  Time  Persistent  Card 

For  each  continuous  time  persistent  statistic  one  card  must  be  provided. 
The  function  of  this  card  is  the  same  as  the  DTPS  card  except  this  is 
for  the  continuous  variables,  i.e.  DD(I)  or  SS(I). 


Cols 

Description 

Fon 

1-4 

CTPS 

10  -  11 

Number 

12 

20  -  31 

Name 

6A2 

40  -  42 

+i  SS (i) 

-i  DD(i) 

13 

50  -  51 

Unit  number 

12 

Another  difference  between  discrete  and  continuous  time  persistent 
statistics  is  that,  in  discrete,  all  the  changes  are  recorded,  but  for 
continuous,  every  DTSAV  the  values  are  recorded.  In  the  case  of  nega¬ 
tive  DTSAVi  the  values  of  continuous  variables  are  recorded  every 
ABS(DTSAV),  before  every  event  and  after  the  event.  The  latter  if  the 
values  have  changed  during  the  event. 

10.  Counter  Card 


This  is  for  defining  the  counters.  For  each  counter  one  card  must  be 
prepared  with  specifying  number,  name,  and  upper  limit  of  the  counter. 

If  at  the  end  of  an  event  the  value  of  one  of  the  counters  is  larger  than 
its  corresponding  limit,  the  simulation  would  be  terminated.  Counters 
can  also  be  used  in  subroutine  STATE.  At  which  case  after  completion 
of  one  step,  the  value  of  one  of  the  counters  might  cause  the  simulation 
to  be  stopped. 
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Cols 

Description 

Format 

1-4 

COUN 

10  -  11 

Number 

12 

20  -  31 

Name 

6A2 

40  -  45 

Limit 

16 

11.  Stream  Cards 

This  is  for  specifying  the  Initializing  seeds  of  the  random  number 
generators.  There  are  a  maximum  of  10  streams  allowed.  The  language 
provides  the  default  values.  Therefore,  only  for  the  streams  whose 
initial  seed  is  to  be  different  from  the  default  values  cards  must 
be  prepared.  One  can  specify  antithetic  random  numbers  by  specifying 
a  negative  seed. 


Cols 

Description 

Format 

1-4 

STRM 

10  -  11 

Number 

12 

20  -  30 

Seed  value 

Ill 

Simulate 

Card 

Number  of 

runs  to  be  performed  is  specified  on  this  card. 

Cols 

Description 

Format 

1-4 

SIMU 

10  -  11 

Number  of  runs 

12 

13.  Initialization  Card 

This  is  for  indicating  the  starting  and  finishing  time  of  each  run  of 
the  simulation,  whether  or  not  to  initialize  the  seeds  between  the  runs, 
to  suppress  the  echo  report,  to  suppress  the  summary  report,  to  clear 
the  files  between  runs,  to  clear  variables  between  runs,  and  finally 
to  clear  the  statistics  and  if  so,  the  time  to  clear. 

Cols  Description  Format 


1-4  INTL 

10  -  19  Starting  time  of  the  simulation 

20  -  29  Finishing  time  of  the  simulation 

..  1  to  initialize  the  seeds  between  the  runs 

4() 

0  otherwise 

^  1  to  suppress  the  echo  report 

**"  0  otherwise 

1  to  suppress  the  summary  report 

0  otherwise 

^  1  do  not  clear  the  files  between  the  runs 

0  otherwise 


F10.3 

F10.3 


I 
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II 
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Cols 

Description 

Format 

60 

1 

do  not  clear  the  variables  between 

the  runs  ^ 

0 

otherwise 

65 

1 

do  not  clear  statistics 

11 

0 

otherwise 

70  -  79 

Time  to  clear  the  statistics 

F10.3 

APPENDIX  B 

DATA  INPUT  AND  SAMPLE  OUTPUT  FROM 
NAVAL  SYSTEMS  SIMULATION  PROGRAM 
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0004  CALL  SCHD  ( 3, T , ATB I B) 

0005  WRITE  (NPRTR,  101)  TNOW ,  (SS  (I)  ,  1=  1  r  4) 

0006  101  FORM  AT  (  1 0X ,  5  (  F 1 0.  5,  1  OX)  ) 

0007  RETURN 

0008  END 
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1. 

DIMENSION  RS ET  (5000) 

2. 

CCMMON/LOC/LFI,LLR,LFFB,LFTB,LFDSB,LFCSE,LFCTE, LFAE, LCAL 

3. 

CC'  ON  IS  ET  (5000) 

4. 

CCMMON/GSCI/NCHDR.NPRTR, SS  (99)  ,DD(99)  , ATRIB(99)  ,SSL(99)  ,JJ, 

5. 

-DDL (99) ,TNOH,  XX (99 ) , DTNON, ISTCF 

6. 

EQUIVALENCE  (IS ET { 1 ) , RS ET  (  1) ) 

7. 

NCHDR=5 

8. 

NFRTR=6 

9. 

IFI=  1 

10. 

LLR=500 1 

11. 

CALL  SOAP 

12. 

STOP 

13. 

END 

14. 

SUBROUTINE  SOAP 

15. 

DIMENSION  BS ET  {  1 ) 

16. 

CCMMON/LOC/LPI ,LLR,LFFB, LFTB,  LFDSB,LFCSB,LFCTB,LFAE,LCAL 

17. 

COM HON/GEN 1/NRUNS,NFILE,TBEG,TFIN,ICLF,ICLV, ICLS,TCL EAR, IRAN, 

18. 

-NTSX, LTHX (99)  , NTAL, NDTPST, NCTPST, NCOUN , N ATR , NENT, I  ECO, ISOM 

19. 

COMMON  ISET(1) 

20. 

EQUIVALENCE  (ISET ( 1) , HSET (  1) ) 

21. 

DATA  STARS/2H**/ 

22. 

CALL  DEFALT 

23. 

CALL  DATAIN 

24. 

IF  (IECO.  EQ.  1)  GO  TO  10 

25. 

CALL  ECHO 

26. 

10 

DO  20  1=1 , NSUNS 

27. 

IF (ICLF. EQ. 1)  GO  TO  11 

28. 

CALL  CLEARF 

29. 

1 1 

IF (ICLV. EQ. 1)  GO  TO  12 

30. 

CALL  CLEARV 

31. 

12 

CALL  INTAL 

32. 

CALL  EXEC 

33. 

IP (ISUM. EQ. 1)  GO  TO  13 

34. 

CALL  SUMRY  (I) 

35. 

13 

IF (NDTPST. EQ.O)  GO  TO  15 

36. 

DO  14  J=1, NDTPST 

37. 

L=LFDSB+  (J-1)  *9 

38. 

NU=ISET  (L+8) 

39. 

IF(NO.LE.O)  GO  TO  14 

40. 

iFITE(NU)  STARS 

41. 

14 

CONTINUE 

42. 

15 

IF( NCTPST. EQ.O)  GO  TO  19 

43. 

DO  16  J=1, NCTPST 

44. 

L  =  LFCSB+  (J-1)  *9 

45. 

NU=ISET  (L+8) 

46. 

IF(NU.LE.O)  GO  TO  16 

47. 

NRITE(NU)  STABS 

48. 

16 

CONTINUE 

49. 

19 

CALL  OOTPUT 

50. 

20 

CONTINUE 

51. 

STOP 

52. 

END 

53. 

SUBROUTINE  DEFALT 

54. 

DIMENSION  RSET { 1 ) , LL ( 10) 

55. 

DATA  LL/1274321477,  2135124613,  174325 154 1 ,  162421 7675 , 20 14632579 

A  SIMULATION  LANGUAGE 


71 


1  2  3  4  5  6 

1 2 34567890 1234567890 12 345 6789 0123456789012345678901 234 5678 9012 3456789 


56. 

57. 

58. 

59. 

60. 
61. 
62. 

63. 

64. 

65. 

66. 

67. 

68. 

69. 

70. 

71. 

72. 

73. 

74. 

75. 

76. 

77. 

78. 

79. 

80. 
81. 
82. 

83. 

84. 

85. 

86. 
8/. 
98. 

89. 

90. 

91. 

92. 

93. 

94. 

95. 

96. 

97. 

98. 

99. 
100. 
101. 
102. 

103. 

104. 

105. 

106. 

107. 

108. 

109. 

110. 


-20  36774231  ,  1 4 5231  3  57 1 ,  1 2  542406 57 , 1  4  10 1 4 336 3 , 2 1  356 21  895/ 
CCMMON/LOC/LFI, LLR, LFFB , LFTE , LF ES B, LFCS E , LFCTB, LFAE , LC AL 
CCMMON/GEN1/NRONS,N  FILE, TBEG, TF IN , ICLF, ICLV, ICLS , TCLEAR , IRAN , 
-NTMX,LTMX (99) , NT AL, N DTPS! , NCTPST, NCCUN, NAT R, N ENT, I ECO, ISOM 
COMMON/GEN4/IS  EED ( 1 0) ,LSEED  (10) 

CCMMON/GEN8/NIR , VALH, VALI 
COMMON  ISET(1) 

EQUIVALENCE  (ISET  ( 1 ) , BSET ( 1) ) 

C* 

C*  THIS  SOBROOTINE  ASSIGNS  THE  DEFAULT  VALUES 
C*  TO  THE  SOAP  VARIABLES. 

C* 

N  IR=  1 

VALH=1. E20 

VALL=-1.E20 

ICLF=0 

ICLV=0 

ICLS=0 

CALL  CLEARV 

NFILE=0 

LFFB=0 

LFTB=0 

LFDSB=0 

LFCSB=0 

LFAE=0 

LCAL=0 

NTMX=0 

NFNT=0 

N ATR=0 

NT AL=0 

NDTPST=0 

NCTPST=0 

N=LLR-1 

DO  1  1=1, N 

ISET  ( I)  =0 

1  CONTINUE 

DO  2  1=1,  10 

ISEED  (I)  =LL  ( I) 

LSEED(I)  =LL  (I) 

2  CCNTINUE 
RETURN 
END 

SUBROUTINE  ECHO 
DIMENSION  RSET  (1)  ,  INFO  (3) 

DATA  INFO/1HN, 1HN,1HF/ 

CGHHON/LOC7LFI , LLR, LFFB , LFTB,LFDSB, LFCSB, LFCTB , LFAE , LC AL 
COMMON/GEN 1/NRUNS.N FILE, TBEG, TFIN, ICLF, ICLV, ICLS, TCLEAR, IRAN, 
-NTMX,LTMX (99) , NT AL, NDTPST, NCTPST, NCOUN , NATR, NENT , I  ECO, ISUM 
CCNMON/GEN2/NEQD,NEQS, NEQT, DTSAV, DTMI N, DTMAX , AERB, RERR, I  ERR, NS  El 
-ISEES,ICSV  (25)  ,  ICGV  (25)  ,ICDIR(25)  ,ICDV(25)  ,VAL(25)  ,TTOL(25)  , 
-ISCD (25) 

CCMHON/GEN4/ISEED ( 1 0) ,LSEED(10) 

COMMON/GSC 1/NCRDR, NPRTB , SS  (99) ,DD(99) ,ATRIB(99)  ,SSL(99)  ,JJ, 

-DDL (99)  ,TNOW,  XX  (99)  ,DTNCV,ISTCP 
CCMMON  ISET  ( 1) 
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111. 

EQUIV ALENCE {ISET  (1)  ,HSET(1)) 

112. 

WRITE (NPRTR ,95) 

1  13. 

95 

FORMAT (  181 , 52X , 29HT  HE  ECHO  REPORT) 

114. 

IF  (NENT. GT.  0. A ND. N EQT. GT. 0)  GO  TO  11 

115. 

IF  (NENT.GT.O)  GO  TO  12 

1  16. 

WRITE (NPRTR, 100)  NEQT 

117. 

100 

FORMAT (// 10X,31HTHIS  IS  A  CONTINUOUS  MODEL  WITH  ,14, 

1  18. 

-21H  CONTINUOUS  EQUATIONS) 

119. 

GC  TO  20 

120. 

1 1 

WRITE (NPRTR,  101)  NEQT 

121. 

101 

F CRM AT (// 10X,29HTHIS  IS  A  COMBINED  MODEL  WITH, 14, 

122. 

-21H  CONTINUOUS  EQUATIONS) 

123. 

GC  TO  13 

124. 

12 

WRITE (NPRTR, 102) 

125. 

102 

FORMAT (//10 X, 24HTHI S  IS  A  DISCRETE  MODEL) 

126. 

13 

WRITE(NPRTR, 103) 

127. 

103 

FORMAT (//I OX, 16HFILE  I NFOR BATIC N/ 10 X , 16 ( 1 H-) ) 

128. 

IF  (NFILE. GT. 0)  GO  TO  16 

129. 

WRITE (NPRTR, 110) 

130. 

1 10 

PORMAT (10X,23HNO  FILES  ARE  DESIGNATED) 

131. 

GO  TO  15 

132. 

16 

DO  10  1=1, NFILE 

133. 

WRITE (NPRTR, 104)  I 

134. 

104 

FORMAT (//I 0X,9HFILE  NO.:  ,14) 

135. 

L=LFFB>  (1-1)  *6 

136. 

WHITE  (NPRTR  ,  105)  L 

137. 

105 

FORMAT ( 101, 23HLOCATION  OF  FILE  BLOCK:, 16) 

138. 

J=L*3 

139. 

WRITE  (NPRTR,  106)  ISET(J) 

140. 

106 

FORMAT (10X,22HRANKING  DISCIPLINE  IS:, 12/ 

141. 

- , 20X , 32H1 :  FCFS  2:  LCFS  3:  HVF  4:  LVF) 

142. 

J=L*4 

143. 

IF  (ISET  (J)  -EQ.O)  GO  TO  14 

144. 

WRITE (NPRTR, 107)  IS ET ( J) 

145. 

107 

FORMAT (10X,21HRANKING  ATRIBUTE  IS  :,I3) 

146. 

14 

J=L*5 

147. 

IF  (ISET  (J)  .EQ.O)  GOTO  15 

148. 

WRITE  (NPRTR,  108)  ISET(J) 

149. 

108 

FORMAT ( 10  X , 29HLOCAT 10  N  OF  STATISTICS  BLOCK:, 16) 

150. 

10 

CONTINUE 

151. 

15 

WRITE (NPRTR, 109)  LCAL 

152. 

1C9 

FGHNAT (//10X,28H****LOCATION  OF  CALENDER  IS:, 16) 

153. 

C* 

154. 

c* 

PRINT  CONTINUOUS  INFORMATION 

155. 

c* 

156. 

20 

IF  (NEQT. LE.O)  GO  TO  21 

157. 

WBITE(NPRTR, 120)  NEQD, NEQS 

158. 

120 

FORMAT (///10X, 22HCONTINUOUS  INFORMATION/ 10X, 22 (1H-) // 

159. 

- 1CX, 37HNO.  OF  DIFFERENTIAL  EQUATIONS  (NEQD) :, 10X, 15/1  OX, 

160. 

-37HNO.  OF  DIFFERENCE  EQUATIONS  (NEQS):  ,10X,l5) 

161. 

WRITE (NPRTR, 121)  DTHIN, DTH AX , AERR , R ERR , DTSAV 

162. 

121 

FORMAT ( 101, 26HHINIHUM  STEP  SIZE  (DTH IN)  : , 1 6X , F 1 0. 5/1  OX , 

163. 

-26HMAXIHUH  STEP  SIZE  (DTBAX) : , 16X,F 10. 5,/10X, 

164. 

-28HABSOLUTE  ERROR  LIMIT  (A ERR)  : ,  1 4X , FI  0. 5/ 1 0X , 

165. 

-28HRELATIVE  ERROR  LIMIT  (RER R) : , 1 4X , FI 0 . 5/10 X , 
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166. 

167. 

168. 

169. 

170. 

122 

171. 

172. 

173. 

130 

174. 

175. 

176. 

177. 

178. 

131 

179. 

180. 
181  . 

132 

182. 

183. 

23 

184. 

185. 

133 

186. 

24 

187. 

134 

188. 

189. 

190. 

191. 

135 

192. 

193. 

25 

194. 

195. 

136 

196. 

197. 

26 

198. 

137 

199. 

27 

200. 

138 

201. 

28 

202. 

21 

203. 

204. 

150 

205. 

206. 

207. 

208. 

151 

209. 

210. 

152 

211. 

212. 

153 

213. 

214. 

215. 

216. 

154 

217. 

218. 
219. 

155 

220. 

30 

-33HTIME  BETWEEN  SAVE  POINTS  (CTSAV)  : ,9X , FI  0. 5) 

IF  (IERR.EQ.-1)  WRITE  (NPRTR  ,  122)  INFO(1) 

IF (IERR.EQ.O)  WRITE (NPRTR, 122)  INFC(2) 

IF  (IERR.EQ. 1)  WRITE  (NPRTR,  122)  INFO  (3) 

FC3MAT(10X,36HACCURACY  ERROR  SPECIFICATION  ( IEBR)  :  1  5X  ,  A  1 ) 

IF  (NSEV.LE.  0)  GO  TO  21 
WRITE (NPRTR , 130) 

FORMAT (///10X, 12RSTATE  EVENTS/ 1 0X, 1 2  (  1 H-) //I  0 X, 68 NUMBER , 1  OX , 5  HE 
-T, 11X, 8 HC BOSSING, 10X,9HDIRECTION, 10X,7HCROSSED, 10X ,9  HTOL  ER  ANCE/ 
-, 8HVARI ABLE, 10X, 11HCF  C ROSS ING, 8X , 5 H VA LUE, 12X,11HOF  CROSSING//) 
DC  28  1=1 ,NSEV 
WRITE  (NPRTR,  131)  I,ICSV(I) 

FORM  AT ( 1 2  X , 1 2 , 1 3  X , I  3) 

IF  (ICGV(I)  .LT.O)  GO  TO  23 
WRITE (NPRTR, 132)  ICGV(I) 

FORMAT ( 1H+,42X,3HSS (,I2, 1H) ) 

GO  TO  24 
J=-ICGV  (I) 

WRITE (NPRTR,  133)  J 

FORMAT (1H*,42X,3HDD  (,I2,1H) ) 

WRITE (NPRTR ,134)  ICEIR(I) 

F0RMAT(1H+',64X,I2) 

IF  (ICDV  (I)  .  LT.  0)  GO  TO  25 
IF(ICDV(I)  .GT.  0)  GO  TO  26 
WBITE(NPRTR,135)  VAL(I) 

FORMAT (1H*,79X,F10. 5) 

GO  TO  27 
J=-ICDV (I) 

WBITE(NPRTR, 136)  J 

FORMAT (1H>,79X,3HDD  (,I2, 1H) ) 

GO  TO  27 

WRITE  (NPRTR,  137)  ICDV(I) 

FORMAT ( 1H*,79X,3HSS  (,I2,1H) ) 

WBITE(NPRTR,  138)  TTOL(I) 

FORMAT ( 1H*,96X,F10. 5) 

CONTINUE 

IF  (NTAL.EQ. 0)  GO  TO  31 
WRITE(NPRTB, 150) 

FORMAT (///10X, 17HTALLY  I NFOR M AT ION/ 10X , 1 7 ( 1 H-) ) 

DO  30  1=1 , NTAL 
L=LFTB*  (1-1)  *8 
WRITE (NPRTR, 151)  I 
FORMAT(//10X,9HTALLY  NO.,  13) 

WRITE(NPRTR,152)  L 

FORMAT { 10 X , 28 HLOCAT ION  OF  THE  TALLY  BLOCK:, 16) 

WRITE  (NPRTR,  153)  ISET(L) 

FCRMAT(10X, 30 HLOCAT ION  OF  STATISTICAL  ARRAY:, 16) 

J1=L*1 

J2=L*6 

WRITE (NPRTR , 154)  (ISET(J) ,J=J1,J2) 

FOR8AT(10X,7HLABEL:  ,6A2) 

J=L*7 

WRITE  (NPRTR,  155)  ISET(J) 

FORMAT( 10X, 18HOUTPUT  DEVICE  NO.:, 14) 

CONTINUE 
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i  21* 

31 

IF  (NDTPST.EQ. 0)  GO  TO  41 

I  22. 

WRITE (NPRTR, 160) 

23. 

160 

FORMAT (///10X,33HDISC.  TIME  PERSISTANT  I N FOR M ATIO 4/ 1  OX , 3 2 ( 1 H-) ) 

24. 

DC  40  1  =  1  , NDTPST 

i  25. 

L  =  LFDSBMI-1)  *9 

26. 

WRITE (NPRTR ,  161)  I 

27. 

161 

FORMAT (//I 0X,27HDISC.  TIME  PERSISTANT  NO.:  ,13) 

28. 

WRITE(NPRTR, 162)  L 

29. 

162 

FCRMAT(10X,1 8 H LOCATION  OF  BLOCK:, 16) 

30. 

WRITE  (NPRTR,  163)  ISET(L) 

:  31. 

163 

FCRMAT(10X,30HLOCATION  OF  STATISTICAL  ARRAY:  ,16) 

32. 

J  1  =  L+1 

33. 

J2=L+6 

34. 

WRITE (NPRTR,  154)  (ISET(J)  ,J=J1,J2) 

.35. 

J=L*7 

36. 

IF  (ISET (J) . GT. 0)  GO  TO  42 

37. 

J=-ISET  (J) 

38  - 

WRITE (NPRTR, 165)  J 

39. 

165 

FORMAT( 10X,45HTIME  PERSISTANT  STATISTICAL  COLLECTION 

ON  XX{ 

'40. 

-  ,  1 2,  1  H)  ) 

241. 

GO  TO  43 

42. 

42 

WRITE (NPRTR ,166)  ISET(J) 

4  3- 

16b 

FORMAT (10X,47HTIME  PERSISTANT  STATISTICAL  COLLECTION 

ON  QUEUE, 13) 

244. 

43 

J  =  L*8 

'45. 

WRITE  (NPRTR,  155)  ISET(J) 

.46. 

40 

CONTINUE 

2  47. 

41 

IF  (NCTPST.EQ.O)  GO  TO  51 

'43. 

WRITE(NPRTR, 170) 

'49. 

170 

FORMAT(///10X,33HCONT.  TIME  PERSISTANT  INFORM AT ION/IOX, 33 ( 1 H — ) ) 

.50. 

DC  50  1=1, NCTPST 

251 . 

L=LFCSB+(I-1)*9 

.52. 

WRITE (NPRTR, 17 1)  I 

53. 

171 

FORMAT (//10X,27H CO NT.  TIME  PERSISTANT  NO.:  ,13) 

>54. 

WRITE(NPRTR, 172)  L 

’55. 

172 

FCRMAT (10X, 18HLOCATION  OF  BLOCK:, 16) 

’56  . 

WRITE  (NPRTR,  173)  ISET(L) 

.57. 

173 

FCRMAT ( 10X, 30 H LOCATION  OF  STATISTICAL  ARRAY:  ,16) 

>58. 

J1=L«-1 

:  59  - 

J2=L+6 

:60. 

WRITE  (NPRTR, 154)  (ISET(J) ,J=J1,J2) 

261. 

J=L*7 

>62. 

IF  (ISET  (J)  .GT.O)  GO  ",0  52 

:63. 

J=-ISET (J) 

;64. 

WRITE  (NPRTR,  175)  J 

265. 

175 

FORMAT (10X,45HTIME  PERSISTANT  STATISTICAL  COLLECTION 

ON  DD( 

266. 

-,  12, 1H)  ) 

67. 

GO  TO  53 

268. 

52 

WRITE(NPRTR,176)  ISET(J) 

'69. 

176 

FOaHAT(10X,45HTlME  PERSISTANT  STATISTICAL  COLLECTION 

ON  £S  { 

270. 

-,I2,1H)> 

271. 

53 

J=L+8 

272. 

WBITE(NPBTR,  155)  ISET(J) 

273. 

50 

CONTINUE 

274. 

51 

IF (NCOUN. EQ.  0)  GO  TO  61 

275. 

WBITE (NPRTR , 180) 
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76. 

180 

77. 

73. 

79. 

|  30. 

181 

*31. 

82. 

182 

|93. 

1  34. 

35. 

.  36. 

I  37. 

1  83. 

183 

39. 

60 

1  90. 

61 

1  91. 

92. 

190 

93. 

94. 

191 

9  5. 

62 

96. 

140 

97. 

98. 

99. 

141 

00. 

70 

'01. 

(02. 

03. 

200 

i  04. 

05. 

80 

06. 

201 

i  07 ’ 

8  1 

1  38. 

202 

i09. 

10. 

111. 

12. 

C* 

13. 

C* 

14. 

C* 

15. 

16. 

'17. 

18. 

19. 

:20. 

21. 

22. 

23. 

'24. 

25. 

5 

>2  6. 

6 

i27  . 

28. 

29. 

:  30. 

C* 

FC  31  AT  (///  1  0  X  ,  19HCOUNTER  I  NFC  R  P  ATIC  N/  10  X  ,  1  9  (  1  H-)  ) 

DC  60  1=1 , NCOON 

L  =  LFCT3«- (1-1) *8 

WRITE (NPRTR , 181)  I 

FORMAT {//10X, 12HCOUNTER  NO.:, 13) 

WHITE (NPRTR , 182) L 

F CHS AT ( 10  X, 26 HLOCAT 10 N  CF  COUNTER  BLOCK:  ,16) 

J1=L+1 

J2=L+b 

WBITE(NPRTR, 154)  (ISET(J) ,J=J1,J2) 

J=  L+7 

WRITE  (NPRTR,  183)  ISET(J) 

FCR1AT  (  10X,7HLIMIT  :,I7) 

CONTINUE 

IF  (NENT.EQ.O)  GO  TO  62 
WRITE(NPRTR,  190)  LFAE 

F CBN AT (///I OX ,34HLO CATION  OP  FIRST  AVAILABLE  ENTRY:, 16) 
NFITE(NPRTR, 191)  NENT 

FORMAT (// 10X , 25H NO.  OP  ENTRIES  ALLOCATED: , 16) 

WHITE (NPRTR , 140) 

FC  RM AT  (// /10X,25HRANDQM  STREAM  INFORMAT ICN/1 OX , 25 { 1 H-) //) 

DC  70  1=1,10 

WRITE  (NPRTR, 141)  I,ISEED(I) 

FCHHAT ( 10X,5HSEED(, 12 , 2 H)  =  , 2X , 1 10) 

CONTINUE 

IF  (IRAN. EQ.  1)  GO  TO  80 
WFITE(NPRTR,200) 

FORMAT (/10X,42HSEEDS  WILL  NOT  BE  INITIALIZED  BETWEEN  RUNS) 

GO  TO  81 

WRITE (NPRTR, 201) 

FORMAT (/10X.38HSEEDS  WILL  BE  INITIALIZED  BETWEEN  RUNS) 

WRITE (NPRTR, 202)  NBUNS 

FORMAT (///10X , 8HTOTAL  OF,I4,23H  RUNS  WILL  BE  PERFORMED) 

RETURN 

END 

SUBROUTINE  CLEAFF 

THIS  ROUTINE  CLEARS  THE  FILES  BETWEEN  THE  RUNS 
DIMENSION  fiSET(l) 

COMMON/LOC/LFI, LLR, LFFB,LFTB, LFDSB, LFCS B ,LFCTB , LFAE , LCA L 
CCMMON/GEN1/NRUNS,NFILE,TPEG,TFIN,ICLF, ICLV,  ICLS,  TCLE AR,  IR AN, 
-  NT  MX, LT  MX (99) , NT  AL, N DT PST, NCT PST, NCOON, N  ATR, NENT, I ECO, ISDN 
COMMON  IS ET (1) 

EQUIVALENCE  (ISET ( 1 ) , RS FT ( 1 ) ) 

IF  (LCA L.  EQ.  0)  RETURN 

IF  (NFILE.EQ.O)  GO  TO  5 

L=LFFB 

GO  TO  6 

L=LCAL 

N=NFILE+1 

DC  10  K=1 , N 

I  =  L»(K-1)  *6 
ISET (I) =0 
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31. 

C* 

JCIN  ALL  THE  SNTBIES  TC  THE  FCCL  BY  TAKING  THE  LAST 

32. 

C* 

ENTRY  OF  THE  FILE  AND  SETTING  ITS  SUCCESSOR  AS  LFAE 

33. 

c* 

AND  UPDATING  LFAE 

34. 

c* 

1  35. 

11=14-1 

36. 

12=14-2 

37. 

LLFE=ISET  (11) 

i lg* 

LLLE= IS  ET (12) 

»  39. 

IF (LLFE. EO. 0)  GO  TO  10 

40. 

J  =  LLLE4- 1 

i  41* 

ISET (J) =LF A  E 

|  -*2. 

LFA  E=LLFE 

43  . 

ISET(II) =0 

44. 

ISET  (12)  =0 

1  45. 

10 

CONTINUE 

»  46. 

RETURN 

47. 

END 

48. 

S  C9R0UT IN  E  CLEAFV 

49. 

CC4i'10N/GSC1/NCRDR,NPBTR,SS  (99)  ,  DD  (99)  ,  ATRIP  (  9  9)  ,  S3  L  (99)  ,  J.J, 

50. 

-ICL  (99)  ,TNOW,  XX  (99)  , DTN0S,ISTCP 

51. 

DC  10  1=1,99 

52. 

SS  (I)  =0 

53. 

DD  ( I)  =0 

54. 

SSL  (I)  =0 

55. 

DDL  (I)=0 

1  56. 

XX  (I)  =0 

57. 

ATRIB  (I) =0 

58. 

10 

CONTINUE 

1  59. 

JJ=0 

'  60. 

RETURN 

61. 

END 

j  52. 

SUBROUTINE  INTAl 

I  o3. 

c* 

64  . 

c* 

THIS  IS  FOR  INITIALIZATION  BEFORE  EACH  RUN 

I  *5* 

c * 

•  66. 

DINSNSION  RSET  ( 1 ) 

67. 

CC  NflON/LOC/LFI,  LLR,  LF  FB,  LFT8  ,  LF  CSB,  L  FCSB  ,  LFCTB  ,  LFA  E  ,  LC  A  L 

68. 

CC  -INON/GE  N1/NRUNS,NFILE,TEEG,TFIN,  ICLF  ,ICLV,  ICLS,TCLEAB,  IRAN 

,  59. 

- NTM  X , LT HX (9  9) , NT  AL, NDTPST, NCTPST, NCOON, NATR, N ENT,  I ECO, IS UN 

'  70. 

CC'M10N/GEN4/ISEED(1Q)  ,LSEED(10) 

71. 

CC1.N0N/GSC1/NCRDR,NPRTR,  SS  (59)  ,  DD  (9  9)  ,  AT  PI  B  (99)  ,  SSL  (99)  ,  JJ, 

72  . 

-DDL  (99)  ,T  NOW  ,  XX(99)  ,DTNO»,ISTOP 

73. 

CCinON  ISET  ( 1 ) 

74. 

EQUIVALENCE  ( IS ET ( 1 )  , RS ET  ( 1 ) ) 

75. 

IF  (MCOUN.  EQ.O)  GO  TO  2 

76. 

DC  1  1=1, NCOUN 

77. 

LCTB=LFCTB4-  (1-1)  *8 

79. 

ISET(LCTB) =0 

79. 

1 

CONTINUE 

80. 

2 

ISTOP=0 

81. 

IF  (LCAL.EQ. 0)  GO  TO  3 

32. 

IF  (ICLS.EQ.  1)  GO  TO  3 

93. 

CALL  SCHD  (-1 ,TCLEAR, ATRIB) 

94. 

3 

IF  (IRAN. EQ.O)  RETURN 

95. 

DC  4  1= 1, 10 
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38b  . 

ISEED(I)  =  LS  EEC { I) 

387. 

4 

CONTINUE 

388  . 

RETURN 

3  39. 

END 

390  . 

SUBROUTINE  EXEC 

391  . 

DIMENSION  A2(100)  ,  A3J130)  ,  A4  (109)  ,A5(100) 

392. 

DIMENSION  RSET  ( 1) 

39  3. 

CCMMON/LOC/LFI, LLR, LFFP , LFT3, LFDSB, LFCSBrLFC 

i  394  . 

CCMMON/GEN1/NRUNS,  N FIL E,  T B EG , T F I N , I CLF , I CLV , 

395. 

-NTMX.LTHX (99) , NTAL, NDTPST, NCTPST,  NCOUN  ,  NATR  , 

396. 

CCMMO  N/GEN  2/NEQD ,N  EQS ,  NEQT,  DTSAV,  DT  *IN ,  DTM  AX 

13  97. 

-I  SEES,  ICS  V  (25)  ,  ICGV  (2  5)  ,ICDIR(2  5)  ,TCDV  (25)  ,V 

*398. 

-ISCD (25) 

399  . 

CCMMON/GEN8/NIR,  VALH,  VALL 

|*00. 

CCMMON/GSC1/NCRDR,NPRTR, SS  (99)  ,DD(99)  ,ATFIB( 

1-01. 

-D  CL ( 99)  ,T  NOW ,  XX  (99)  ,D1N0W,ISTCP 

402  . 

COMMON  IS  ET  (1) 

.  *03. 

EQUIVALENCE  (IS  ET  (1)  ,  RSET  (1)  ) 

I  404. 

WBITE(NPRTR,200) 

1  405. 

209 

FORMAT (//I  OX, 2 2H EXECUTION  IS  ATTEMPTED) 

4  06  . 

C* 

I  407. 

C* 

INITIALIZATION 

*  408. 

c* 

4  09. 

TTSAV=VALH 

.4  10. 

RESLS=. 0 1  *  ABS ( DTSAV) 

§411. 

RESLS  =  AfllNI  (DTHIN,  RESLS) 

*412. 

DTACC=DTN A  X 

413. 

TNOW=TBEG 

|4  14. 

FRMUL=. 5 

•  415. 

T  BUL=. 9 

4  16. 

ISEES=0 

|l  17. 

SCALL= 1.0 

|4ia. 

IF (8ERR.GT.0.)  SCA  LL=2. /REE  B 

*419. 

AEM=SCALL*AERR 

4  29  . 

CALL  INTLC 

1.21. 

IF  (NDT'  EQ.O)  GO  TO  10 

"•22. 

DC  5  I-  . NDTPST 

•  3. 

LDBL= LFDSB* (1-1 ) *9 

|- 2-4. 

CALL  DCLCT(LDBL) 

I  -5  . 

5 

CONTINUE 

1  0 

I F  (NEQT .EQ.O)  GO  TO  30 

T  L  AST=T  NOW 
DT  MOW  =  0 . 0 

IF (NSEV.GT.O)  CALL  SSCND 
~C  15  1  =  1  ,  N  EQT 
DHL  (I)  =DD  (I) 

3 ">  L  (I)  =  SS  (I) 

■<  '  '  11  f 
i  "ATE 

'  <  1.  E  .  .  9)  GOTO  16 

*  »  »  - 
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TB,LFAE,LCAL 
ICL3,TCLEARr IRAN, 
NENT,IECO,ISUM 
,  ASRR,  REP.R,  I  ERR,  NSEV, 
A L  (25)  ,  TTOL  (25)  , 


99) , SSL (99) ,JJ, 


•  f  :>TSA  ?) 


78 


I 

SIMULATION  LANGUAGE 

II  2  1  4  5  6  7 

1234567890  123456  78901 2  34567890  1234567890 1234567(390  1234  5678  99  12  34567890  123  4 


*41  . 

IF  (ISCD  (I)  .  EQ.O.OR.  ISCE  (T)  *ICDI!>  (I)  .  LT.O)  GO  TO  24 

442. 

IX=ICSV  (I) 

443  . 

CALL  SCHD  (IX,TNCW, ATRIB) 

*44. 

24 

ISCD(I) =0 

445. 

25 

CONTINUE 

4  46  . 

C* 

447. 

C* 

N  E  4  ITERATION 

448. 

c* 

449  . 

30 

T  N  3XT=V  ALH 

450. 

IF  (LCAL.EQ.O)  GO  TO  40 

451. 

IF  (ISET  (LCAL)  .E0.0)  GO  TC  35 

452. 

J=  LCAL*  1 

453. 

1 1 FE  =  ISET (J) 

454  . 

LATR=ISET ( LLFE+-3) 

455. 

TNEXT=RSET (LATH) 

456. 

35 

IF (NEQT.GT.O)  GO  TO  40 

457. 

IF  (TNEXT.GT.TFIN.  OR. ISTOP. EQ. 1 )  GO  TO  1000 

458. 

TNOW=TNEXT 

459. 

CALL  DISC 

460. 

GC  TO  30 

461. 

40 

IF(TNEXT.GT.TFIN.AND.TNCW.EC.TFIN)  GO  TO  1000 

462. 

IF  (TNEXT.GT.TFIN)  T  NEXT  =  TFIN 

463. 

45 

IF  (ISTOP.  EQ- 1)  GOTO  1000 

464. 

C* 

465. 

C* 

SAVE  CCNTINUOOS  VAR IABELS  IF  IT  IS  SAVING  TINE 

466. 

C* 

467. 

IF  (ABS(TTSAV-TNOW) -BESLS)  50,  50,55 

468. 

50 

CALL  SSAVE 

469. 

TTSAV=TTS AV+ABS (DTSAV) 

4  70. 

55 

IF  (TNEXT-TNOW)  60,65,70 

471. 

60 

WRITE (NPRTR, 960)  TNEXT, TNOW 

472. 

CALL  ERROR (9,1) 

473  . 

65 

IF  (DTSAV.  LT.O)  CALL  SSAVE 

474. 

CALL  DISC 

475. 

IF  (DTSAV.  GE.  0)  GO  TC  30 

476. 

.  c* 

477. 

c* 

SAVE  THE  VALUE  OF  THE  CTPST  VAR.  IF  THEY  HAVE  CHANGED 

478. 

c* 

IN  THE  EVENT 

479. 

c* 

480. 

DC  68  1=1 , NCTPST 

481. 

J=LPCSB*(I-1) *9 

482. 

LSTAT=ISET (J) 

483. 

ITYPE=ISET (J+7) 

484. 

IF  (ITYPE.  LT.  0)  GO  TO  66 

485. 

X=SS  (ITYPE) 

486. 

GC  TO  67 

487. 

66 

I1YPE=- ITYPE 

488. 

X=DD (ITYPE) 

489. 

67 

IF (X. EQ.RSET(LSTAT) )  GO  TO  68 

490. 

NU=I SET (J+8) 

491  . 

CALL  CTPST (X,LSTAT, NU) 

492. 

68 

CONTINUE 

493. 

GO  TO  30 

*94. 

70 

IF (NEQD.GT.O)  GO  TO  75 

4  95. 

DT?UL=DTMAX 

t 
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496. 

GC  TO  80 

497. 

75 

DTPUL=DTACC 

498. 

80 

IF  (TLAST+ DTFUL . GT.  TTS A V)  DTFUL=TTSA V-TLAST 

499. 

IF  (T L AST ♦  DTFUL.  GT.  TNEXT)  DT FfJ  L=  TN  EX T— TL AST 

500. 

IF  (NEQD.GT. 0)  GO  TO  240 

501. 

85 

DTMOW=DTFUL 

502. 

TNOW=TLAST+DTNOW 

503  . 

CALL  STATE 

504. 

90 

IF  (NSEV.GT. 0)  CALL  SSCND 

505. 

IF  ( IS  EES )  380,95,450 

506. 

C* 

1  507. 

c * 

NC  STATE  EVENT  HAS  CCCURED 

508. 

c* 

509. 

95 

TLAST=TNOW 

i  510. 

DC  100  1=1, NEQT 

1  511. 

DDL  (I)  =DD  (I) 

512. 

SSL  (I)  =  SS  (I) 

■  513. 

100 

CONTINUE 

1  514. 

IF  (TNOW.EQ.  TFIN)  GO  TO  30 

1  515. 

GO  TO  45 

516. 

C* 

1517. 

c* 

STATE  EVENT  HAS  PASSED;  REDUCE  THE  STEP  SIZE 

4  1  518. 

c* 

519. 

380 

IF  (DTFUL-DTHIN)  381,381,382 

.520. 

381 

IF  (IERR)  450,401,401 

I  521. 

4  Cl 

IF  (ISEES+1000)  410,400,400 

)  522. 

4  CO 

I=-TSEES 

523. 

WRITE (NPRTR, 980)  I,TNOW 

1  524. 

GC  TO  420 

■  525. 

410 

I=-ISEES- 1000 

i  526. 

WHITE (NPRTR, 990)  I,TNOW 

.  527. 

420 

IF  (IERR.EQ. 0)  GO  TO  450 

!  |528. 

CALL  ERPOR  (17,  1) 

1  529. 

382 

DTFUL=FRNUL* DTFUL 

530. 

IF  (DTFUL.  LT.  DTNIN)  DTFUL=DTHIN 

■  531. 

I£EES=0 

*  532. 

GO  TO  2  40 

533. 

C* 

1  534. 

C* 

AT  LEAST  ONE  STATE  EVENT  HAS  OCCUR  E  D ;  SCHEDULE  THEM  ON 

§  5  3  5  . 

c* 

THE  CALENDER  ;  ACCEPT  THE  SS  AND  DD  VALUES;  START  A  NEW  STEP 

536. 

c* 

1  .537. 

450 

ISESS=0 

1  538. 

DC  451  1=1, NSEV 

1  539. 

IF(ISCD(I)  .EQ.O.OE.ISCD  (I)  *ICDIE  (I)  .LT.O)  GO  TO  453 

I  540. 

IX=ICSV  (I) 

i  j  5  4  1 . 

CALL  SCHD  (IX,TNOW,  ATRIB) 

'  542. 

453 

ISCD(I)  =0 

1  543. 

451 

CCNTINUE 

■  544. 

TLAST=TNOW 

‘  c  / »  C 

DC  452  1=1, NEQT 

|  546. 

SSL (I) =SS  (I) 

|  547. 

DDL  (I)  =DD  (I) 

|  1  548* 

452 

CONTINUE 

•  ’  549. 

GC  TO  30  | 

|  550. 

;  1 

i 

i 

C* 

l 

i 

i- 
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SSI . 

C* 

RDNGE-KUTA  INTEGRATION 

SS2. 

C* 

553. 

240 

DTNOW=. 25*DTFUL 

554. 

DC  250  1=  1 , NEQD 

555. 

250 

S£  (I)  =  SSL  (I)  *DTNOW*DDL  (I) 

5  56. 

TNOW=TLAST*DTNOW 

557. 

CALL  STATE 

558. 

CH=3.0*DTFUL/32. 0 

559. 

DC  260  1=1, NEQD 

560  . 

A2  (I)  =DD  (I) 

561. 

260 

SS  (I)  SSL  (I)  *CH*  (DDL  (I)  +3.  0*A2  (I)  ) 

I  562. 

DTNOW=. 375*DTFUL 

563. 

TNOW=TL AST+DTNOW 

564  . 

CALL  STATE 

|  565. 

CH= DTFUL/ 2197.0 

f  5  6  6  . 

DO  270  1=1, NEQD 

567. 

A3  (I)  =DD  ( I) 

.  568. 

270 

SS  (I)  =SSL  (I)  *CH*  (1932.0*DDL  (I)  *  (729  6. 0*  A  3  (I)  -  7200.  0*  A2  (I)  )  ) 

I  569. 

DTNOW=12.  0*DTFUI/13.0 

1  570. 

TNOH=TL AST+DTNOH 

571. 

CALL  STATE 

I  572. 

CH=DTF  UL/4 104. 

1  573. 

DO  280  1=1, NEQD 

574. 

A 4  (I)  =DD  ( I) 

g  575. 

280 

SS  (I)  =SSL  (I)  *CH*  ((8  341.0*DDL(I)  -845.0*A4  (I)  )  +  (2944  0. 0* A 3  (I) 

1576. 

—  32832.  0*A2(I)  )  ) 

*577. 

DTNOW=DTPUL 

578. 

TN04=TL AST+DTNOR 

1579. 

CALL  STATE 

1580. 

CH=DTFUL/20520. 0 

531  . 

DC  290  1=1, NEQD 

■  582. 

A5  (I)  =DD  { I) 

|  593. 

290 

SS  (I)  =SSL  (I)  *CH*  ( (-6080.  *DDL  (I)  ♦  (9  29  5.  0*A4  (I)  -5643. 0*A5  (I)  )  1 

584. 

-♦  (41040. 0*A2 (I) -  28352.0* A3 (I)  ) ) 

585. 

DTNO«=. 5*DTFUL 

1  536  . 

TNOW=TLAST*DTNOW 

■  587. 

CALL  STATE 

588. 

CR=DTFUL/ 7618050.0 

■  539. 

E  EOET=0 . 0 

1590. 

DO  300  1=1, NEQD 

591. 

SS (I)  =SSL (I) +CH* ((902880. 0*rCL  (I) ♦ ( 38 55 735. 0*A 4  (I) -1371249.0* 

592. 

-A5  (I) ) )  *(39 53664. 0*A3 (I)  *277020. 0*DD(I)  )  ) 

(  593  . 

IF  (DTFUL.  LT.  DTHIN)  GO  TO  300 

»  594. 

T  EBR  =  A  BS (SSL (I)  )  +ABS (SS (I) ) ♦ A  EH 

595. 

IF  (TERR.  LE.0. 0)  GO  TO  300 

t  596. 

EEBR  =  ABS { (—2090. 0*DDL  (I) ♦ (2 1970.0* A4 (I)  - 15043. 0*A5 (I) )  )  ♦ 

f  597. 

-  (22528. 0*A3 (I)  -2736  0. 0*DD (I) )  ) 

598. 

IF  (EEOET.GE. EERR/TERR)  GO  TO  300 

599. 

EEOET=EE2B/TERR 

i  600. 

iir=i 

601. 

300 

CONTINUE 

602. 

IF  (DTFUL. LT. DTHIN) GO  TO  85 

1  603. 

ESTOL=DTFUL*EEOET*SCALL/7524CO.  0 

1  oQ4  . 

IF  (ESTOL.  LE. 1.0)  GO  TO  310 

605. 

C* 
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6  06. 

C* 

ACCURACY  IS  NOT  MET;  MUST  REEUCE  THE  STEP  SIZE 

607. 

C* 

'■>03. 

IE  (DTFUL.LE.DTMIN)  GO  TC  330 

609. 

FR ACI=. 1 

.->10. 

IF  (ESTOL. LT. 59049. 0)  FR ACI=T MUL/ESTOL**. 2 

oil. 

DTACC=FRACI*DTFUL 

612. 

IF  (DTACC.LT. DTMIN)  DTACC=  DT  MI N 

1  3  . 

DTFUL=DTACC 

*•>14. 

GC  TO  240 

615. 

C* 

.•>16. 

C* 

ACCURACY  IS  ACCPTABLE ;  INCREASE  THE  STEP  SIZE 

|6  17. 

C* 

6  13. 

3  10 

FRACI=5.0 

619. 

IF  (ESTOL. GT.1.889568E-4)  F R ACI=T MU L/ESTCL** . 2 

16  20. 

DT ACC=FRACI*DTFUL 

1621. 

IF  (DTACC. GT. DTMAX)  DTACC=CTM AX 

622. 

IF  (DTACC.LT. DTMIN)  DTACC=DTMIN 

|  623. 

GC  TO  85 

I  b  24  . 

330 

IF(IERR)  90,350,340 

625. 

340 

WRITE (NPRTR , 970)  ILR, TNOW 

626. 

CALL  ERROR  (18,1) 

I  627. 

350 

WHITE (NPRTR, 970)  ILR, TNOW 

'  628. 

GO  TO  85 

629. 

C* 

1  630. 

C* 

RON  IS  COMPLETED 

1  631. 

C* 

632. 

1000 

IF  (NDTPST.EQ.O)  GO  TO  1003 

.  633  . 

1001 

DO  1002  1=1, NDTPST 

1  634. 

LDBL=LFDSB*(I-1) *9 

■  635. 

CALL  DCLCT(LDBL) 

636. 

1002 

CONTINUE 

■  637. 

1003 

IF  (NCTPST.EQ.O)  RETURN 

1  638. 

CALL  SSAVE 

639. 

RETURN 

.  640. 

960 

FORMAT  (/// 36X,6HTNEXT=,E17.9,5X,5HTNOW=,E17.9) 

1  641. 

970 

FORMAT  (/2X,38HSPECIFIED  LOCAL  ERROR  EXCEEDED  FOR  SS(,I3, 

1  642. 

-9H)  AT  TIME,  El  2.  4) 

643. 

980 

FORMAT  (/2X,38HSPECIFIED  TOLERANCE  EXCEEDED  FOR  SS(,I3, 

1  644. 

-9H)  AT  TIME  ,  E 1 2.  4) 

1  645. 

990 

FORMAT  (/2X,38HSPECIFIED  TOLERANCE  EXCEEDED  FOR  DD(,I3, 

646. 

-9H)  AT  TI  ME,  El  2. 4) 

■  6  47  . 

END 

I  648. 

SUBROUTINE  DISC 

’  649. 

DIMENSION  RSET  ( 1 ) 

650. 

CCMMON/LOC/LFI , LLR, LFFB , LFTB , LFES B , LFCS B , LFCTP , LFAE , LC AL 

1  651. 

CCMMON/GE  N1/NRUNS,NFILE,TEEG,TFIN,ICLF,ICLV, ICLS, TC LEAR , I R Al 

!  652. 

-NT MX, LTMX (99) ,NTAL, NDTPST, NCTPST, NCOUN,NATP, NENT , I  ECO, ISUM 

653. 

COMMON/GSC1/NCRDR, NPRTR, SS (99)  ,DD  (99) , ATRIB (99) ,SSL (99) , JJ, 

:  6  54. 

-DDL (99) ,TNOH ,  XX(99) , DT NOW , ISTOP 

|  655. 

COMMON  IS  ET  ( 1) 

656. 

EQUIVALENCE  (IS ET  (1)  ,RSET(1)) 

6  57. 

C* 

j  653. 

c* 

TAKE  THE  FIRST  ENTRY  OF  THE  CALENDER 

’  659. 

c* 

660. 

J  =  LC  A  L*  1 
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o  6  1  . 

LE  MT=IS  E? ( J) 

o62. 

L  =  LC  AL 

b63. 

CALL  ULIN  K (L, LENT) 

664. 

C* 

t>6  5  . 

C* 

TAKE  ITS  EVENT  CODE  AND  PROCESS  IT 

666  . 

C* 

667. 

J=LENT+2 

668. 

ICODE=ISET (J) 

669. 

IF  (ICODE.  NE.-1)  GO  TO  11 

o70. 

CALL  CLEARS 

671. 

RETURN 

672. 

C* 

673. 

C* 

COPY  THE  ATRI3UTES  INTO  ATRIE  ARRAY 

o  74. 

c* 

675. 

1 1 

IF  (NATR.EQ.O)  GO  TO  16 

676. 

J=LENT+  3 

677. 

LATR=ISET (J) 

678. 

DC  15  K=1 , NATR 

679. 

J=LATR ♦K 

680. 

ATRIB(K) =RSET (J) 

681  . 

15 

CONTINUE 

682. 

16 

CALL  EVENT (ICODE) 

683. 

C* 

684. 

C* 

FOR  ALL  THE  IX  DTPS!  SEE  IF  THEIR  VALUE  IS  CHANGED 

685. 

c* 

IF  SO  COLLECT  STATISTICS 

636. 

c* 

637. 

IF  (NTHX.EQ.O)  RETURN 

688. 

DO  10  K  =  1  ,  NT  MX 

6  89. 

L=LTHX (K) 

690. 

J=L*7 

691. 

I=-  IS  ET  (J) 

692. 

J=L+8 

693. 

N0=ISET(J) 

694. 

LSTAT  =  ISET  { L) 

695. 

J=LST AT+2 

6  96. 

IP(XX  (I)  .  EQ.  RSET(J)  )  RETURN 

697. 

X=XX(I1 

693. 

CALL  DTPST  (X,LSTAT,  NU) 

699. 

10 

CONTINUE 

700. 

RETURN 

701. 

END 

702. 

SUBROUTINE  SSCND 

703. 

CCSNON/GEN2/NEQD,NEQS,  NEQT, DTS A  V, DTMIN, ETMAX, AERR, RERR,  I  ERR,  NSEV, 

704. 

-ISEES,ICS  V  (25)  ,  ICGV  (2  5)  ,ICDIB(25)  rICDV(25)  ,VAL(25)  ,TTOL(25)  , 

705. 

-ISCD (25) 

7  06. 

DO  10  J=1  ,  NSEV 

707. 

ISCD  (J)  =KKRSS  (ICGV  (J)  ,ICDV  (J)  ,  VAL  (J)  ,  ICDIR  (J)  ,TTOL(J)  ) 

708. 

10 

CONTINUE 

709. 

RETURN 

710. 

END 

711. 

FUNCTION  KKRSS  (IKRSG,IKBSD,CArD,LDIR,TCL) 

712. 

CCNHON/GSCI/NCRDR.NPBTR, SS  (99) ,DD(99) , ATRIB(99)  ,SSL(99)  ,JJ, 

713. 

•DDL (99)  ,TNOW,  XX  (99) ,DTNOW,ISTCP 

714. 

CCN1ON/GEN2/NEQD,NE0S, NEQT, DTS A  V, DT NIN , DTNAX , AERR, R ERR, I  ERR , NS EV , 

715. 

-ISEES,ICS  V  (25)  ,  ICGV  (25)  ,ICDIR(25)  ,ICDV(25)  ,VAL(25)  ,TTOL(25)  , 
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.716. 

*717. 

718. 

719. 
1720. 

721  . 
722. 

1723. 
724. 
725. 

1  726. 
727. 
728. 

1  729. 
730. 
731  . 
732. 
733  . 
734. 
735. 
736. 
737. 
738. 
739. 
740. 
741. 
742. 
743. 
744. 
745. 
746. 


1  747. 
!  748. 

749. 

750. 
i  751. 

752. 

753. 
I  754. 
!  755. 

756. 
,  757. 

758. 

759. 

760. 
1  761. 
!  762. 

763. 
I  764. 

765. 

766. 

1  767. 
768. 
769. 
770. 


I 


-ISCD (25) 

DIMENSION  IKRSG  (1)  ,  IKR  S  D  ( 1)  ,  C  ADD  ( 1)  ,LDIR(1)  fTOL(1) 

JKRSG  =  IKRSG  (1) 

JK  RSD=I KRS  D  ( 1 ) 

IF  (JKRSG)  10,260,20 
10  J  KRSG=- JKRSG 

CRSGL=DDL (JKRSG) 

CFSGN=DD (JKRSG) 

GC  TO  30 

20  CRSGL=SSL (JKRSG) 

CR  3GN=SS (JKRSG) 

30  IF  ( JKRSD)  40,50,60 
40  JKRSD=- JKRSD 

CRSDL=DDL(JKRSD) 

CRSDN=DD( JKRSD) 

GC  TO  70 

50  CRSDL=CADD(1) 

CHSDN=C  AD  D  (  1) 

GC  TO  70 

60  CRSDL=SSL  (JKRSD) 

CRSDN=SS (JKRSD) 

70  IF  (CRSGL— CRSDL)  80,260,90 
80  IF  (CRSGN-CHSDN)  260,100,100 
90  IF  (CRSGN-CRSDN)  130,130,260 
100  IF  (LDIR( 1) )  260,110,110 

110  IF  (CRSGN-CRSDN-TOL  (1) )  120,120,200 

120  KKRSS= 1 
GO  TO  160 

130  IF  (LDIR ( 1) )  140,140,260 

140  IF  (CRSGN-CRSDN+TOL  ( 1) )  210,150,150 

150  KKRSS=-1 

160  IF  (ISEES)  270,170,270 
170  IF  (IKRSG ( 1) )  180,270,190 

180  ISS2S=JKRSG* 1000 
GO  TO  270 
190  ISEES=JKR  SG 
GC  TO  270 
200  KKRSS=2 

GC  TO  220 
210  KKRSS=-2 

220  IF  (ISEES)  270,230,230 
230  IF  ( IKRSG (1 ) )  240,270,250 
240  ISEES=-JKRSG- 100  0 
GO  TO  270 
250  ISEES=— JKRSG 
GO  TO  270 
260  KKRSS=0 
270  RETORN 
END 

SOBROOTINE  SSAVE 
DIMENSION  RSET ( 1) 

COMNON/LOC/LFI, LLR,LFFB,LFTB,LFDSB,LFCSB,LFCTB,LFAE,LCAL 
CCNMON/GE  N 1/NRONS,NFILE,TEEG, TFIN ,ICLF , I CLV, ICLS , TCLEAR, IR AN , 
-NTMX, LT  MX (99) , NT  AL, NDTPST, NCTPST, NCOUN, NATF, KENT , I  ECO , ISOM 
CONMON/GSC 1 /NCRDR, NPRTR,3S  (99)  ,DD  (99) , ATPIB (99)  , SSL  (99)  ,JJ, 
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771. 

772. 

773  . 

774. 

775. 

7  76. 

777. 

779  . 

779. 

790. 

781. 

15 

}7  32  . 

16 

783. 

784. 

10 

17  85. 

§796. 

787. 

.788. 

C* 

1789. 

C* 

790. 

C* 

791. 

c* 

1 7  9  2  . 

c* 

1793. 

794. 

■795. 

■796. 

797. 

798. 

■799. 

1300. 

801. 

■  302. 

§30  3- 

304. 

905. 

■  306. 

■307. 

808. 

■  309. 

|310. 

c* 

811. 

c* 

-812. 

c* 

■  313. 

*314. 

315. 

1316. 

1317. 

818. 

c* 

,819. 

c* 

1  320. 

c* 

821 . 

822. 

1 323. 

f  324. 

325. 

-PEL (99) ,TNOW,  XX(99) , DT KCK , IS1C F 
COMMON  ISET(1) 

EQUIVALENCE  ( ISET  ( 1)  ,RSET(1)) 

DC  10  N=  1  , NCTPST 
LCBL=LFCSB*  (N-1)  *9 
LSTAT=ISET(LCBL) 

ITYPE=IS ET  (LCBL+7) 

IF (ITYPE.LT. 0)  GO  TO  15 
X=SS  ( ITYPE) 

GO  TO  16 
X=DD (-IT YPE) 

NT1=ISET  (LCBL+-8) 

CALL  CTPST (X, LSTAT, NO) 

CONTINUE 

RETURN 

END 

SUBROUTINE  DATAIN 

THIS  SUBROUTINE  READS  THE  INPUT  CARDS  AND  SETS 
UP  THE  'ENVIRONMENT'  FOR  THE  MODEL.  FOR  THE  TIMEBEING 
THE  CARDS  HAVE  FIXED  FORMAT  AND  ORDER. 

DIMENSION  K ARD  (  1  3)  ,  RSET  ( 1)  ,  NAME  (6)  ,  INFO  (4) 

DATA  INFO  /1HN,1HW,1H  ,  1HF/ 

DATA  KARD  /4HPR0J, 4 HDISC, 4  HR ANK , 4 HCCNT , 4HSEV N ,4HSTAT , 4HTALY, 
-4HDTPS, 4HCTPS, 4HCOUN, 4HSTRM, 4HSIMU, 4 HINT!/ 

COHMON/LOC/LFI, LLR, LFFB, LFTB, LFDSB , LFCS E, LFCTB, LFAE , LCAL 
CCMMON/GEN 1/NRUNS , N FILE, TBEG, TFIN , ICLF, ICLV, ICLS,TCLEAR,IBAN, 
-NTMX,LTNX (99) , NT  AL, NDTPST, NCTPST, NCOUN , NATR, NENT,IECC,ISUM 
CCMMON/GEN2/NEQD,NEQS,NEQT, CTSAV, DT»IN,DTMAX, AERR , R ERfi , I  ERR , NSE V, 
-ISEES,ICSV  (25)  ,  ICGV  (2  5)  ,ICDIK(25)  ,ICDV(25)  ,VAI(25)  ,TTOL{25)  , 

-ISCD (25) 

CCMMON/GEN3/ITITLE (15) , I A N AL ( 1 0) , ID  ATE ( 5) 

CCMMON/GEN 4/1 SEED (10) ,LSEED(10) 

CCMNON/GEN8/NIH, VALH,  VALL 

COMMON/GSC1/NCRDH, NPRTR, SS (99)  , DO  (99) , ATRIB(99)  ,SSL  (99) , JJ, 

-DDL (99)  ,TNQW,  XX  (99) , DTBOW , ISTOP 
COMMON  ISET  { 1) 

EQUIVALENCE  (ISET  (1)  ,  RSET  (1)  ) 

BEAD  THE  PROJECT  CARD 

READ  (NCRDR,  1001)  LABEL,  (ITITLE(I)  ,1  =  1,15)  ,  (IANAL(I)  ,1=1,10)  , 

-  (IDATE(I)  ,1=1,5) 

WRITE  (NPRTR,  2001)  LABEL,  (ITITLE  (I)  ,  1=1 ,  1  5)  ,  (IANAL  (I)  ,  1=1  , 10)  , 

-  (IDATE  (I)  ,1=1,5) 

IF  (LABEL. NE.KARD(I) )  GO  TO  100 

BEAD  THE  DISCRETE  CARD 

READ (NCRDR,  1002)  LA  BEL, NFIL E , NENT, N ATF 
WRITE (NPRTR, 2002) LABEL, NFILE, NENT, N ATR 
IF  (LABEL.  NE. KARD (2) )  GO  TC  101 
LFEB=LFI 

IF  (NFILE. EQ. 0. AND. N ENT. EQ. 0)  GO  TO  19 
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<3  26. 

IF  (WENT. EQ.O)  GO  TO  117 

132  7. 

IF  (NFILF.  EO.  0)  GO  TO  11 

328. 

DC  10  1  =  1  , NFILE 

.329. 

C* 

■3  30. 

C* 

READ  SANK  CARDS 

*331. 

C* 

6  32. 

READ (NCR DR,  1003)  LAB  EL, IFL , I E A NK , NR AT 

133  3. 

WRITE  (NPRTR, 200  3) LABEL, IFL, IF ANK, NRAT 

•  3  34. 

IF  (LABEL.  NE.  KARD  (3)  )  GO  TO  102 

335. 

C* 

*336. 

C* 

SET  THE  FILE  BLOCK 

■  337. 

C* 

8  33. 

CALL  SETFB (IRANK, NRAT) 

339. 

10 

CONTINUE 

(  340. 

c* 

*841. 

c* 

SET  UP  THE  CALANDEB  AND  THE  ENTRY  POOL 

342. 

c* 

34  3. 

11 

CALL  SETCAL 

344. 

c* 

845. 

c* 

READ  THE  CONTINUOUS  CARD 

846. 

c* 

347. 

19 

READ  (NCRDR,  100  4)  LABEL,  NEQD,NEQS,  NSEV,DT/1IN,  DTflAX,  AERR,  RERR,  IER 

348. 

WRITE (NPRTR, 2004) LA B EL, NEQD, NEQS, NS EV, DTH IN , DTM AX, AERR , R ERR , I ER 

849. 

IF (LABEL. NE.KARD  (4)  )  GO  TO  103 

1350. 

IF(NEQD.LT.O.OR.  NEQS.LT.O)  GO  TO  188 

1  351. 

NEQT=NEQD«-NEQS 

852. 

IF (NEQT.EQ.O.AND.NENT.EQ.O)  GO  TO  189 

.853. 

IF (NEQT.EQ. 0)  GO  TO  29 

354. 

IERR=2 

355. 

IF (IER. EQ. INFO  (1) )  IERR=-1 

356. 

IF (IER. EQ. INFO (2) . OR. IER. EQ. I N FC ( 3) )  IERR=0 

1  357. 

IF  (IER. EQ.INFO(4) )  IERR=1 

!  <58. 

IF (TERR. EQ. 2)  GO  TO  122 

859. 

C* 

.  360. 

C* 

READ  STATE  EVENT  CARDS 

1361. 

c* 

362. 

IF (NSEV.LT. 0)  GO  TO  104 

363. 

IF(DTMIN.LT.O.OR.DTHAX.LT.O.OB.AERR.LT.O.OR.IERR.LT.O)  GOTO  123 

1  364. 

IF (NSEV.EQ. 0)  GO  TO  29 

1  365. 

DC  20  1=1 , NSEV 

366. 

READ (NCRDR,  1005)  LABEL, ICSV (I) , ICGV (I)  , ICDIR (I)  , ICDV  (I) , 

,367. 

VAL  (I)  ,  TTOL  (I) 

368. 

WRITE  (NPRTR,  200  5)  LABEL,  ICSV  (T)  ,  ICGV  (I)  , ICDIR  (I)  ,ICDV  (I)  , 

869. 

VAL  (I)  , TTOL (I) 

370. 

IF(LABEL. NE.KARD  (5) )  GO  TO  104 

(371. 

IF (ICSV (I) . EQ.O)  GO  TO  190 

’  372. 

R=ICGV (I) 

873. 

IF (ABS (R) .GT.25.OR. ICGV (I) . EQ.O)  GO  TO  190 

1374. 

IF (ICDIR (I) .EQ. 1. OR. ICDIR (I) .EQ.O. OR. ICDI P (I) . EQ.-1)  GO  TO  22 

375. 

GO  TO  191 

876. 

22 

IF(ICDV (I) .EQ.O)  GO  TO  20 

.877. 

R=ICDV  (I) 

1  378. 

IF  (ABS  (R)  .GT.25)  GO  TO  192 

*379. 

20 

CONTINUE 

880. 

C* 

I 
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C* 

•2. 

c* 

3. 

29 

;•*. 

5. 

*6  . 

C* 

7. 

C* 

I 

C* 

1  1. 

9. 

i'1. 

I  2. 

"  3. 

■4. 

I5* 

I  o  . 

•7  . 

30 

9. 

C* 

9. 

C* 

0. 

c* 

>1. 

39 

2. 

3. 

)  4  . 

1  5* 

17. 

'8. 

40 

1  9* 

C* 

Vi. 

C* 

1 1 . 

C* 

1 2- 

49 

I  3. 

14. 

l"5* 

6. 

’.7. 

18. 

1  9* 

50 

1  0. 

C* 

21. 

C* 

,2. 

C* 

3* 

59 

:  4  . 

25. 

1  6’ 

1  7. 

28. 

r9* 

0. 

60 

ii. 

69 

)2. 

C* 

1  3’ 

C* 

64. 

C* 

25. 

C* 

1  2  3  4  5  6  7 

123456789) 123 4567890 1234567390 12 34567890 123456 7390 12 3456 789 1123456739012345^ 

C*  HEAD  ST AT  CARD 

C* 

29  READ (NCRDR, 1006)  LAB EL , NT AL , NDTPST , NCTFST , NCOUN, DTSAV 

WRITE (NPRTR, 2006) LABEL, NT AL, NDTPST, NCTPST , NCOUN, DTSA V 
IF  {LABEL. NE. KARD  (6)  )  GO  TO  105 

C* 

C*  READ  TALY  CARDS 

C* 

IF (NTAL.LT.O. OR. NDTPST. LT.O.CR. NCTPST. IT . 0. OR . NCCUN . LT . 0) GO  TO  124 
IF  (NTAL.EQ.O)  GO  TO  39 
LFT3=  LF I 
DC  30  1=1 , NTAL 

READ (NCR DR,  1007) LABEL, N,  (NAME(J) ,J=1,6) ,NU 
WRITE  (NPRTR  ,  2007)  LABEL,  N,  (N  AM  E  ( J)  ,  J=  1 , 6)  ,NU 
IF(LABEL.NE.KARD(7)  )  GO  TO  106 
CALL  SETTAL  (NAME,NU) 

30  CONTINUE 

C* 

C*  READ  DTPS  CARDS 

C* 

39  IF  {NDTPST. EQ.O)  GO  TO  49 

LF  DSB=LFI 
DO  40  1=1, NDTPST 

READ (NCRDR,  1008) LABEL, N,  (NAME (J)  ,J=1,6) , ITYPE,NU 
WRITE  (NPRTR, 2008)  LABEL, N, (NAME ( J) ,J=1,6)  ,ITTPE,NU 
IF(LABEL.NE.KABD(8)  )  GO  TO  107 
CALL  SETDS  (NAME, ITYPE, NO) 

40  CONTINUE 

C* 

C*  READ  CTPS  CARDS 

C* 

49  IF  (NCTPST. EQ.O)  GO  TO  59 

LFCSB=LFI 

DC  50  1=1, NCTPST 

READ(NCRDR,  1008)  L AB EL, N,  (  N A ME ( J) , J=  1 , 6)  , ITYPE, ND 

WRITE (NPRTR, 2008)  LABEL, N,  (NAME  (J) ,J=1,6)  , ITYPE, NU 
IF  (LABEL.  NE.  KARC  (9)  )  GO  TO  108 
CALL  SETCS  ( NAM E, IT YPE , NU) 

50  CONTINUE 
C* 

C*  READ  COUN  CARDS 

C* 

59  IF  (NCOUN.  EQ.O)  GO  TC  69 

LFCT3=LFI 

DC  60  1=1, NCOUN 

READ< NCRDR, 1009) LABEL, N, (NAME (J) , J= 1 ,6)  , LIMIT 
WRITE  (NPRTR, 2009  LAB  EL , N,  (N AH E ( J) , J= 1 , 6)  , LIMIT 

IF(LABEL. NE.KARD (10) )  GO  TO  109 
CALL  SETCT(NAME, LIMIT) 

60  CONTINUE 

69  IF  (NBNT.EQ. 1)  CALL  SETENP 


READ  STRM  CARDS 
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'  6  . 

i  7. 

n. 

i  •). 

!  i . 

32. 

i1* 

1^- 

70 

.5. 

i6. 

7  1 

I7* 

*.3. 

*9. 

1  °* 

C* 

'  1. 

c* 

>2. 

c* 

'3. 

c* 

4. 

80 

-5. 

36  . 

7. 

C* 

3. 

C* 

39. 

c* 

0. 

81 

1. 

>2. 

i  3. 

'  4. 

5. 

1 6  . 

7. 

3. 

>9. 

n. 

1  . 

2. 

’  3  . 

'4. 

5. 

1001 

'6. 

1002 

’7. 

1003 

8. 

1004 

9. 

1005 

30. 

1006 

1  . 

1  C07 

2. 

1008 

33  . 

1009 

.  *  4. 

1010 

i  5. 

1011 

.6. 

1012 

<7. 

2001 

j  8. 

2002 

'  <1. 

2003 

>0. 

2004 

I 


DC  70  1=1  ,  10 

R2AD{NCHDE,  1010)  LABEL, IE, IS 
IF  (LABEL.  NS.  KABE  (1 1)  )  GC  TO  71 
WRITE (NPETR, 2010)  L ABFL , IR , IS 
IF (IR. LT- 1.0R. IR.GT. 10)  CALL  ERROR{13,IR) 
LSEED (IR) =IS 
IF(IS.LT.O)  IS=-IS 
ISEED  (IR)  =IS 
CONTINUE 
GC  TO  80 

WRITE  (NPRTR  ,2011)  LABEL,  I R 
IF  (LABEL.  NE.  KARD(12)  )  GO  TO  111 
NHUNS=IR 
GO  TO  81 


READ  SI  NO  CARD 

READ (NCRDR , 1011) LABEL, NRUNS 
WRITS (NPRTR, 201 1) LABEL, NRU NS 
IF (LABEL.  NE.KARD (12) )  GO  TO  111 

READ  INTL  CARD 

READ (NCRDR, 1012)  LABEl,TBEG , TFIN , I RA N ,IECO, IS  0 H, ICLF, ICLV,ICLS, 
-TC  LEAR 

WRITE  (NPRTR,  20  12)  LABEL,  TBEG,  TFIN,  IRAN,  I  ECO,  ISUfi,  ICLF,  ICLV,ICLS, 
-TCLEAR 

IF  (LABEL.  NE. KARD (13)  )  GO  TO  110 
IF (IRAN-NE- 1 -AND. IRAN. NE.O)  GO  TO  121 
IF  (IECO.NE.  1-AND-IECO.NE.O)  GO  TO  119 
IF  (ISUM.NE.  1.  AND. ISUfi. NE.O)  GO  TO  120 
IF (TFIN. LE. TBEG)  TFIN=V  ALH 
IF  (ICLF. NE. 1. AND. ICLF. NE.O)  GO  TO  113 
IF (ICLV.NE. 1. AND. ICLV. NE.O)  GO  TO  118 
IF  (ICLS.NE. 1. AND. ICLS. NE.O)  GO  TO  114 
IF  (ICLS.EQ. 1) RETURN 

IF  (TCLEAR. GT.TFTN.OR.TCLEAR.IT. TBEG)  GO  TO  115 
RETURN 

FORMAT (A4,5X, 15A2, 10A2, 5A2) 

FORMAT ( A4,5X, 12, 8X, 11,9 X,I 2) 

FORMAT (A4,5X,I2,8X, I1,9X,I2) 

FORMAT (A4 ,5X,I2,3X, I2,3X,I2,8X, 4F10. 5, A  1 ) 
FCRMAT(A4,5X,I2,8X,I3,7X,I2,8X,I3,7X,2F10.5) 

FORMAT (A4,5X, 4 (12, 8X) , F 10. 5) 

FORMAT(A4,5X,I2,8X,6A2,8X,I2) 

FORM AT ( A4 , 5X, 12 , 8X , 6A2, 8X, 13, 7X ,12) 

FORMAT (A4 ,5X,I2,8X,6A2, 8X, 16) 

FORHAT(A4,5X,I2,8X,I11) 

FORMAT (A4,5X, 12) 

FORMAT(A4,5X,2P10.3,5X,6  (I1,4X)  ,5X,F10.3) 

FORMAT (10X, A4,5X, 15A2, 10A2, 5A2) 

FORMAT (10X,A4,5X,I2,8X,I  1, 9X , 1 2) 
FORMAT(10X,A4,5X,I2,8X,I1,9X,I2) 

FORMAT (10 X,A4, 5 X, 12,3 X, T 2,3 X, 12 , 8X,  4F 1 0. 5,  A 1 ) 
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■>1. 

2005 

FCRMAT(10X,A4,5X,I2,8X,I3, 

7X, 12 ,8X, I3.7X.2F1 0. 5) 

,»2. 

2006 

FORMAT(10X,A4,5X,4(I2,8X)  , 

FlO. 5) 

)  3. 

2007 

FORMAT, 10X,A4,5X,I2 ,  8X , 6  A2 

,  8  X ,  1 2) 

)4. 

2008 

FORMAT (10 X, A4,5X,I2,8X,6A2 

,8X,I3,7X,I2) 

1*5. 

2009 

FORMAT, 10 X, A4.5X.I2 ,8X,6A2 

,  8 X  ,  1 6) 

2010 

FORMAT ( 10  X, A4.5X, 12 ,8X, 1 1 1) 

)7. 

2011 

FORMAT (10X,A4,5X, 12) 

g^8. 

2012 

FORMAT, 10 X, A4.5X.2F 10.3 ,5X 

,6  (11 ,4X)  ,5X, F10.3) 

r  * 

100 

WRITE  (NPRTR  ,  200) 

V). 

200 

FORMAT , 10 X, 50H***** ERROR: 

THE  PCRJ  CARD  IS  MISSING 

) 

)1 . 

GO  TO  1000 

1)2. 

101 

WRITE (NPRTR, 201) 

1)3. 

201 

FORM  AT (10X.50H***** ERROR: 

THE  DISC  CARD  IS  MISSING 

) 

)4. 

GO  TO  1000 

1)5. 

102 

WRITE (NPRTR, 202) 

f  )6. 

202 

FORMAT (1  OX, 50H***** ERROR: 

NOT  ENOUGH  RANK  CARDS 

) 

)  7  . 

GO  TO  1000 

)8. 

103 

WRITE (NPRTR, 20 3) 

)9  . 

203 

FORMAT, 1  OX, 5 OH** *** ERROR: 

THE  CONT  CARD  IS  MISSING 

) 

10. 

GO  TO  1000 

1  1. 

104 

WRITE (NPRTR ,204) 

12. 

204 

FORMAT (1 OX, 50H***** ERROR: 

NOT  ENOUGH  SEVN  CARD  OR  NSEV.LT.O 

) 

13. 

GO  TO  1000 

14. 

188 

WRITE  (NPRTR. 288) 

.15. 

2  68 

FORMAT ( 10X,34H*****ERROR: 

NECD.LT.O.OR. NEQS. LT.O) 

16. 

GO  TO  1000 

17. 

189 

WRITE (NPRTR, 289) 

18. 

289 

FORMAT (101,38 ERROR: 

THIS  IS  NEITHER  CONTINUOUS, 

(19. 

-20H  NOR  DISCRETE  MODEL  ) 

ho. 

GO  TO  1000 

21. 

190 

WRITE (NPRTR, 290) 

I  '2. 

290 

FORMAT, 10 X, 50 H***** ERROR: 

IMPROPER  VALUE  FOR  THE  INDEX  OF 

9 

3. 

-22X.20HCROSSING  VARIABLE 

) 

34. 

GO  TO  1000 

25. 

191 

WRITE  (NPRTR, 291) 

p6. 

291 

FORMAT (1 OX, 50H***** ERROR: 

IMFFOFER  DIRECTION 

) 

:7 . 

GO  TO  1000 

29. 

192 

WRITE (NPRTR, 292) 

,29. 

292 

FORMAT ( 10 X, 50H***** ERROR: 

IMPROPER  CROSSED  VARIABLE 

) 

1 10. 

GC  TO  1000 

M. 

105 

WRITE (NPRTR, 205) 

'2. 

205 

FORMAT ( 10X.50H*****  ERROR : 

STAT  CARD  IS  MISSING 

) 

1  3. 

GO  TO  1000 

4. 

106 

WRITE (NPRTH , 206) 

) 

15. 

206 

FORMAT (1 OX, 50H***** ERROR: 

NOT  ENOUGH  TAIY  CARDS 

*  . 

GC  TO  1000 

17. 

1  07 

WRITE  (NPP.TP ,  207) 

19. 

207 

FORMAT (10X.50H***** ERROR: 

NOT  ENOUGH  DTPS  CARDS 

) 

19. 

GO  TO  1000 

0. 

108 

WRITE  (NPRTR, 208) 

.1. 

208 

FORMAT (10X.50H***** ERROR: 

NOT  ENOUGH  CTPS  CARDS 

) 

12. 

GO  TO  1000 

*  3. 

109 

WRITE (NPRTR, 209) 

*  4. 

2  09 

FORMAT ( 10X.50H***** ERROR: 

NOT  ENOUGH  COUN  CARDS 

) 

i5. 

GC  TO  1000 
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4b  . 

1  10 

WRITE (NPRTR, 210) 

47. 

210 

FORMAT (10 X,50H***** ERROR  : 

INTI  CARD  IS  MISSING 

) 

48. 

GO  TO  1000 

49. 

1  1 1 

HRTTE (NPRTR  , 21 1) 

,30  . 

211 

FORMAT (10X,50H*****ERROR: 

SI MU  CARD  IS  MISSING 

) 

hi. 

GO  TO  1000 

i  2. 

113 

WRITE (NPRTR ,213) 

33. 

213 

FORMAT ( 10 X, 50H***** ERROR: 

II I  EG  AL  CODE  FOR  CHARING  THE 

t 

1  * 

-22X, 5HEILES) 

-5. 

GC  TO  1000 

1  14 

WRITE(NPRTH ,214) 

|7. 

214 

FORMAT (10X»50H*****ERROR: 

ILLEGAL  CODE  FOR  CLEARING  THE 

t 

938. 

-22X,20H STATISTICAL  ARRAYS 

) 

39. 

GO  TO  1000 

.0. 

115 

WRITE(NPRTR,215) 

f  1. 

215 

FORMAT (1  OX, 50 H***** ERROR: 

TCIEAR.GT.TFIN  .OR. TCLEAB. LT.TBEG 

) 

>2. 

GC  TO  1000 

-3. 

1 17 

WRITE  (NPRTR, 217) 

;  '  4  . 

217 

FORMAT (10 X,50H*****ERROR: 

DISCRETE  EVENT  BUT  NO  ENTITY 

> 

>5. 

GC  TO  100  0 

>  6  « 

1 18 

WRITE  (NPR TR ,218) 

,7. 

218 

FORMAT ( 1  OX, 50 H***** ERROR: 

ILLEGAL  CODE  FOR  CLEARING  VARIABLES 

) 

■8. 

GO  TO  1000 

>9. 

1  19 

WRITE  (NPR TR ,219) 

>7. 

219 

FORMAT (1  OX, 50H***** ERROR: 

ILLEGAL  CODE  FCR  DELETING  ECHO  REPORT 

) 

j  1. 

GC  TO  1000 

*  2. 

120 

WRITE (NPRTH, 220) 

r3. 

220 

FORMAT (1 OX, 50H*****EFHOR: 

ILLEGAL  CODE  FOR  DELETING  SUMMARY  REPO 

P- 

- 

■ , 2HRT) 

■5. 

GC  TO  1000 

'6. 

121 

WRITE(NPRTR,221) 

P- 

221 

FORMAT (10 X, 4 9H***»* ERROR: 

ILLEGAL  CODE  FCR  INITIALIZING  RANDOM 

-  ,24HSO.  STREAMS  BETWEEN 

RONS) 

>  9. 

GC  TO  1000 

!7. 

122 

WRITE  (NPRTR, 222) 

1  U 

222 

FORMAT ( 10 X, 44 H***** ERROR: 

II I  EG AL  CODE  FCR  ERROR  INDICATOR) 

I  2. 

GC  TO  1000 

3. 

123 

WRITE  (NPRTR  ,  223) 

,  4. 

223 

FORMAT (1 OX, 50H***** ERROR: 

ETMIN.LT. O.OR.DTMAX.LT.O.OF. AEFR.LT.O. 

9 

5. 

-12HOR.RERR.LT.O) 

'  ,6  . 

GC  TO  1000 

7. 

124 

WRITE (NPRTR, 224) 

9  . 

224 

FORMAT (1  OX, 50H, *****ERROR 

:  NTA L. LT.O.OR. NDTPST. LT. 0. OF. NCTPST. LTO 

9. 

-14H.OR.  NCOON.LT. 0) 

). 

1000 

I=SQRT  (-1.) 

1  . 

STOP 

2. 

END 

■3. 

SUBROUTINE  SETCAL 

•4. 

DIMENSION  RSET  ( 1) 

5. 

COMMON/LOC/LFI, LLR, LFFB , LFTE , LF DS B, LFCSB , LFCTB, LFA E, ICAL 

6. 

COMMON/GEN8/NIR, VALH , VALI 

'7. 

COMMON  IS ET  ( 1) 

•9. 

EQUIVALENCE  (ISET ( 1) , RSET  (  1 ) ) 

>. 

I I=LFI*6 

7. 

1=  (II*NIR-  1)  /NIP 
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101. 

102. 

10  3. 
104. 

i  105. 
*1  06. 
107. 
1  OB. 
109. 
1  10. 

11  1. 
1  12. 
113. 


114. 
1115. 
[116. 
1  17. 
1  IB. 

119. 

120. 
1121. 

122. 

123. 

124. 
'125. 

126. 
1 1 27. 

1  128. 
129. 
1  30. 
:  1 3 1 . 


1  35. 

1  36. 
1  37. 
1  38. 
139. 


'  1  40. 

141  . 
'  142. 
1  43. 
1  44. 

145. 

146. 
1  47. 

1148. 
'  149. 

150. 

151. 
1152. 

153. 

154. 

155. 


C* 

C* 

C* 


IF  (I.GT.LLE)  CALL  ERRCF  (1,1) 

LCAL=LFI 

LFI=II 

ASSIGN  THE  RANKING  RULE 

J  =  LC  A  L  ♦  3 
I  SET  (J)  =4 
J  =  LCAL+  4 
I  SET  (0)  =1 
RETURN 
END 

SUBROUTINE  SETFB  (IR A NK, NR A T) 

DIMENSION  RS  ET  ( 1 ) 

CCMMON/LOC/LFI ,LLR,LFFB,LFTB,LFDSB,IFCSB, LFCTE , LF AE , LC  AL 
CCMNON/GEN1/NRUNS, NFILE,TBEG, TFIN,ICLF, ICLV , ICLS , TCLEAR , I R A N , 
- NTMX , LT MX (9  9)  , NTAL, NDTPST , NCI  PST, NCOUN, NATR, N ENT, I  ECO, IS UN 
CCMMON/GE  N8/NIR, VALH,VALL 
COMMON  ISET(1) 

ECOIVALENCE  (ISET  {  1 )  ,  RSET  (  1 )  ) 

II=LFI+6 

1= (IH-NIR-1)  /NIR 
IF  (I.GT.LLR)  CALL  EBRCR  (1,1) 

IF (IRANK.  LT. 1. OR. IRANK. GT. 4)  CALL  ERROR  (3,1) 

I  =  LFI 
J=I+3 

ISET  (J)  =IR  ANK 
IF  (IRANK.  LT.  3)  GO  TO  1 

IF  (NRAT.GT.  NATR.OR.NRAT.LT.  1)  CALL  ERROR  (2,1) 

J=I  +  4 

ISET  (J)  =N  R AT+  1 

LFI=I*6 

RETURN 

END 

SUBROUTINE  SETENP 
DIMENSION  RSET  ( 1 ) 

CCMMON/LOC/LFI,  LLR,  LFFB,LFTB,  LFES  B,  LFCS B  ,  LFCTB  ,  L>’A  E  ,  LC SL 
CC MMON/GE N 1/NR UNS,N FI IE, TB EG, TF IN, ICLF, ICLV, ICLS, TCLEAR, I: AN, 
-NT MX, LT  MX (9  9) , NTAL, NDTPST, NCTPST, NCOUN, NATR, NFNT, I  ECO, IS UM 
C  CMMON/GE  N8/NIR, VALH,VALL 
CO  MMON  IS  ET  ( 1) 

EQUIVALENCE  (ISET  ( 1 ),  RSET  (  1 )  ) 

K  K=0 
I=LFI+4 

J=LLR-  (  NATR*  1) 

II=(I*NIR-1) /NIR 

IF  (II. GT.  J)  CALL  ERROR  (1,1) 

LFAE=LFI 
K  K  =  KK* 1 
L=LFI+1 
ISET  (L)  =1 
L  =  LFI ♦3 
ISET  (L)  =.I 
LFI=I 
L I  R=  J 


10 
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.156. 
|157. 
1  58. 

159. 

160. 
161. 
162. 

163. 

164. 

165. 

166. 
167. 

1  b8  . 

169. 
I  170. 
!  171. 

172. 

173. 

174. 

175. 

176. 

177. 

178. 

179. 

180. 
!  131. 

132. 
183. 
,  184. 
1  195. 
1  196. 
137. 
'  188. 
199. 
190. 


i  . 

193. 
j  194. 

195. 

196. 
.  197. 
I  198. 

199. 

200. 
201. 

s  202. 

203. 

204. 
;  205. 

2  06. 
207. 
f  209. 
•  209. 
210. 


( 


I=LFI*4 

J=LLB- (NATR  +  1) 

I 1= (I  +  NIR-1)/NIR 

IF (II. LE. J)  GO  TO  10 

1  =  1-3 

ISET(I)  =0 

NENT=KK 

RETURN 

END 

SUBROUTINE  SETT  A  L  (  N  A  HE  ,  IU) 

DIMENSION  NAME  (6) 

DIMENSION  RSET ( 1) 

COMMON/LOC/LFI ,LLR,LFFB,LFTB, LF DSB , LFCS B , LFCT E , LFA E ,  LCAL 
CCMMON/GEN8/NIR, VALH,VALL 
COMMON  ISET  ( 1) 

EQUIVALENCE  (ISET  (1)  ,RSET(1)) 

I=LFI*8 

J=LLR-5 

II=(I*NIR-1) /NIR 

IF  (II .  GT.  J)  CALL  ERROR  (1,1) 

ISET  (LFI)  =J 
I I=LFI+ 1 
K  K=LFI  +  6 
DC  10  L=I I ,  KK 
LL=L-II*1 
ISET(L)=NAHE(LL) 

10  CONTINUE 
II=LFI*7 
ISET ( II) = IU 
I I=J+3 

RSET (II) =VALH 
I I=J ♦ 4 

RSET (II) =VALL 

LFI=I 

L  IR= J 

RETURN 

END 

SUBROUTINE  3ETCS  (N AME, ITYPE , NU) 

DIMENSION  NAME (6) 

DIMENSION  RSET  ( 1 ) 

CCMHON/LOC/LFI  ,  LLR,  LFFB  ,  LFTB,  LFDSB,  LFCSE  ,LFCTB  ,  LFA  E  ,  LC  A I. 
CC.MMON/GE  N1/NRUNS,NFILE,TBEG,  TF  IN,  ICLF,  ICLV,  ICLS,TCLEAR,  IR  AN, 
-NTMX,LTMX (99)  , NT AL, NDTPST, NCTPST, NCOUN, NATR, NENT, I  ECO, ISOM 
CCMMON/GEN8/NIR, VALH, VAIL 
COMMON  ISET  (1) 

EQUIVALENCE(ISET  (1)  ,RSET(1)) 

I=LFI*9 

J=LLR-7 

11= (I ♦NIB— 1 ) /NIR 

IF  (II. GT.  J)  CALL  ERROR  (1,1) 

ISET  (LFI)  =J 
II=LFI+ 1 
KK=LFI»6 
DO  10  L=I I , KK 
LL= L-IH-1 


I 
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211. 

212. 

213. 

214. 

215. 

216. 

217. 

218. 

219. 

220. 
221. 
222. 

223. 

224. 

225. 

226. 
2  27. 
228. 

229. 

230. 

231. 
2  32. 
233. 
234  . 
2  35. 

236. 

237. 

238. 

239. 

240. 
241  . 

242. 

243. 

244. 

245. 

246. 

247. 

248. 
249  . 
2  50. 

251. 

252. 

253. 

254. 

255. 
,  256. 

257. 

258. 

259. 

260. 
261. 
262. 

263. 

264. 

265. 

I 


ISET(L) =  N  A  M  E { L  L ) 

10  CONTINUE 
R= ITYPE 

IF  (ITYPE.  EO.  O.OR.ABS  (R)  .GT.99.)  CALI  FRROR  (19,  HYPE) 

I I=LFI+7 
I SET (II) =ITYPE 
I T-LFT+8 
ISEm (II) =NU 
1 1  =  .J  +  3 

R  SET (II ) =  V  ALH 
I I =J*4 

R  SET (II) =VALL 

LFI=I 

L 1 11=  J 

RETURN 

END 

SUBROUTINE  3ETDS  (NAME, ITYFE, NU) 

DIMENSION  NAME  (6) 

DIMENSION  RSET ( 1 ) 

CCMMON/LOC/LFI, LLR, LFFB, LFTB, LF DS B , LFCS P , LFCTB , LFA E , LC A L 
CC  1M0N/GEN1/NRUNS,NFILE,TBEG,TFIN,ICLF,ICLV, ICLS,TCLEAR, IRAN, 
-NT MX, LTMX (99)  , NT  AL, NDTPST, NCTPST, NCOUN, NATR , KENT, T ECO , IS UM 
C  CMMON/GE  N8/NI R, VALH,VAIt 
COMMON  ISET  ( 1) 

ECUIVALF11CE(ISET  (1)  ,RSFT(1)) 

I=LFI*9 
J  =  LLR-7 

II=(I*NIR-1)/NIR 

IF (II. GT. J)  CALL  ERROR  (1,1) 

ISET  (LFI)  =J 
I I  =  LFI*  1 
KK=LFI+6 
DO  10  L=I I , K  K 
LL=L-II+  1 
ISET(L)  =  NA  ME  (LL) 

10  CONTINUE 

IF  (ITYPE.  GT.O.  AND.  ITYPE.  LE.  NFIIE)  GOTO  11 
IF  (ITYPE. EQ.O. OR. ITYPE. LT. -99 . CR . IT YPE.  GT.  NFIIE) 

-CALL  ERROR  (19, ITYPE) 

NTMX=NTMX+  1 
IT  MX (NTMX) =LFI 
GC  TO  12 

11  LF8L=LFFB+ (ITYPE-1) *6 
II=LFBL+5 

ISET (II) =LFI 

12  II  =  LFI«-7 

ISET (II) =IT YPE 
II  =  LFI«-8 
ISET (II) =  NU 
II=J»3 

RSET (II) =VALH 
I I=J  +  4 

R  SET (II ) = V  ALL 
LFI  =  I 
L  LR= J 
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i  266. 

|  267  . 

268. 

.  269. 

I  270. 

*  271. 

272. 

|  273. 

1  274. 

275. 

100 

*  276. 

I  277. 

101 

*  i  278  . 

201 

'279. 

I  280. 

’  <281. 

102 

232. 

202 

2  93  . 

284. 

103 

1235. 

203 

1236. 

287. 

104 

1238. 

204 

1239. 

j  290. 

I1  291. 

105 

1292. 

205 

J  293. 

I  294. 

'  295. 

106 

1296. 

206 

1  297. 

1  298. 

1299. 

107 

.  300. 

2  07 

1301. 

'.302. 

!  303. 

108 

|  304. 

2  08 

1  305. 

1306. 

109 

,1307. 

209 

It  308. 

*1309. 

1  10 

1310. 

210 

1 1  3  1 1  . 

11312. 

1313. 

1  11 

/I  3  14. 

211 

1315. 

1316. 

1317. 

1  12 

1 1 3 18  . 

212 

*1319. 

1  320. 

RETURN 

END 

subroutine  erp.of  (Iccde,i) 

CC  MMON/GS  Cl /NCR  D  R, NPRTR, 3S  (99)  ,DD(99)  ,ATRIB(99)  ,SSL(99)  ,JJ, 
-DDL  (99)  ,TNOW,  XX  (99)  ,  DTNOW  ,ISTOP 
I  £=0 

IF  (ICODE. LT.  1000)  GOTO  100 
ICODE=1 OOO-ICODE 
I  £=  1 

GC  TO  (101, 102, 103 ,104,  105, 106, 107,108,109, 110,111, 112, 113, 
-114,115,116,117,118),  ICCDE 
WRITE  (NPRTR,201) 

FORMAT  (  10X,  37H***** ERROR:  TRFRE  IS  NOT  ENOtJGH  SPACE) 

IF  (IS-  EQ.  1)  GO  TO  998 

GC  TO  999 

WRITE (NPRTR, 202) 

FORMAT ( 10X, 44H*****ERR0R:  RANKING  ATRIBDTE  IS  OUT  OF  RANGE) 

GC  TO  999 

WRITE (NPRTR, 203) 

FORMAT (10X,33H*****ERROR:  RANKING  CODE  IS  WRONG) 

GO  TO  999 

WRITE (NPRTR, 204)  I 

FORNAT(10X,44H*****ERROR:  THE  ENITY  THAT  IS  TC  BE  REMOVED, 
-, 27H  COPIED,  OR  LOCATED  IN  FILE,I4,15H  DOES  NOT  EXIST) 

GC  TO  998 

WRITE (NPRTE ,205)  I 

FORMAT (10X,29H*****ERROR:  THERE  IS  NO  CTPST, 14) 

IF  (IS. EQ.  1)  GO  TO  998 
GO  TO  999 

WRITE (NPRTR, 206)  I 

FORNAT(10X,29H*****ERROR:  THERE  IS  NO  DTPST,I4) 

IF  (I3.EQ.  1)  GO  TO  998 
GC  TO  999 

WRITE (NPRTR, 207)  I 

FORMAT ( 1QX,29H*****  ERROR:  THERE  IS  NO  TALLY, 14) 

IF  (IS. EQ.  1)  GO  TO  998 

GO  TO  999 

WRITE (NPRTR, 208) 

FORMAT (10X,36H*****ERROR:  NO  CONTINUOUS  NO  DICRETE) 

GO  TO  999 
WRITE (NPRTR, 209) 

FORMAT ( 10 X,36H**^**ERRORl  TNEXT<TNOW  ) 

GC  TO  998 

WRITE (N  PRTR ,210)  I 

FORMAT ( 10X, 29H***** ERROR:  THERE  IS  NO  FILE  ,13) 

IF  (IS.  EQ.  1)  GO  TO  998 
GC  TO  999 

WRITE (NPRTR, 21  1)  I 

FCRMAT(10X,33H,*****EBROR:  THEBE  IS  NO  COUNTER  ,14) 

IF  (IS. EQ. 1)  GO  TO  998 
GO  TO  999 

WPITE(NPRTR,212)  I 

FORMAT( 10X,33H,*****ERROR:  THERE  IS  NO  ATRIBUTE , 14) 

IF  (IS.  EQ.  1)  GO  TO  998 
GC  TO  999 
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11321. 

1  13 

W  BITE  (NPRTR  ,  2  1  3)  I 

'  1  322. 

213 

FORMAT (10X,38H*****ERROR:  STRFAM  NO.  IS  OUT  OF  RANGE, 14) 

1323. 

IF  (IS. EQ. 1)  GO  TO  998 

. 1324. 

GC  TO  999 

1325. 

1 14 

WRITE (NPRTR , 21 4) 

1326. 

214 

FORMAT ( 10X,25H*****ERROR:  CRANE  IS  ZERO) 

1327. 

GC  TO  998 

1328. 

1  15 

WRITE (NPRTR, 215) 

1329. 

215 

FORMAT (10X,45H*****ERR0H:  INCORRECT  PARAMETER  SPECIFICATION) 

1330. 

GC  TO  998 

1331. 

1  16 

WRITE  (NPRTR, 216) 

1332. 

216 

FORMAT ( 10X,45H*****  ERROR :  INCORRECT  CUMULATIVE  PROBABILITY  , 

1  333. 

- 1 3HSPECIF ICATION) 

1334. 

GO  TO  998 

1335. 

1  17 

WBITE(NPRTR,217) 

1336. 

2  17 

FORMAT ( 10X,45H*****ERROR:  FATAL  INTEGRATION  TOLERANCF  ERROR) 

1337. 

GO  TO  998 

1338. 

1  18 

WRITE(NPRTR,218) 

1339. 

218 

FORMAT (10X,45H*****ERROR:  FATAL  INTEGRATION  ACCURACY  ERROR) 

1340. 

GC  TO  998 

1341. 

1  19 

WRITE (NPRTR, 219)  ITYPE 

1342. 

219 

FORNAT(10X,11H*****ERROR:,I4,22H  IS  CUT  OF  RANGE  INDFX  ) 

1343. 

GO  TO  999 

1344. 

998 

CALL  SUMHY(O) 

1345. 

999 

I=SQRT (-1.) 

1346. 

STOP 

1347. 

END 

1348. 

SUBROUTINE  SETCT (NA HE, LIMIT) 

1349. 

DIMENSION  NAME  (6) 

1350. 

DIMENSION  RSET ( 1 ) 

1351. 

CCMMON/LOC/LFI ,LLR,LFF3,IFTB,LF  DSB , LFCS  E , LFCT  E, LFAE,LCAL 

1352. 

CCMMON/GEN8/NIR, VALH,VALL 

1353. 

COMMON  IS  ET  (  1) 

13  54. 

EQUIVALENCE  (IS  ET  (1)  ,RSET(1)) 

1355. 

I=LPI+8 

1356. 

11= (I+NIR-1) /NIR 

1357. 

IF  (II.GT.  LLR)  CALL  ERROR  (1,1) 

1358. 

II=LPI+1 

1359. 

KK=LFI+6 

1360. 

DC  10  L=II,KK 

1361. 

LL=L-II* 1 

1362. 

ISET(L) =NAME (LL) 

1363. 

10 

CONTINUE 

1364. 

II=LFI+7 

1365. 

I  SET (II)  =  LIMIT 

1366. 

LFI  =  I 

1367. 

RETURN 

1368. 

END 

1369. 

SUBROUTINE  FILER (IFILE, A) 

i 1370. 

DIMENSION  A (99) 

1371. 

DIMENSION  RSET ( 1 ) 

1372. 

COMMON/LOC/LFI, LLR, LFFB , LFTD, LFDS B, LFCSB ,LFCTB, LFA E , LCAL 

1  1 373. 

COM.MON/GEN1/NRONS,NPILE,TBEG,TFIN,ICLF,ICLV,  ICLS,TCL EAR,  IRAN, 

11374. 

-  NT  NX ,LT  MX (99)  , NT  AL, NDTPST, NCTPST, NCOUN , N ATR , NENT, IECO, ISOM 

1375. 

CC1MON/GEN 5/IE VENT, TIME 
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, 1376. 

1377. 

1378. 

1379. 

.  1380. 

1381. 

1392. 

1383. 

1384. 
1  385. 

13  36. 

1387. 

1388. 

1389. 

1390. 

1391. 

1392. 

1393. 

1394. 

1395. 

1396. 

1397. 

1398. 

1399. 

1400. 

1401. 

1402. 

1403. 

1404. 

1405. 
1  406. 

1407. 

1408. 

1409. 

1410. 

1411. 

1412. 

1413. 

1414. 

1415. 

1416. 

1417. 

1418. 

1419. 

1420. 

1421. 

1422. 

1423. 

1424. 

1425. 

1426. 

1427. 

1428. 

1429. 

14  30. 

I 


C* 

c* 

c* 


c* 

c* 

c* 


c* 

c* 

c* 


10 

c* 

c* 

c* 

c* 

c* 

c* 

12 


c* 

c* 

c* 

1 

c * 
c* 
c* 

2 


c* 

c* 

c* 

3 


COMMON  ISFT(I) 

EQUIVALENCE  (ISET (1)  , RSET  (  1 ) ) 

TAKE  ONE  ENTRY  FROM  THE  FCOL 

IF  (LFAE.EQ.O)  CALL  ERROR  (1001,1) 

LENT=LFAE 
J=LENT  *1 
L£  =  ISET  (J) 

LFAE=LS 

FIND  THE  ATTRI BITES 

J=LENT+3 
LATR=ISET (J) 

LOAD  A  INTO  THE  ATTRIBUTE  ZCNF 

IF  (NAT R.EQ.O)  GO  TO  12 
DO  10  K=1 , NATR 
J=LATR»K 
RSET(J) =A (K) 

CONTINUE 

FIND  THE  FILE  BLOCK;  FIND  THE  RANKING  ATTRIBUTE; 

COPY  IT  IN  THE  LATE  POSITION; 

IF  NRAT  IS  1  THEN  FILE  IS  CALANDEB  AND  TIME  AND  EVENT 
MUST  BE  COPIED 

IF  (TFILE. GT.  (NFILE+ 1) .OR.IFILE. LT.1)  CALL  ERROR  ( 1 0 1 0 , IFTLE) 
LF3L=LFFB* (IFILE-1) *6 
J=LPBL*3 
IF ANK=ISET  (J) 

GO  TO  (1,2, 3,3)  , IR ANK 

FIPO 

L  =  0 

GC  TO  40 

LIFO;  ADD  IT  TO  THE  TOP 

J=LFBL+1 
L  =  ISET(J) 

GC  TO  40 

COPY  THE  VALUE  OF  THE  RANKING  ATTRI EUTE  IN  ISET(LATR) 

J-LFBL+4 
NBAT=ISET (J) 

IF  (NRAT. EQ. 1)  GO  TO  5 
K=LATB* (NRAT-1) 

RS3T (LATR) =RSET (K) 

GC  TO  6 
CONTINUE 


5 
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.1431. 

C* 

1 1  432 . 

C* 

*1433. 

C* 

1434. 

1  1  435. 

■1436. 

1437. 

6 

-  1438. 

C* 

■  1439. 

C* 

"  1440. 

C* 

1441  . 

I  1442. 

1 1443. 

1  444. 

I  1445. 

1  1446. 

7 

1447. 

1448. 

1449. 

C* 

1450. 

C* 

1451. 

c* 

1452. 

t  \ 

•  1453. 

1454. 

. 1455. 

1456. 

'  1457. 

1458. 

| 1459. 

c* 

| 1460. 

c* 

1461. 

c* 

. 1462. 

4 

I  1463. 

' 1464. 

1465. 

11466. 

1 1467. 

1468. 

,1469. 

40 

11470. 

'1471. 

1472. 

11473. 

1 1474. 

1475. 

|  1 476  . 

1  1477. 

C* 

14  78. 

c* 

, 1479. 

C* 

[ 1480. 

c* 

'  1481. 

c* 

1482. 

c* 

1  1483. 

c* 

I  1484. 

1485. 

THIS  IS  EVENT  CALANDER 

RSET(LATR) =T IMS 
J=LENT+2 
ISET(J) =IE VENT 
CONTINUE 

FIND  THE  POSITION 

J^LFBL+1 
LLFE=ISET (J) 

IF(LLFE.NE.O)  GO  TO  7 
L-0 

GO  TO  40 
CONTINUE 
L=LLFE 

IF (IRANK.EQ.4)  GO  TO  4 

THIS  IS  HVF 

I*L  +  3 

J=ISET(I) 

IF  (RSET (J) . LT.RSET (LATR) )  GO  TO  40 
I  =  L*1 
L=  IS  ET  (I) 

IF  (L.EQ.O)  GO  TO  40 
GO  TO  11 

THIS  IS  LVF 


I  =  L+3 
J -ISET  (I) 

IF  (RSET  (J) . GT. RSET (LATR) )  GO  TO  40 
I  =L+  1 
L=ISET (I) 

IF(L.EQ.O)  GO  TO  40 
GO  TO  4 

CALL  LINK (LFBL.L, LENT) 

RETURN 

END 

SUBROUTINE  LINK ( LFB L, L, LENT) 

DIMENSION  RS ET ( 1 ) 

CCNHON/LOC/LFI • LLR, LFFB , LFTB , LFDSB , LFCS B , LFCTB , LFAE, LC AL 
CCMMON  ISET  ( 1) 

EQUIVALENCE  (ISET(1)  ,RSET(1)  ) 

THIS  SUBROUTINE  LINKS  THE  ENTITY  AT  LENT  TO  THE  FILE  AT 
LFBL,  AT  ADDRESS  L. 

IF  L=0  THE  ENTITY  IS  LINKED  TO  THE  BOTTOM  OF  THE  FILE 
ALSO  THEN  IT  HOST  BE  CHECKED  IF  THER  IS  ANY  ENTITY  IN 
ALREADY  IN  THE  FILE. 

IF  (L. NE.O)  GO  TO  1 
J=LPBL* 1 
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1  486. 

IF  (ISET(J)  .GT.O)  GO  TO  2 

1487. 

C* 

1  438. 

c* 

THERE  IS  NO  ENTITY  IN  THE  FILE 

1489. 

c* 

t 1 490 . 

ISET  (J)  =LENT 

*1491. 

J=LFBL*  2 

1492. 

ISET  (J)  =LENT 

■1493. 

ISET (LENT) =0 

11494. 

J=LENT4l 

1495. 

ISET  (J)  =0 

1496. 

GC  TO  4 

1 1 497  . 

2 

CONTINUE 

■1498. 

C* 

1499. 

c* 

LINK  THE  ENTITY  TO  THE  BOTTOM  OF  THE  FILE 

| 1500. 

c* 

1 1501. 

J=LFBL*2 

1502. 

LLLE=ISET (J) 

1503. 

ISET ( J) =L  ENT 

I  1504. 

ISET (LENT) =LLLE 

*  1505. 

J=LENT*1 

1506. 

ISET  (J)  =0 

I  1507. 

J=LLLE+  1 

| 1508. 

ISET  (J)  =LENT 

1509. 

GO  TO  4 

.1510. 

1 

CONTINUE 

1 1511. 

c* 

*1512. 

c* 

ADD  IT  TO  THE  LIST  BUT  SEE  IF  IT  IS  TOP  OF  THE  LIST 

1513. 

c* 

11514. 

J=LFBL*1 

1 1 5  1 5. 

LLF3=ISET (J) 

1516. 

IF  (L.  NE.LLFE)  GO  TO  3 

,1517. 

c* 

11518. 

c* 

THIS  IS  TOP  OF  THE  LIST 

1519. 

c* 

1520. 

ISET  ( J)  =LE NT 

11521. 

ISET (LENT) =0 

*1522. 

J=LENT*1 

1523. 

ISET  (J)  =LLFE 

1 1 524. 

ISET(LLFE) =LENT 

1 1 525 . 

GC  TO  4 

1526. 

3 

CONTINUE 

1527. 

c* 

1528. 

c* 

LINK  IT  TO  THE  FILE  AND  DON’T  WORRY  ABOUT  THE  LFE  6 

» 1529. 

c* 

1530. 

LF=ISET  (L) 

'1531. 

ISET (LENT) =LP 

1 1 532  . 

J=LENT* 1 

1533. 

ISET ( J) =L 

1534. 

J=LP*1 

i 1535. 

ISET  (J)  =LENT 

1536. 

ISET  (L)  =L ENT 

15  37. 

4 

CONTINUE 

1  15  38. 

c* 

1 1539. 

c* 

COLLECT  STATISTICS  ON  THE  FIIE,  IF  NECESSARY 

1540. 

c* 

I 
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1541. 

T£ET  { LF3L)  =ISET  (LFBL)  +1 

1542. 

J=LFBL*5 

1543. 

LET3=IS ET  (J) 

1544. 

IF  (LDTB.EQ.  0)  RETURN 

1545. 

LSTAT= ISET (LDTB) 

1546. 

X=ISET (LFBL) 

1547. 

J=LDTB+8 

1548. 

NO=ISET  (J) 

1549. 

CALL  DTPST (X,LSTAT, NO) 

1550. 

RETURN 

1551. 

END 

1552. 

SUBROUTINE  RMOVE  (  I, IFILE, A) 

1  553. 

C* 

1554. 

C* 

TrilS  SUBROUTINE  REMOVES  THE  ITH  ENTRY  OF  THE  FILE  IFILE 

1555. 

C * 

AND  LOADS  ITS  ATRIBUTES  INTC  ARRAY  A. 

1556. 

C* 

1557. 

DIMENSION  A  C  1 ) 

1558. 

DIMENSION  RS ET ( 1 ) 

1559. 

CO  NMON/LOC/LFI ,LLR,LFFB, LFTB, LFDSB , LFCS B, LFCTB, LF AE , LC AL 

1560. 

CO MNON/GEN6/ LFBL, LENT 

1561. 

COMMON  ISET  ( 1) 

1562. 

EQUIVALENCE  (IS ET ( 1 ) , RSET  (  1) ) 

1563. 

CALL  COPY (I, IFILE, A) 

1564. 

CALL  OLINK (LFBL, LENT) 

1565. 

RETURN 

1566. 

END 

1567. 

SUBROUTINE  COPY ( I, IFILE, A) 

1568. 

C* 

1569. 

C* 

THIS  SUBROUTINE  COPIES  THE  ATRIBUTES  OF  THE  ITH  ENTRY 

1570. 

C* 

OF  FILE  IFILE  INTO  A. 

1571  . 

c* 

1572. 

DIMENSION  A ( 1) 

1573. 

DIMENSION  RS  ET ( 1 ) 

1574. 

CCMMON/GEN1/NRONS,NFILE,TEEG,TFIN,ICLFrICLV, ICLS,TCLEAR, IRAN, 

1575. 

- NTMX, LT  MX (99) , NT AL, NDTPSI, NCTPST, NCOUN, NATR, NFNT, IECO, ISOM 

1576. 

COMMON  IS  ET  ( 1) 

1577. 

EQUIVALENCE  (ISET  ( 1 ),  RS  ET  (  1)  ) 

1579. 

IF  (NATR.  LE.O)  RETURN 

1579. 

CALL  LOCATE  (I, IFILE, LENT) 

1580. 

J  =  LENT*  3 

1581. 

LATR=ISET (J) 

1582. 

DO  6  J= 1 , NATR 

1583. 

K=LATR ♦J 

1584. 

A  (J)  =RSET  (K) 

1585. 

6 

CONTINUE 

1586. 

RETURN 

1587. 

END 

1588. 

SUBROUTINE  LOCATE (I , IFILE , L ENT) 

1589. 

C* 

1590. 

c* 

THIS  SUBROUTINE  FINDS  THE  LOCATION  OF  THE  ITH  ENTRY 

1591  . 

c* 

IN  PILE  IFILE.  THE  LOCATION  IS  RETURNED  IN  LENT. 

1592. 

c* 

1593. 

DIMENSION  A ( 1) 

1594. 

DIMENSION  RSET  ( 1) 

1595. 

CCMMON/LOC/LFI, LLR , LFFB , LFTB , LFDS B , LFCSB , LFCTE , LFAE , LC AL 

I 
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1596. 

1597. 

1598. 

1599. 

1600. 
1601. 
1602. 

1603. 

1604. 

1605. 

1606. 
16  07. 
1608. 

1609. 

1610. 
1611. 
1612. 

1613. 

1614. 

1615. 

1616. 

1617. 

1618. 

1619. 

1620. 
1621. 
1622. 

1623. 

1624. 

1625. 

1626. 

1627. 

1628. 

1629. 

1630. 

1631. 

1632. 

1633. 
.1634. 

1635. 

1636. 

1637. 
1 1 638 . 
•1639. 

1640. 

(1641. 

1642. 

1643, 
.1644, 

1645, 
'  1 646. 
1647. 
11648. 
1 1 649 . 
1650. 


CCMMON/GEN1/NRUNS, NFILE,TBEG, If  IN, ICLF, ICLV, ICL5, TCLEAR , T R AN , 
-NTMX,LTMX  (99)  ,  NT  AL,  NDTPST,  NCTPST,  NCOUN,  NATR,  NENT  ,  I  ECO,  ISIJN 
CCMMON/GEN6/LFBL,L 
COMMON  ISET  ( 1) 

EQUIVALENCE  (ISET  ( 1 ) , RS ET { 1) ) 

IF  (IFILE.GT.  (NFILE*  1)  .OR.IFILE.LT.  1)  CALL  ERROR  ( 10  10,  IFILE) 
LFBL=LFFB+ (IFILE-1) *6 

IF  (I.GT.I5ET  (LFBL)  .OR.I.LT.  1)  CALL  ERROR  ( 4 ,  I  F ILE) 

J=ISET (LFBL) /2 
IF(I.GT.J)  GO  TO  2 

C* 

C*  START  FROM  THE  TOP 

C* 

J=LFBL* 1 
LENT=ISET (J) 

K=  1 

4  CONTINUE 

IF(K.EQ.I)  GO  TO  6 
J=LENT+ 1 
LENT=ISET (J) 

K=K  ♦1 
GO  TO  4 

2  CONTINUE 

C* 

C*  START  FROM  THE  BOTTOM 

C* 

J=LFBL*-2 
LENT=ISET (J) 

K=ISET (LFBL) 

5  CONTINUE 

IF  (  K.  EQ.  I)  GO  TO  6 
LENT=ISET  (LENT) 

K=K-1 
GO  TO  5 

6  L=LENT 
RETURN 
END 

SUBROUTINE  ULINK  (LFBL, LENT) 

C* 

C*  THIS  SUBROUTINE  UNLINKS  ENTRY  AT  LENT  OF  FILE  AT  LFBL 
C*  AND  LINKS  IT  TO  THE  ENTRY  POOL. 

C* 

DT TENSION  RSET(1) 

CCNMON/LOC/LFI , LLR, LFFB ,LFTB, LF DSB, LFCS B , LFCTB,LFA E, LCAL 
COMMON  ISET  { 1) 

EQUIVALENCE  (ISET ( 1 ) , BS ET (  1) ) 

LP=ISET  (LENT) 

J=LENT> 1 
LS=ISET (J) 

C* 

C*  LINK  IT  TO  THE  POOL 

C* 


ISET ( J) =LFAE 
LFAE=LENT 
IS  ET (LENT) =0 
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1651. 

IF(LP.EQ.O)  GO  TO  10 

j  1652. 

IF(LS.EQ.O)  GO  TO  20 

}  1653. 

C* 

1654. 

C* 

UNLINK  FS ON  THE  MIDDLE 

|  1655. 

<:♦ 

1  1656. 

I  SET (LS)  =LP 

1657. 

.1  =  L  P  +■  1 

. 1658. 

TEST  (J)  =LS 

1  1659. 

GC  TO  30 

9  1660. 

10 

CONTINUE 

1661. 

IF(LS.EQ.O)  GO  TO  25 

1  1662. 

c* 

1  1663  . 

c* 

UNLINK  THE  FIRST  ONE 

1664. 

c* 

I  1665. 

I  SET (LS) =0 

I  1666. 

J=LFBL+ 1 

1667. 

I SET  (J)  =LS 

1668. 

GO  TO  30 

I  1669. 

25 

CONTINUE 

!  1670. 

c* 

1671  . 

c* 

UNLINK  THE  ONLY  ITEM 

; 1672. 

c* 

1  1673. 

J=LFBL*1 

1674. 

ISET  ( J)  =0 

1675. 

J=LFB  L*  2 

1 1676. 

ISET (J) =0 

•  1 677 . 

GC  TO  30 

1678. 

20 

CONTINUE 

1  1679. 

c* 

|  1680. 

c* 

UNLINK  THE  LAST  ITEH 

1681  . 

c* 

.  1682. 

J  =  LP  + 1 

11633. 

ISET  (J)  =0 

1 1684. 

J  =  LFBL  +  2 

1685. 

ISET  (J)  =LP 

11636. 

30 

CONTINUE 

| 1687. 

c* 

16e8. 

c* 

COLLECT  STATISTICS  ON  THE  FILE,  IF  NECESSARY 

,1689. 

c* 

Ilb90. 

ISET (LFBL)  =ISET (LFBL) -1 

*1691 . 

J=LFBL+5 

1692. 

LDTB=ISET (J) 

11693. 

IF  (LDTB.EQ. 0)  RETURN 

11694. 

X=ISET (LFBL) 

1695. 

LSTAT=ISET (LDTB) 

.1696. 

J  =  LDTB«-8 

116  97. 

NU=IS ET  (J) 

1698. 

CALL  DTPST(X,LSTAT,NU) 

1699. 

RETURN 

1 1700. 

END 

11701. 

SUBROUTINE  DTPST  (X,  L,  NU) 

1702. 

c* 

11703. 

c* 

UPDATES  THE  INFORMATION  ON  DISCRETE  TIME  PERSISTANT  STATISTICS 

11704. 

c* 

WITH  X.  L  IS  THE  LOCATION  OF  THE  CORRESPONDING  STATISTICS' 

1705. 

c* 

BLOCK.  NU  IS  THE  UNIT  NO  ON  SHICH  THE  STATISTIC 

I 
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1  1706. 
1707. 
1708. 
1709. 

11710. 
1711  . 
1712. 

11713. 
1714  . 
1715. 
1716. 
I  1717. 
1 1718. 
1719. 
1720. 
1721. 
1722. 
1723. 
1724. 
1725. 
1726. 
1727. 
1728. 
1729. 
1730. 
1731. 
1732. 
1733. 
1734. 
1735. 
1736. 
1737. 
1738. 
17  39. 
1740. 
1741. 
17  42. 
1743. 
1744. 
1745. 
1746. 
1747. 
|  17  48. 
1749. 
1750. 
[1751  . 
I  1752. 
1753. 
1754. 
i  755. 
'  1756. 
1757. 
]  758. 
j  759. 
1760. 


C*  MUST  BE  DUMPED. 

C* 

DIMENSION  R  SET  ( 1 ) 

CCMMON/LOC/LFI,LLR,  LFFB,LFTB, LFDSB,  LFCS3,LFCTBr  LFAE,I.CAI. 

CC MMON/GS Cl/NCRDR,NPRTfl,SS  (99)  , DD  (99) , AT RIB ( 99)  , SSL  (99) ,JJ, 

-DDL (99)  ,T  NOW ,  XX  (99 )  , DT NOW , IS TO P 
COMMON  ISET(1) 

ECUIVALENCE(ISET  (1)  ,RSET(1)) 

LE  1  =  L*1 
L  P  2=  L  ♦  2 
L  E  3=L  +  3 
LF4=L+4 
LE5=L+5 
XI=RSET  (LP2) 

RSET  (LP2)  =X 

IF  (X.LT.RSET  (LP3) )  RSET(LP3)=X 
IF  (X. GT.RSET  (LP4)  )  RSET(IP4)=X 
DT=TNOW-RSET (LP5) 

IF  (DT. LE. 0)  RETURN 
XET=XL*DT 

R  SET (L) =RSET  (L)  +  XDT 

RSET (LP 1) =  RSET (LP1) *XL*XDT 

RSET (LP5)  =TNOW 

IF  (NU.LE. 0)  RETURN 

WRITE  (NO)  TNOW , X 

RETURN 

END 

SUBROUTINE  CTPST  (X,  L,  Nn) 

C* 

C*  UPDATES  THE  INFORMATION  ON  CONTINUOUS  TIME  PERSISTANT  STATISTICS 
C*  WITH  X.  L  IS  THE  LOCATICN  OF  THE  CORRESPONDING  STATISTICS' 

C*  BIOCK.  NU  IS  THE  UNIT  NC  ON  WHICH  THE  STATISTIC 

C*  MOST  BE  DUMPED. 

C* 

DIMENSION  RS ET  ( 1 ) 

CC1MON/LOC/LFI ,LLR,LFFB,LFTB,LFDSB, LFCS  B,LFCTE,LFAE,LCAL 
CO MMON/GS Cl/ NCR  DR,  NFRTR,  SS  (99)  ,DD(99)  ,ATFTB(99)  rSSL(99)  ,JJ, 

-DDL (99) ,TNOW,  XX  (99)  ,DTNOW,ISTCP 
COMMON  ISET  ( 1) 

EQUIVALENCE  (ISET  (1)  ,RSET(1)) 

LE  1=L  + 1 
LF2=L+2 
LPJ=L*3 
LP4=L*4 
LP5=t+5 
XL=RSET  (LP2) 

RSET  (LP2)  =  X 

IF  (X.  LT.RSET  (LP3) )  FSET(LP3)=X 
IF  (X.GT.RSET  (LP4)  )  RSET(LP4)=X 
DT=TNOW-BSET  (LP5) 

IP  (DT. LE. 0)  RETURN 
XDT=  (XL*X)  *DT/2. 

X 2 DT= (XL*  XL*X*X) *DT/2. 

R  SET (L) =RS ET  (L)  *XDT 
RSET  (LP  1)  =RSET  (LP1)  +  X2DT 
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1761 . 

F.  £  ET  (  LP  5)  =TNOW 

1762. 

IF  (Nil.  LE.  0 )  RETURN 

1  763. 

WRITS  (NU)  TNOK  ,  X 

1  764  . 

RETURN 

1768. 

ENO 

1  7o6  . 

SfJBFOUTINE  3CHD  ( IEV  ENT  ,  T IM  E  ,  A) 

1767. 

C* 

il  768. 

C* 

SCHEDULES  EVENT  I EVENT  AT  TIMF  TIME  ON  CALANDER 

U  7  6  9  . 

C* 

1770. 

DIMENSION  A  (99) 

1771  . 

CC.1M0N/GEN1/NRUN5,  N  FILE,  TEEG,T  FIN,  ICLF,  ICLV,  ICL3,TCLEAP  ,  IRAN, 

il  772 . 

-  NT  MX , LT  MX (99)  , NT  AL,  NDTPST, NCTPST, NCCUN, NATR, NENT, IECO, ISUM 

1773. 

CCMMON/GEN5/TI,TT 

1774. 

COMMON/GSCl/NCRDR,NPRTR,SS  (99)  , DD  (9 9) , AT  FIB ( 9 9)  , SSL (99)  ,  JJ, 

11775. 

-I  EL  (99)  ,  T  N  0  W  ,  XX  (99  )  ,  DT  NOW  ,  ISTC  F 

Il  776. 

II=IEVENT 

1  777. 

TT=TIME 

*1778. 

N=  NFI LE* 1 

11779. 

CALL  FI  LE  M  (  N , A) 

1730. 

RETURN 

1731  . 

END 

1 !  782- 

SUBROUTINE  TALLY  (N,VAR) 

*1783. 

c* 

1734  . 

c* 

COLLECTS  STATISTICS  ON  TALLY  i  N  USING  VAR 

|1  785. 

c* 

1  786. 

DIMENSION  RS  ET ( 1 ) 

1787. 

COMMON/GEN1/NRUNS,NFILE,TBEG,TFIN,ICLF,ICLV, ICLS, TCL EAR , IRAN 

1788. 

-  NTMX, LT  MX  (99)  ,  NT  AL,  NDTPST, NCTPST, NCOUN, NATR, NENT, IECO, ISUM 

1  789. 

COMMON/LOC/LFI , LLR, LFFB, LFTB, LFDSB, LFCSE ,LFCTB, LFAE, LCAL 

",790. 

COMMON  ISET  ( 1) 

1791  . 

EQUIVALENCE  (ISET ( 1 ) , RSET {  1) ) 

|  792. 

IF  (N.GT.  NTAL.OR.  N.  LT.  1)  CALL  ERROR  (7,N) 

1  793. 

LTBL=LFT3+  (N-1) *8 

1794. 

L=ISET (LTBL) 

1795. 

RFST(L)  =R  S  FT  (L)  +VAR 

I  796. 

J=L«-1 

*,797. 

RSET  ( J)  =R  SET  (  J)  ♦  VAR  *  V  AR 

1798. 

J=  1+2 

1  799. 

R  SET  (J)  =  R  S  ET  (J)  +1. 

1  900. 

J=L*3 

1801. 

IP  (VAR.  LT.  RSET  (J)  )  RSET(J)=VAR 

302. 

J=L*4 

I  303. 

IF  (VAR.  GT.  RSET  (J)  )  RSET(J)=VAR 

i  8  0  4  . 

J=  LTBL+7 

18)5. 

IF  (ISET  (J)  .SO. 0)  RETURN 

1  308. 

NU=IS  ET  (J) 

'  807. 

WRITE  (NU)  VAR' 

1  903. 

RETURN 

'809. 

END 

810. 

SUBROUTINE  COUNT (N, INC) 

1811. 

DIMENSION  RSET  ( 1) 

18  12. 

CCMflON/LOC/LFI,  LLP,  LFFB,  LFTB,  LFDSB,  LFCSB,LFCT  E,  LFAE,  LCAL 

|  313. 

COMMON /GEN 1/NRU NS , N  FILE, TBEG, TF IN , ICLF, ICLV, ICLS, TCLEAR , IRAN 

9  14. 

-NTMX, LTMX (99)  , NT  AL,  NDTPST , NCT PST , NCCUN , N AT R , NE NT , I  ECO , I SU M 

1915. 

COMMON/ GSC1/NCRDR,NPRTR,SS (n9)  ,0^(99)  , AT RIB (99)  , SSL  (99)  ,JJ, 
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1816. 

-D  CL  (99)  ,TNO»,  X  X  (99  )  ,  DT  JiC  W  ,  1ST  C  P 

1817. 

COMMON  ISET(1) 

18  18. 

EQUIVALENCE  (ISET ( 1)  , RSET  ( 1) ) 

1819. 

IF  (N.  GT.  NCOUN.OR.N.  LT.  1)  CALL  ERROR  (1011,1) 

1820. 

LCT3=  LFCT  B-*-  (N-1)  *8 

1321. 

ISET (LCTB)  = I SET ( LCT  B)  *INC 

1822. 

J=LCT3*7 

1823. 

IF (ISET (LCTB) . LT.ISFT  (J)  )  RETURN 

1324. 

ISTOP=  1 

1825. 

RETURN 

1326  . 

ENO 

1827. 

SUBROUTINE  CLEARS 

1828. 

CCMMON/GEN1/NRUNS,N  FILE, T E EG , TF IN , ICLP, ICLV, ICL3,TCLEAR , ISA N , 

18  29  . 

-NTMX,LTMX (99)  , NT  AL, NOT PST, NCTPST, NCCUN, NATR,  N ENT, I  ECO, IS'JM 

1830. 

IF  (NTAL.EQ.O)  GO  TO  100 

1831. 

DO  10  1= 1 , NTAL 

1832. 

CALL  TRONT (I) 

18  3  3. 

10 

CONTINUE 

13  34. 

100 

IF  (NDTPST.EQ.O)  GO  TO  200 

18  35. 

DC  20  1= 1 , NDTPST 

1  836. 

CALL  TRUND  (I) 

1  837. 

20 

CONTINUE 

1  838. 

200 

IF (NCTPST. EQ. 0)  RETURN 

18  39. 

DC  30  1=1, NCTPST 

1840. 

CALL  TRUNC  (I) 

1841. 

30 

CONTINUE 

1842. 

RETURN 

1843. 

END 

844. 

SUBROUTINE  TRUND  (N) 

845. 

C* 

1  846. 

C* 

TRUNCATES  DISCRETE  TIME  PERSISTANT  STATISTICS 

'  347. 

C  + 

848. 

DIMENSION  RSET(I) 

18  49. 

COMMON/ LOC/LFI,LLE, IFFB, IFTB, LFDS3, LFCSB, LFCTR, LFAE, LC AL 

18  50. 

CCM.MON/GEN1/NRUNS,  N  FILE,  T  BEG,  T  FIN,  ICLF,  ICLV,  ICLS,  TCLEA  R  ,  IRA  N  , 

851. 

-NT1X,LTMX (9  9)  , NTAL, NDTPST, NOT  PST , NCCUN , N ATR , N ENT , I  ECO, ISUM 

852. 

CCMMON/GEN8/NIR, VALR, VALl 

1853. 

CONMON/GSC1/NCRDR, NPRTR , SS (99)  ,DD(9  9)  , ATPIB(99)  ,  SSL  (99)  , JJ, 

954. 

-DDL (99) ,TNOS,  XX (99) ,DTNOJrf,ISTOP 

855. 

COMMON  IS  ET  (  1) 

1856. 

EQUIVALENCE  ( ISET { 1 ) , RS ET ( 1 ) ) 

1857. 

IF  (N.GT. NDTPST. OR.  N.LT.  1)  CALL  ERROR  (1006,  N) 

358. 

LCBL=LFDSB+ (N-1) *9 

359. 

LSTAT=ISET (LDBL) 

I860. 

RSET  (LSTAT)  =0 

:  861. 

J  =LSTAT ♦ 1 

862. 

RSET  (J)  =0 

863. 

J=LSTAT+2 

.  864  . 

R  £  ET ( J)  =0 

865. 

J=LSTAT+3 

’  866. 

RSET  (J)  =V ALH 

867. 

J=LSTAT*4 

1  868. 

RSET  (J)  =V A LL 

1  369. 

J=LST  AT  *5 

870. 

RSET  ( J)  =TNCW 

I 
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1371. 

J=L3TAT+6 

.1372. 

RSET  (J)  =TNOW 

1873  . 

RETURN 

18  74. 

END 

11875. 

SUBROUTINE  TRUNC(N) 

*18  76. 

C* 

1  877. 

c* 

TRUNCATES  CONTINUOUS  TIME  PERSISTANT  STATISTICS 

.1878. 

c* 

11879. 

DIMENSION  RSETfl) 

1  830. 

CCMMON/LOC/LFI,LLR, LFPB, LFTE, LF DS 3 , LFCS B, LFCTB , LF AE , LCAL 

1831  . 

CCMMON/GEN 1/ NSU NS , N FILE, TB EG, T F IN , ICLF, ICLV, ICL3 , TCLEAR , IR A N , 

11882. 

-  NT  NX  , LT  MX (9  9)  , NTAL, NDTPST, NCTPST, NC CUN, NATS, NENT, I  ECO, ISUM 

•1833. 

C  C  NMON/GE  N8/NIR ,VALH,VALI 

1884. 

COMMON/GSC1/NCRDR,NPRTR,SS(99)  , DD  (9  9) ,ATPIB(99)  ,SSL(99)  , JJ, 

•  1 8  35  . 

-DDL (99)  ,TNOH,  XX  (99)  ,DTNOH,ISTCP 

1 1  886. 

COMMON  ISET(1) 

1887. 

EQUIVALENCE  (ISET(1) , RSET  (1) ) 

1888. 

IF(N. GT. NCTPST. OR. N. LT.  1)  CALL  ERROR  (1  00  5,  N) 

11889. 

LCBL=LFCS B ♦  (N-1)  *9 

’ 1890. 

LSTAT=ISET (LCBL) 

1891. 

R  SET  ( I.STAT)  =0 

11892. 

J=LSTAT*1 

| 1 893 . 

RSET  (  J)  =0 

1894. 

J=LSTAT+2 

1895. 

RSET (J) =0 

1 1 896  . 

J=LSTAT+3 

*1897. 

R  SET ( J) =V  ALH 

1898. 

J=LSTAT+4 

11899. 

RSET  (J)  =V A LL 

I  900. 

J=LSTAT  +5 

1901. 

R  SET ( J)  =T  NON 

.’9  02. 

J=LSTAT  +6 

I  903. 

RSET ( J) =TNOW 

9  04  . 

RETURN 

1905. 

END 

1  906. 

SUBROUTINE  TRUNT(N) 

1  907. 

c* 

19  08. 

c* 

TRUNCATES  STATISTICS  OF  TALLY  N 

1*9  09. 

c* 

1910. 

DIMENSION  RSET  ( 1) 

1911. 

CCNMON/LOC/LFI , LLR ,LFFB,IFTB, LFDS B , LFCS B , LFCTB, LEA E , LCAL 

1912. 

CCMMON/GE  N  1  /NRUNS,  N'FIL  E,  T  BEG,  T  FIN, ICLF, ICLV,  ICLS,  TCLEAR  ,  IRA  N  , 

1  9  13. 

- NTMX, LT  MX (99)  , NTAL, NDTPST, NCTPST, NCOUN, NATE, NENT, I ECO, ISUM 

1  914. 

CCSMON/GEN8/NIR, VALH, VAIL 

1915. 

COMMON  ISET  ( 1) 

|  9  1 6  . 

EQUIVALENCE  (ISET ( 1 ) , BS ET ( 1) ) 

|  917. 

IF  (N. ST. NTAL. OR. N.LT. 1)  CALL  EEROF (  1007, N) 

1918. 

LTBL=LFTB ♦ (N-1) *8 

1919. 

LSTAT=ISET (LTBL) 

I  920. 

RSET (LSTAT) =0 

*921. 

J=LSTAT*1 

1922. 

RSET  (J)  =0 

1  923. 

J=LSTAT*2 

|  924. 

RSET  (J)  =0 

1925. 

J=LSTAT ♦3 

I 
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1  926  . 

RS2T  (J)  =V  ALH 

1927. 

J=L3TAT*4 

1923. 

R£ET(J)  =V  ALL 

1929. 

RETURN 

,1  930. 

END 

*1931. 

FUNCTION  TAVG(I) 

1932. 

C* 

1  933  . 

C* 

FINDS  THE  AVERAGE  OF  TALLY  I 

11934. 

c* 

1935. 

DIMENSION  R S ET  (  1 ) 

1936. 

CCMMON/LOC/LFI, LLR, LFFB , LFTE , LF DS B, LFCSB, LFCTB, LFAE, LCAL 

11937. 

CCMMON/GEN1/NRUNS, NFILE, T E EG,  T F IN , ICLF , ICLV, ICLS,TCLEAR, IRAN , 

■1938. 

- NTMX, LT  MX (99)  , NT  AL, NDTPST, NOT  PST, NCCUN, NATR, NENT, IECC, ISOM 

1939. 

COMMON  IS  ET  { 1 ) 

|1940. 

EQUIVALENCE  (ISET ( 1 ) , RS ET  (  1 ) ) 

|1941  . 

IF  (I.GT.NTAL.  OR.I.LT.1)  CALL  ER  ROR  ( 1 007  ,  T) 

1942. 

L1BL=LFTB* (1-1) *8 

1943. 

LSTAT=I SET (LTBL) 

j 1 9  44  . 

J=L5TAT+2 

1  1945. 

T A VG=RS  ET (LSTAT) /RSET (J) 

1946. 

RETURN 

J 1 9  47. 

END 

j 1 948  . 

FUNCTION  TSTD(I) 

1949. 

c* 

1950. 

c* 

FINDS  THE  STANDARD  DEVIATION  OF  TALLY  I 

|l951  . 

c* 

1 1 952. 

DIMENSION  RSET  (  1 ) 

1953. 

COMMON/LOC/LFI , LLR, LFFB, LFTB, LFCSB, LFCSB, LFCTB, LFAE, LCAL 

1  954. 

CC MMON/GEN 1 /NRONS, NFILE, TBEG,TFIN, ICLF, ICLV, ICLS,TCL EAR, IRAN, 

1  955. 

-NTMX  ,  LT  MX  (99)  ,  NT  AL,  NDTPST,  NCTPST,  NCOtlN,  NATR  ,  NENT,  IECO,  ISUM 

1956. 

COMMON  ISET  ( 1) 

.’957. 

EQUIVALENCE  (IS ET ( 1 ) , RS ET ( 1 ) ) 

I  953. 

IF  (I.GT.NTAL.OR.I.LT. 1)  CALL  ERROR  ( 1007, I) 

i  959  . 

LTBL=LFTB+ (1-1) *8 

1960. 

LSTAT=ISET (LTBL) 

1  961. 

J=LST AT+2 

1  962. 

B  N=RS  ET (J) 

1963. 

J=LSTAT+1 

j1  9  6  4  . 

TSTD= (RSET (J)-RN*RSET  (LSTAT) )/ (RN-1) 

1  965. 

IF  (TSTD.LT.O)  TSTD=0.0 

*  966. 

TSTD=SQHT (TSTD) 

1967. 

RETURN 

1  968. 

END 

1  969. 

FUNCTION  TNUM  (I) 

970. 

c* 

I  971. 

c* 

FINDS  THE  NUMBER  OF  TALLIES  IN  TALLY  I 

1  972. 

c* 

973. 

DIMENSION  RSET  ( 1) 

.  974. 

COMMON/LOC/LFI, LLR, LFFB, LFTB, LFCSB, LFCSP, LFCTB, LFAE, LCAL 

I  975. 

CCMMON/GEN1/NRUNS , NFILE,TBEG,TFIN,ICLF,ICLV, ICLS.TCLEAR , IRAN, 

1  9  76. 

-NTMX,LTHX (99) , NTAL, NDTPST, NCTPST, NCCUN, NATR, NENT, I  ECO, ISOM 

977. 

COMMON  IS  ET  ( 1) 

1978. 

EQUIVALENCE  (IS ET ( 1 ) , RSET ( 1 ) ) 

1  979. 

IF  (I.GT.NTAL.OR.I.LT.  1)  CALL  EERO  R  {  1 007 ,  T) 

930. 

LTBL=LFTB* (1-1) *8 

I 
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.1981  . 
II  982. 

1983. 

1984. 
j  1 9  8  5  - 
*19  96. 

1937. 

11983. 
1939. 
1990. 

11991. 
1992. 
1993. 
1994. 
11995. 
|1  9  96. 
1997. 
,1999. 
11999. 
2000. 
2001  . 
2002. 

'  2003. 
2004. 
1 2  005. 
rooe. 

2007. 

2008. 

12009. 
2010. 
2011. 
12012. 
f0 13. 
2014. 
7015. 
11016. 
'2017. 
2018. 
12019. 
|  020. 
2021. 
,2022. 

023. 

\024. 

2025. 

1:026. 

-.’027. 


2028. 
,7029. 
j  030. 
20  31. 
2032. 
'  033. 
1  034. 
7035. 

I 


LSTAT=ISET (LTBL) 

J=  LST AT  +  2 
TNUM=RS  ET (J) 

RETURN 

END 

FUNCTION  THIN  (I) 

C* 

C*  FINDS  THE  MINIMUM  OF  TALLY  I 

C* 

DIMENSION  RSET  { 1) 

COMMON/LOC/LFI  ,LLR,  LFFB,  IFTB,  LFCSB,TFCSE,  LFCT  5,LFA’P,  LCAL 
C C MMON/GE N 1/NSUNS, N FILE , TBFG,T FIN, ICLF,ICL7, ICL3.TCL EAR , IRAN , 
- NTMX , LT  MX (99)  ,  NT  AL, NDTPST , NCT PST , NCCUN , N ATP, NENT , I  ECO, I SUM 
COMMON  IS  ET  {  1 ) 

EQUIVALENCE  (ISET ( 1 )  , RSET  (  1) ) 

IF  (I.  GT.  NT  AL.  OR.  I.  LT.  1)  CALL  EB ROR ( 1 007 , I) 

LTBL=LFTB«-  (1-1)  *8 
LSTAT=ISET (LTBL) 

J=LSTAT+3 
TMIN=RSET (J) 

RETURN 

END 

FUNCTION  THAX  (I) 

C* 

C*  FINDS  THE  MAXIMUM  OF  TALLY  I 

C* 

DIMENSION  RSET (1) 

COMNOM/LOC/LFI , LLR , LFFB , LFTB, LFDS  B, L PCS B, LFCT E, LFAE , LC AL 
CO MMON/GE N1/ NR U NS ,NFILF,TEEG,  TFIN,ICLF, ICLV, ICLS, TCLEAR , IR AN , 
-NTMX.LTMX (99)  , NT  AL. NDTPST, NCTPST, NCOUN , N ATR, NENT, IECO, ISUM 
COMMON  IS  ET  ( 1 ) 

EQUIVALENCE  (ISET ( 1 ), RSET  {  1 ) ) 

IF  (I.GT.NTAL.  OP.I.LT.  1)  CALL  ERROR  { 1 007,  I) 

L13L=LFTB4  (1-1) *8 
LSTAT=ISET (LTBL) 

J=LSTAT+4 
TMAX=HS  ET (J) 

RETURN 

END 

SUBROUTINE  DCLCT  (LDBL) 

C* 

C*  THIS  IS  FOP  COLLECTING  THE  DISCRETE  TIME  PERSISTANT 

C*  STATISTICS 

C* 

DIMENSION  RSET(1) 

COMMON/GSC 1/NCRDR, N PRTR , SS (99) ,DD  (99) , ATRIB(99)  ,SSL  (99) , JJ, 
-DDL (99)  ,T NON ,  XX  (99) ,DTNOS, ISTOP 
COMMON  ISET  (1) 

EQUIVALENCE  (ISET ( 1 ), RSET  ( 1 ) ) 

LSTAT=ISET (LDBL) 

I 3YPE=ISET (LDBL+7) 

IF  (ITYPE.  LT.O)  GO  TO  10 
X  =  NQ  (ITYPE) 

GO  TO  20 

10  X=XX (-ITYPE) 
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:  7  0  3  6  . 

20 

12037. 

2038. 

2039. 

12040. 

"2041. 

C* 

2042. 

c* 

12043. 

c* 

|2044. 

c* 

2045., 

.2046. 

■2047. 

*2048. 

2049. 

12050. 

1  2051. 

2052. 

2053. 

2054. 

2055. 

10 

2056. 

1 1 

2057. 

2058. 

2059. 

2060. 

2061. 

c* 

2062. 

c* 

2063. 

c* 

|2064. 

2065. 

2066. 

12067. 

>  2068. 

2069. 

,2070. 

.2071  . 

2072. 

2073. 

12074. 

j  2075. 

2076. 

,2077. 

2078. 

2079. 

2080. 

c* 

'  2081  . 

c* 

2082. 

c* 

2083. 

1084. 

i  2085. 

2086. 

2087. 

’  2088. 

2089. 

2090. 

N  t  =  I S  ET  (LDBH-8) 

CALL  DTPST  (X,LSTAT,  NU) 

return 

END 

SUBROUTINE  CCLCT  (LCBL) 

THI3  13  FOR  COLLECTING  THE  CONTINUOUS  TINE  PERSISTANT 
STATISTICS 


DIMENSION  RS  ET ( 1 ) 

CC  MMON/GSC 1/NCRDR,  N  PRTR,  SS  (9  9)  , DD (99)  , ATR IB  (99)  ,3SL  (99)  ,  JJ, 
-DLL (99)  ,TNOW  ,  XX  (99)  ,DTNC«,ISTCP 
COMMON  IS  ET  ( 1) 

EQUIVALENCE  (I  SET ( 1 )  , RS ET ( 1) ) 

LSTA  f=I SET (LCBL) 

I1YPE=ISET (LDBL  +  7) 

IF  (ITYPE.  LT.O)  GO  TO  10 
X=SS  (ITYPE) 

GO  TO  11 
X  =  DD (-ITYPE) 

NU=ISET (LCBL+8) 

CALL  CTPST  (X,LSTAT, NO) 

RETURN 

END 

FUNCTION  D AVG  (I) 

FINDS  THE  AVERAGE  OF  DISCRETE  TIME  PERSISTANT  STATISTICS 
DIMENSION  RSET ( 1 ) 

CC  MMON/LOC/LFI , LLR, LFFB,LFTB, LFDSB, LFCS  B, LFCTB, LFAErLCAL 
CCNNON/GEN1/NRUNS,NFILE,TEEG,TFIN,ICLF,ICLV,ICLS,TCLEAR,IEAN, 
-  NT  MX, LT  MX (99)  , NT  AL, NDTPST, NCTPST, NCOUN , NATR , NENT, I  ECO, ISUM 
COMMON  ISET(1) 

EQUIVALENCE  {ISET ( 1 ) , RS FT  {  1) ) 

IF (I. GT.NDTPST.OR.I . LT.  1)  CALL  ERROR (1006, 1) 

LEBL=LFDSB*  (1-1)  *9 
LSTAT=ISET (LDBL) 

J=LSTAT+6 
T=TNOH-RS  ET (J) 

CALL  DCLCT  (LDBL) 

DAVG=RSET ( LSTAT) /T 

RETURN 

END 

FUNCTION  CAVG  (I) 

FINDS  THE  AVERAGE  OF  CONTINUOUS  TIME  PERSISTANT  STATISTICS 
DIMENSION  RSET  ( 1 ) 

CCMMON/LOC/LFI, LLR, 1FFB,LFTB,LFESB, LFCSB ,LFCTB, LFAE, LCAL 
CCMMON/GEN1/NRUNS, NPILE , TEEG, TF IN  r ICLF , ICLV, ICLS,TCL EAR, IRAN, 
-NTMX,LT  HI (99) , NT  AL, NDTPST, NCTPST, NCOUN , NATR , NENT , IECO, ISUM 
COMMON  IS  ET  ( 1) 

EQUIVALENCE  (ISET  ( 1)  ,  RSET  ( 1) ) 

IF  (I.GT. NCTPST. OR.I.LT.  1)  CALL  ERRO  R  ( 1 005,  I) 

LCBL=LFCSB* (1-1) *9 
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2091  . 
»2092. 
12093. 
2094. 

12095. 
2096. 
2097. 
2098. 

12099. 
2100. 
2101. 

12102. 
2103. 
2104. 
2105. 
-  2106. 
2107. 
2108. 
2109. 
2110. 
2111. 
3112. 
2113. 
3114. 
2115. 
2116. 
i  117. 
2113. 
.’119. 
120. 
'.121  . 
3122. 
|  123. 
I  124. 

125. 
.  126. 
1127. 
'  128. 

129. 
1130. 
1 1  31. 

132. 
•  133. 
1134. 
135. 
136. 
1 1  37. 
1138. 
139. 

I ’40. 
'41. 
142. 
143. 
I  44. 
*  45. 


LSTAT=ISET (LCBL) 

J=LSTAT ♦6 
T=TN08-RSET (J) 

CALL  CCLCT  (LCBL) 

CA VG=RSET (LSTAT) /T 

RETURN 

END 

FUNCTION  DSID  (I) 

C* 

C*  FINDS  THE  STANDARD  DEVIATION  OF  DISCRETE  TINE  PERSISTANT  STATISTIC 

C* 

DIMENSION  RSET  ( 1) 

COMMON/LOC/LFI  , LLR,LFFB,LFTB,  LF  ESB,  LFCSB  ,  LFCTB  r  LFA  E,  LC  AL 
CCMMON/GEN1/NPUNS, NFILE, TBPG, TF IN ,ICLF , ICLV, ICLS f TCLF AR, IR AN , 
-NTMX, LT  NX (99)  , NT  AL, NDTPST, NCTPST, NCCON, NATP , NENT, IECO,ISUM 
COMMON  IS ET  (  1) 

EQUIVALENCE  (IS ET ( 1 ) , RS ET (  1) ) 

IF(I. GT. NDTPST. OR. I. LT.  1)  CALI  EPROP (1006, I) 

LEBL=LFDSB+ (1-1) *9 
CALL  DCLCT  (LDBL) 

LSTAT=ISET(LDBL) 

J= LSTAT +6 
T^TNOH— RSET  (J) 

AV=RSET (LSTAT) /T 
J=LSTAT+ 1 

DSTD=RSET ( J) /T-AV*AV 
IF (DSTD.LT.  0)  DSTD=0 
DSTD=SQRT (DSTD) 

RETURN 

END 

FUNCTION  CSTD  (I) 

C* 

C*  FINDS  THE  STANDARD  DEVIATION  OF  CONTINUOUS  TIME  PERSISTANT  STATISTIC 

C* 

DIMENSION  RSET  ( 1) 

CCMMON/LOC/LFI,LLR, LFFB, tFTB, LFCSB, LFCSB, LFCTB, LFAE,LCAL 
CCMMON/GE  N1/NRUNS, NFILE, TBEG, TFIN ,ICLF , I CLV, ICLS,TCLEAB,IRAN, 
-NTMX,  LT  MX  (99)  ,  NT  AL,  NDTPST,  NCTPST,  NCOUN,  NATR,  NENT ,  I  ECO,  ISU.M 
COMMON  ISET  (1) 

EQUIVALENCE  (ISET ( 1 ), RSET  ( 1) ) 

IF  (I.GT. NCTPST. OR. I.LT.  1)  CALL  E  RRO  R  ( 1  005,  I) 

LCBL=LFCSB ♦ (I— 1) *9 
CALL  CCLCT  (LCBL) 

LSTAT=ISET (LCBL) 

J=LSTAT*6 
T=TNOW-RSET (J) 

A V=RSET (LSTAT) /T 
J=LSTAT+1 

CSTD=RSET (J) /T-A V*A  V 
IF  (CSTD.  LT.  0)  CSTD=0 
CSTD=SQ8T (CSTD) 

BETURN 

END 

FONCTION  DPRD  (I) 

C* 
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'  1  46. 
'147. 
I' 148. 
21  49. 
’  1  50  . 

1151. 
152. 
M53. 

1154. 
155. 
1156. 
.157. 


C*  FINDS  THE  PERIOD  OF  DISCRETE  TIME  PERSISTANS  STATISTIC 

C* 

DIMENSION  RSET  (  1 1 

CCMMON/LOC/LFI, LLB,  LFFB, LFTB, LFESB, LFCSB ,LFCTR, LFAE, LCAL 
COMMON/GEN 1/NEUNS,N FI LE,TBEG,TF IN, ICLF,ICLV, ICLS, TCLFAR, IRAN, 
-NTMC.LTMX (99)  ,  NT AL, NDTPST, NCTPST, NCCUN , N ATR, N FNT , I  ECO, ISUM 
CC1.10N  ISET{1) 

EQUIVALENCE  (ISET ( 1 ) , RS ET ( 1 ) ) 

IF  (I.GT.NDTPST.OR.  I.LT.  1)  CAIL  ERROR  ( 1 006, 1) 

LEBL=LFDSB+  (1-1) *9 
LSTAT=ISET (LDBL) 

J=LSTAT*6 

DFRD=TNOW-RSET(J) 

RETURN 

END 

FUNCTION  CPRD  (I) 

C* 

C*  FINDS  THE  PERIOD  OF  CONTINUOUS  TIME  PERSISTANT  STATISTIC 

C* 

DIMENSION  RSET  ( 1) 

CCMMON/LOC/LFI, LLR, LFFB , LFTB, LF DSB, LFCSB , LFCTB , LFAE , LCAL 
CC  MMON/GE N1/NRUNS,NFIL£, TBEG, TFIN,ICLF , ICLV, ICLS, TCL EAR, IRAN, 
-N1MX,LTMX (99) , NT  AL, NDTPST, NCTPST, NCOUN,NATR,NENT, I  ECO, I SUM 
COMMON  ISET  { 1) 

EQUIVALENCE  (ISET ( 1 ), RSET ( 1) ) 

I F ( I. GT. NCTPST. OR. I .IT. 1)  CAIL  EBROB ( 1005,1) 

LCBL=LFCSBMI-1)  *9 
LSTAT=I SET (LCBL) 

J=LSTAT+6 
CPBD=TNOW-RSET (J) 

RETURN 

END 

FONCTION  DMIN  (I) 

C* 

C*  FINDS  THE  MINIMUM  OF  DISCRETE  TIME  PERSISTANT  STATISTIC 

C* 

DIMENSION  RSET  (  1 ) 

COMMON/LOC/LFI , LLB, LFFB, LFTB, LFDS B, LFCSB ,LFCTB, LFA E, LCAL 
CCNNON/GEN1/NRONS,NFILE,TBEG,TFIN,ICLF,ICLV, ICLS, TCLEAR, IRAN, 
-NTNX,LTNX (99) , NT  AL, NDTPST, NCTPST, NCOUN, NATR , NENT, IECO, ISUM 
COMMON  ISET  { 1) 

EQUIVALENCE  ( IS ET ( 1 ) , RS ET  ( 1) ) 

IF  (I.GT.NDTPST.OR. I.LT.  1)  CALL  ERROR  ( 1006,  I) 

LEBL  =  LFDSBMI-1)  *9 
CALL  DCLCT  (LDBL) 

LSTAT=I SET (LDBL) 

J=LSTAT«-3 
DMIN=RSET (J) 

RETURN 

END 

FUNCTION  CHIN  (I) 

C* 

C*  FINDS  THE  MINI  HUN  OF  CONTINUOUS  TIME  PERSISTANT  STATISTIC 

C* 

DIMENSION  RSET  ( 1 ) 


r-  o 
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2201. 

CCMMCN/LOC/LFI ,  LLR,  LFFB, LFTE,  LF  IS  E,  L FCSB  ,  LFCTB,  L?A  E ,  LC  AL 

1-202. 

C  C  MMON/GE  N 1/NRUNS,NFILE,TBEG,TFIN, ICLF , ICLV, ICL5, TCL EAR , IRAN, 

2203. 

-  NT  MX  , LT MX (99)  , NT AL, N DTP ST, NCTPST, NCOUN, NATR, NENT, I ECO, I SUN 

’204. 

COMMON  IS ET  ( 1) 

1.205. 

EQUIVALENCE  (ISET  (  1 ),  RSET  ( 1)  ) 

*2206. 

IF  (I.  GT.  NCT PST. OF. .  I  .  LT.  1)  CALL  EPBO F  ( 1 00 5,  I) 

2  207. 

LCBL=LFCSB+  (1-1)  *9 

12208. 

CALL  CCLCT  (LCBL) 

122  09  . 

LSTAT=ISET (LCBL) 

2210. 

J=LSTAT+3 

-2211. 

CMIN=RSET  (J) 

■  2212. 

RETURN 

*2213. 

END 

2214. 

FUNCTION  DMAX  (I) 

12215. 

C* 

’2216. 

C* 

FINDS  THE  MAXIMUM  OF  DISCRETE  TIKE  PERSISTANT  STATISTIC 

2217. 

C* 

2218. 

DIMENSION  RSET  {  1 ) 

2219. 

CC MMON/LOC/LFI, LLR, LF FB, LFTB,LFESB,L FCSB, LFCTB, LF A E.LCAL 

2220. 

CC  MMON/GE  N1/NRUNS,NFILE,TBEG,TFIN,ICLF,I CLV ,ICLS,TCLEAR, IRAN, 

2221. 

-  NT  MX, LT  MX (99)  , NT AL, NDTPST, NCTPST, NCOUN , NATR , N ENT , IECO, I  SUM 

2222. 

COMMON  ISET(1) 

2223. 

EQUIVALENCE  (ISET ( 1 ), RSET  (  1) ) 

2224. 

IF  (I. GT. NDTPST. OR. I . LT. 1 )  CALL  ERROR ( 1 006, I) 

2225. 

LEBL=LFDSB+ (I- 1) *9 

2226. 

CALL  DCLCT  (LDBL) 

2227. 

LSTAT=ISET (LDBL) 

,2228. 

J=LSTAT»4 

2229. 

DM AX=RSET (J) 

'2230. 

RETURN 

2231. 

END 

12232. 

FUNCTION  CMAX  (I) 

1 2233 . 

c* 

2234. 

c* 

FINDS  THE  MAXIMUM  OF  CONTINUOUS  TIME  PERSISTANT  STATISTIC 

.2235. 

c* 

1  2236. 

DIMENSION  RSET  ( 1 ) 

*2237  . 

COMMON/LOC/LFI, LLR, LFFB, IFTB, LFESB,LFCSE,LFCTB,LFAE, ICAL 

2238. 

CCMMON/GEN 1/NRUNS, N FILE, TBEG,TF IN, ICLF, ICLV, ICLS, TC LEAR, IRAN, 

12239. 

-NTMX,LTMX (99) , NT AL, NDTPST, NCTPST, NCOUN, NATR, NENT, IECO, ISUM 

12240. 

COMMON  IS  ET  (1) 

2241. 

EQUIVALENCE  (ISET { 1 ), RSET  ( 1) ) 

,2242. 

IF(I. GT. NCTPST. OR.I. IT.  1)  CALL  ERROR  (1005,  I) 

12243. 

LCBL=LFCSBMI-1)  *9 

*2244. 

CALL  CCLCT  (LCBL) 

2245. 

LSTAT=ISBT (LCBL) 

12246. 

J=LSTAT*4 

12247. 

CBAX=RSET (J) 

2248. 

RETURN 

■7249. 

END 

■  22  50 . 

FUNCTION  ATP (LENT, I) 

2251. 

c* 

2  252. 

c* 

FIDS  THE  VALUE  OF  ITH  ATRIBUTE  OF  THE  ENTRY  AT  LOCATION  LENT 

|2253. 

c* 

*2254. 

DIMENSION  RSET  (1) 

2255. 

COMMON/LOC/LFI, LLR, LFFB , LFTB, LF DSB, LFCSB , LFCTB , LFAE, LCAL 

I 
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2256. 

CCM30N/GEN1/NFUNS,NFILE,TEEG,  TFIN,  ICLF,  TCI.V,  ICLS,  TCLEAR  ,  1 3 AN  , 

|257. 

-NTMX, LT  MX (99)  , NT  AL, NDTPST, NCTPST, NCOUN, NATR, NENT , I ECC, ISU M 

1258. 

COMMON  ISET(1) 

1259. 

EQUIVALENCE  (I SET ( 1 )  , RS ET  ( 1 ) ) 

*260. 

IF  (I.  GT. NATH. OR.  I.  LT.  1)  CALL  ERROR ( 10 12, I) 

*261. 

J=LENT+  3 

1262  . 

L ATR=IS ET  (J) 

1:263. 

J=LATR*I 

1.264  . 

ATR=  RSET (J) 

1265. 

RETURN 

-1266. 

END 

11267 . 

FUNCTION  LPRED  (LENT) 

*2268. 

C* 

2269. 

c* 

GIVES  THE  PREDECESSOR  OF  THE  ENTRY  LOCATED  AT  LENT 

Il270. 

c* 

?  1271  . 

DIMENSION  RSET  (  1 ) 

2272. 

COMMON/LOC/LFI ,LLR,LFFB, LFTB, LPDSB , LPCS B , LFCTB , LFA E, LCAL 

, 2273. 

COMMON  IS  ET  (  1) 

1  1274. 

EQUIVALENCE  (ISET ( 1 ) , HS ET  (  1 ) ) 

2275. 

LPRED=ISET(LENT) 

2276. 

RETURN 

1  277. 

END 

! '278. 

FUNCTION  LSUCC  (LENT) 

2279. 

c* 

f.280. 

c* 

GIVES  THE  SUCCESSOR  OP  THE  ENTRY  LOCATED  AT  LENT 

1  281. 

c* 

2282. 

DIMENSION  RSET  (1) 

2283. 

CCMH0N/L0C/LFI,LL8,LFFB, LFTB, LPDSB, LFCSB, LFCTB, LFAE, LCAL 

1:284. 

COMMON  ISET  ( 1) 

>.235. 

EQUIVALENCE  (ISET  ( 1 ),  RSET  { 1)  ) 

2236. 

J=LENT+ 1 

1  287. 

LSUCC=ISET (J) 

|  288. 

RETURN 

2239  . 

END 

.2290. 

FUNCTION  NQ  (IFILE) 

291  . 

c* 

'2292. 

c* 

GIVES  NO.  OF  ENTITIES  IN  FILE  I 

2293. 

c* 

1.294. 

DIMENSION  RSET ( 1) 

|:295. 

COMMON/LOC/LFI,LLR ,LFF3, IFTB,LFDSB, LFCSB, LFCTB, LFAE, LCAL 

2296. 

CCMHON/GEN 1/NRUNS,N FILE, TBEG, TFIN, ICLF, ICLV, ICLS, TCLEAR, IRAN, 

,2297. 

- NTMX, LT  MX (99)  , NT AL, NDTPST, NCTPST, NCCUN, NATR, NENT, I ECO, I SUM 

1298. 

COMMON  IS  ET  (  1) 

*2299. 

EQUIVALENCE  (ISET  ( 1 ), RSET  ( 1 ) ) 

2300. 

IF(IPILE.GT.  (NFILEH). OR. IFILE.LT. 1)  CALL  ER  ROR  ( 10  1 0  ,  IFILE) 

. 2301. 

LFBL=LFFB* (IFILE-1)  *6 

1  2302. 

NQ=ISET  (LFBL) 

2303. 

RETURN 

1 7304. 

END 

'305. 

FUNCTION  LFE  (I) 

2306. 

c* 

2307. 

c* 

GIVES  THE  LOCATION  OF  FIRST  ENTRY  OF  FILE  I 

| 2308. 

c* 

1 2309. 

DIMENSION  RSET  ( 1) 

2310. 

CCMMON/LOC/LFI,LLR,LFFB,IFTB,LFDSB, LFCSB, LFCTB, LFAE, LCAL 

I 
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i.311. 

*312. 

313. 

.314. 

§315. 

316. 

317. 

1  318. 

*3  19. 

1320. 

|  321. 

1  322. 

c* 

3  23. 

c* 

’324. 

c* 

1  325. 

’.326. 

!327. 

329. 

’  329. 

330. 

>331. 

’  332. 

.333. 

>334. 

(  335. 

1  336. 

3337. 

j’338. 

I  339. 

:  3  40  . 

.341  . 

1  342. 

1  343. 

>344. 

•  3  45. 

1  346. 

i347. 

>348. 

1349. 

■-350. 

>351. 

■  352. 

|  353. 

!354. 

10 

>355. 

■  356. 

■  357. 

1 1 

>358. 

1  359. 

1  360. 

12 

>361. 

.’362. 

12 

1.363. 

■>364. 

>365. 

15 

I 


CC  MMON/GE  N1/NRUN5,NFILE,1EEG,TFIN,ICLF,ICLV,  ICLS, TCL E AE , IP AN , 
-NTMX,LTMX  (99)  ,  NT  AL,  NDTPST,  NCTPST,  NCOUN  ,  N  ATE  ,  N  ENT ,  I  SCO ,  ISUM 
COMMON  IS  FT  {  1) 

EQUIVALENCE  (ISET  (  1 )  ,  RSET  (  1)  ) 

IF  (IFILE. GT.  (NFILE+I)  .OR. IFILE. LT. 1)  CALL  E? PGR ( 1 0 1 0 , I  FI LE) 
LFBL=LFFB  * (T-1) *6 
J=LF3L* 1 
LFE=I3ET(J) 

R  ETT7B  N 
END 

FUNCTION  LLE  (I) 

GIVES  THE  LOCATION  CF  LAST  ENTRY  OF  FILE  I 
DIMENSION  RSET  ( ' 

COMMON/LOC/LFI, LLR, LFEB, LFTB, LFESB,LFCSB,LFCTB, LFAE, LCAL 
CCMMON/GEN1/NRUNS,NFILE,TBEG,  TFIN,ICLF,ICLV,  ICLS, TCLEAR,  IRAN, 
“NT.MX,  LT  MX  (99)  ,  NT  AL,  NDTPST, NCTPST, NCOUN, NATR , NENT, I  ECO , ISUM 
COMMON  ISET(1) 

EQOIVALENCE  (ISFT ( 1 ) , RSET ( 1) ) 

IF (IFILE. GT.  (NFILE+1)  .OR. IFILE. LT. 1)  CALL  ERROR ( 1 0 1 0 , IFILE) 

LFBL=LFFB* (1-1) *6 

J=LFBL*2 

LLB=ISET(J) 

HETURN 

END 

SUBROUTINE  OUTPUT 

RETURN 

END 

SUBROUTINE  SUMRY  (NNN) 

DIMENSION  RSET  (  1 ) 

CCMMON/LOC/LFI ,LLR,LFFB,LFTB,LFDSB, LFCSB,LFCTB, LFAE, LCAL 
CCNM0N/GEN1/NRUNS,NFILE,TBEG,TFIN,ICLF,ICLV,  ICLS,  TCL EAR, IRAN, 
-N1MX,  LT  MX  (9  9) , NT AL, NDTPST, NCTPST, NCOUN, NATH, NENT, I  ECO, ISOM 
CCMMON/GEN2/NEQD ,NE0S,NEQT, DTS A V, DT MIN, DTMAX , AERR, R ERR, I  ERR, NS EV, 
-IS  EES, ICS  V  (25)  ,  ICGV  (25)  ,ICDIR{25)  ,ICDV(25)  ,VAL(25)  ,TTOL(25)  , 

-ISCD (25) 

CCMMON/GEN3/ITITLE ( 15) ,IANAL(10)  ,IDATE(5) 

CCMMON/GSC1/NCRDR, NPRTR, SS  (99) ,DD(99) ,ATRIB(99)  ,SSL(99) ,JJ, 

-DDL (99)  ,TNOi,  XX  (99)  ,DTNOW,ISTCP 
COMMON  IS ET  ( 1) 

EQUIVALENCE  (ISET  (1)  ,RSET(1)) 

WRITE  (NPRTR,  10) 

FORMAT (1H1,///50X,28HS  OHMART  REPORT) 

WRITE  (NPRTR  ,  11)  (ITITLE(I)  ,1=1,15)  ,  (IANAL(I)  ,1=1,  10)  ,  (IDATE(I)  , 
-1=1,5) , NNN , NRUNS 

FORMAT (//25X,9HPROJECT:  , 1 5 A2 , 1 5X , 9HA NA L YST:  ,10A2/ 

“ , 25 X,6HDAT E:  ,5A2,38X,3HRUN,I3,3H  OF, 13) 

WRITE (NPRTR, 12)  TBEG,TCLEAR,TNOW 

F0HMAT(//25X,28HSIM0LATI0N  STARTED  AT  TIME:  ,E11.2/,25X, 
-28HSTATISTICS  CLEARED  AT  TIME:  ,E11.2, 

-/25X, 14HCURFENT  TIME:  ,E11.2) 

IF  (NCOUN.  EQ.O)  GO  TO  19 
WRITE  (NPRTR,  15) 

FORMAT (////50X, 28H****  COUNTER  INFORMATIN  ****///15X, 5HINDEX, 
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-1366. 

-10X,5HLABEL,20X,5HLIMIT, 10  X  ,  1  3 H COT 3  ENT  VALUE/ 1 5X, 5  { 1 P-)  , 

■13  67. 

-10X,5 (1H-)  ,20X,5  (1H-) ,  10X, 13 (1H-) ) 

1363. 

DC  16  1=1 , NCOUN 

1369. 

LCT3=LFCTB+ (1-11  *8 

11370. 

J  1=LCTB ♦ 1 

■  1371. 

J2=LCTB+6 

2372. 

J3=LCTB  +7 

•1373. 

WRITE (NPRTR, 17)  I, (ISET (J)  ,J=J1 ,J2)  ,ISET (J3)  ,  ISET(LCTB) 

11374. 

17 

FORMAT {/16X,I3, 11X,6A2, 12X,I6,13X,I6) 

2375. 

16 

CONTINUE 

1376. 

19 

IF  (NTAL.EQ.  0)  GO  TO  30 

| 1377. 

WBITE(NPRTR,20) 

1378. 

20 

FCRMAT(/////50X,26H****  TALI!  STATISTICS  ****///15X, 5HINDEX, 8X, 

2379. 

-5HLABEL,9X, 1 1HNO.  OF  OBS. , 1 2X , 4 HM EA N , 1 1 X , 9 HSTD.  DEV.,10X, 

1330. 

-7EMININUM,  10X,7HMAXIHUM) 

1381. 

WRITE (NPRTR, 201) 

1382  . 

201 

FORMAT ( 15  X ,  5 ( 1 H-) , 8 X, 5 (  1 H-)  , 9X , 1 1 ( 1 H-) ,  1 2X , 4  ( 1H-) , 1 1 X , 9  ( 1 H-)  , 

2383. 

-10 X,  7  ( 1  H— )  ,  10X,  7  (1H-)  ) 

1384. 

DC  21  1=1 , NT AL 

2385. 

LT3L=LFTB* (1-1) *8 

2386. 

J1=LTBL*1 

'387. 

J2=LTBL+6 

1388. 

WRITE (NPRTR, 22)  I,(ISET(K) ,K=J1,J2) 

2389. 

22 

FORMAT (/16X,I3,6X,6A2) 

’390. 

LSTAT=ISET (LTBL) 

1391. 

J=LSTAT*2 

2392. 

H=HS  ET  (J) 

2393. 

A V=RSET (LST AT) /RSET (J) 

394. 

STD=  (RSET  (LSTAT+1)  -RSET  (J)  *AV*A  V)  /  (RSET  (J)  -1.) 

'  395. 

IF  (STD. LT. 0)  STD=0 

1396. 

STD=SOHT(STD) 

,'397. 

RMN=R  SET(LSTAT*-3) 

1  398. 

RHX=RSET (LSTAT*4) 

.399. 

WRITE (NPRTR ,23)  N, A  V, STD , RUN, BN X 

1400. 

23 

FORHAT (1H  +  ,43X,I8,8X,E1 1. 2, 6X, E 1 1. 2 , 7X, El  1.2, 6X, El  1.2) 

j  401. 

21 

CONTINUE 

1  402. 

30 

IF  (NDTPST.EQ.O)  GO  TO  40 

1403. 

WRITE (NPRTR, 31) 

i'  404. 

31 

FORMAT (/////43X,45H  ****  DISCRETE  TIME  PERSISTANT  STATISTICS  **** 

405. 

-///15X,5HINDEX,8X,5HIABEL,  10X,4HNEAN, 1 1X,9HSTD.  DEV. , 10X, 

1406. 

-7HMINIM0H ,10X,7HHAXIMUH/15X,5  { 1 H-) , 8X, 5  { 1 H-)  , 10 X , 4 ( 1 H-)  , 1 1X , 

1407. 

-9  |1H-)  ,10X,7(1H-) ,10X,7(1H->) 

|  408. 

DC  32  1  =  1  ,NDTPST 

1.409. 

LEBL=LFDSB> (1-1) *9 

1410. 

J1=LD8L+1 

1  411. 

J2=LDBL*6 

|  412. 

WRITE (NPRTR, 22) I, (ISET(K) ,K=J1,J2) 

1413. 

LSTAT=ISET (LDBL) 

>414. 

T=TNOW-RSET (LSTAT+6) 

1  415. 

AV=RSET (LST&T) /T 

>416. 

STD=RSET ( LSTAT+ 1) /T-AV*  A V 

■417. 

IF(STDvLT.O)  STD=0 

1  418. 

STD=SQfiT(STD) 

\  *19. 

RHN=RSET ( LST AT *3) 

1420. 

RHX=RSET (LSTAT*4) 
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■421. 

WRITE (NPRTR,  34) AV,STD,EMN,FMX 

2422. 

34 

FORMAT  (1H  +  ,38X,  Ell.  2,6X,El1.2,7X,F11.2,fcX,El1.2) 

2423. 

32 

CONTINUE 

2  424. 

40 

IF  (NCTPST.  EQ.O)  GO  TO  50 

14  25. 

WRITE (N PR TR, 41) 

2426. 

41 

FORM AT (/////43  X,47H****  CONTINUOUS  TIME  PERSISTANT  STATISTICS 

2427. 

-///I 5X, 5HINDEX.8X, 5HLABEL, 10X,4HMEAN, 1 1X,9HSTD.  DEV. , 10 X, 

2429. 

-7HMINIMUM,10X,7HMAXTMUN/15X,5{1H-)  ,  8X,5  (1H-)  ,  10X, 4  ( 1  H-)  ,  1 1 X , 

2429. 

-9  ( 1H-) , 1 0  X ,  7 ( 1 H— )  ,  1  0  X ,  7 (1H-)  ) 

2430. 

DC  42  1=1 , NCTPST 

2431  . 

LCBL=LFCS  B  + (1-1)  *9 

2432. 

J  1  =  LCBL  + 1 

2433. 

J2=LCBL+6 

24  34  . 

WRITE (NPRTR, 22) I, (ISET(K)  , K=J1, J2) 

2435. 

LSTAT=ISET (LCBL) 

>436. 

T=TNOW-RSET(LSTAT+6) 

2437. 

A V=RSET  (LSTAT)  /T 

2438. 

STD=RSET(LSTAT+1)/T-AV*AV 

'439. 

IF (STD. LT. 0)  STD=0 

2440. 

STD=SQRT  (  STD) 

2441  . 

RMN=RSET  (LSTAT ♦■3) 

442. 

RHX=RSET(  LSTAT  *-4) 

443. 

WRITE (NPRTR, 44) A  V, STD , RMN , RMX 

2444. 

44 

FORMAT (1H+,38X, Ell. 2, 6X ,  E 1 1 .2, 7Xr E 1 1.2, 6X, El  1.2) 

>445. 

42 

CONTINUE 

446. 

50 

IF(NEQT.EQ.O)  GO  TO  60 

.447. 

WRITE (NPRTH,130) 

448. 

130 

FORMAT (////43X,35HFINAL  VALOE  CF  CONTINUOUS  VARIABLES/, 43X, 

449. 

-35(1  H-)  //) 

450. 

DO  51  1=1 , NEQT 

2451  . 

WRITE  (NPRTR,  131)  I,  SS(I),I,  DD  (I) 

'452. 

131 

FORMAT (35X,3BSS ( , 12 , 2H) =,E10. 2, 18X, 3HDD (,I2, 2H) =,K10.2) 

453. 

51 

CONTINUE 

.454. 

60 

RETURN 

2455. 

END 

4  56. 

FUNCTION  DR AND ( ISTRM) 

457. 

CCMMON/GEN4/IS EED ( 10)  ,LSEED  (10) 

4  58. 

DATA  MULT/65539/ 

459. 

IF (ISTRM.  LT. 1.0R.ISTRM.GT. 10)  CALL  ERROR ( 1 0 1 3 , ISTRM) 

460. 

ISEED (ISTRM) =IS EED ( ISTR M) *  MULT 

461. 

IF (ISEED( ISTRM) )  90,10,100 

462. 

10 

CALL  ERROR  (1014,  1) 

463. 

90 

ISEED (ISTRM) =ISEED ( ISTRM) +2147483647+1 

464. 

100 

DFAND=ISEED  (ISTRM) 

465. 

DR AN D= DR AND*.  46 5661 3 E-9 

466. 

IP  (LSEED  (ISTRM)  .  LT.  0)  DRAND=  1.-DP  AND 

467. 

RETURN 

468. 

END 

469. 

P UNCTION  UNIF(A,B,IST) 

470. 

IF(A.GT.B)  CALL  ERROR  ( 1 0  15  ,  1) 

471. 

ONIP=A+  (B-A)  *DRAND  (1ST) 

472. 

RETURN 

473. 

END 

474. 

FUNCTION  ERLNG  (B,A, 1ST) 

475. 

K=A 

SIMULATION  LANGUAGE 
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.  24 76. 
I  1477. 
'2478. 
2479. 
2  4  80. 
2491  . 
.492. 
.481. 
2  4  94. 

2485. 

2486. 

2497  . 

2498  . 
2499. 
2490. 
249  1  . 
2492. 
>493. 

494. 

2495. 

2496. 

497. 

498. 
2499. 
'500. 

501. 

.502. 

2503. 

504. 

505. 
.506. 
'507. 

508. 

.509. 

2510. 

511. 

512. 
2513. 

514. 

515. 
.516. 
2517. 

518. 

519. 
2520. 

521. 

522. 
'523. 
>524. 

525. 

526. 
2527. 
529. 

‘  529. 
2530. 


1  2  3  4  5  6  7 

1  234567890  1234567890  12  345  67390  12  34567890  123456789'J  1  2  345673901  2  34567390  1  23 

IP  (K.LT.I.OR.B.LT.O)  CALL  ERROR  (101^,1) 

8=1 

DC  10  1=1 , K 

R=R*DFAND(IST) 

10  CONTINUE 

ERLNG=-B* ALOG  (R) 

RETURN 

END 

FUNCTION  EX  PON { A  V, I  ST) 

IF(AV.LT.O)  CALL  ERROR  ( 10  15,  1) 

Y  =  DHAND  (1ST) 

EXPON=—AV* ALOG (Y) 

RETURN 

END 

FUNCTION  FLOGN (QM,QS,IS) 

IF  (QS.LT.  O.OR.QN.LE.O)  CALL  ER ROR  { 1  0  15 ,  1 ) 

QES=QS*QS 

QMS=QM*QM 

XSS=ALOG( OSS/QMS*  1. ) 

XS  =  SQRT  (XSS) 

X !*=A LOG  (QM)  5*XSS 
V A  =8 NORN (XH,XS,IS) 

RLQGN=EXP  (VA) 

RETURN 

END 

FUNCTION  RNORM (XMN,  STD,  ISTRM) 

DIMENSION  ENORM (  10)  ,  IEVEN(IO) 

DATA  IEVEN/ 1 0*1/ 

IE  (STD.  LT . 0)  CALL  ERROR  ( 10  1 5,  1) 

IF  (IEVEN(ISTRM)  .GT.  1)  GC  TO  20 
10  UA  =  2.*DRAND  (ISTRM)  -  1. 

UE  =  2.  *DRAND(  ISTRM)  -  1. 

IEVEN (ISTRM) =2 
W=UA*0A*UB*UB 
IF  (W.GT.1.0)  GO  TO  10 
W=SQRT (-2. *ALOG (H) /H) 

RNORH=UA*« 

ENORM  (ISTRM)  =fJB*W 
GO  TO  30 

20  IEVEN (ISTRM) = 1 

RNORM=ENORN  (ISTRM) 

30  RNORM=RNORN*STD*  XMN 

RETURN 
END 

FUNCTION  TRIAG  (XL, XM,XH, ISTRM) 

IF  (XL.GT.  XM.OR.  XL.  GT.  XH.OR.XM.GT.XH)  CALL  ERROR  {  10  1  5 , 1) 

RN=DRAND( ISTRM) 

BBA=XM-XL 
CM A=  XH-  XL 

IF  (RN-BMA/CM  A)  10,10,20 
10  TPIAG=XL4SQRT(BHA*CHA*RN) 

RETURN 

2  0  TFIAG=XH-SQRT  (CMA*  (  1.-RN)  *  (XH-Xf!)  ) 

RETURN 

END 
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I 


1  2  3  4  5 

123456764 0123456789012345b739012 3 4567 89 01234567390 


6  7 

12345678901234567390123 


*’531. 

■  1532. 

*'533. 

1534. 

11535. 

11536. 

15  37. 

,’538. 

1539. 

1540. 

2541  . 

1542. 

1543. 

2  544  . 

'545. 

10 

'546. 

1547. 

20 

2543. 

'549. 

30 

.550. 

2551  . 

40 

'552. 

50 

553. 

2554. 

60 

>555. 

556. 

70 

-557. 

2558. 

559. 

80 

560. 

90 

2  56  1  . 

'562  . 

100 

563. 

.564. 

1565. 

566. 

'  567. 

1568. 

10 

569. 

1  570. 

.571. 

1572. 

i  573. 

20 

'  574. 

:575. 

30 

,  576. 

1  577. 

-578. 

1579. 

580. 

.581  . 

40 

1582. 

583. 

584. 

'585. 

I 


FUNCTION  KEIIUL(B,A,ISTRM) 

»EIBUL=  (-  B  *  A  LOG  (DP  A  ND  (I  STR  M)  )  )  **(1./A) 

RETURN 

END 

FUNCTION  BETA (A,B,ISTRW) 

X  =  GA!!A(1. ,A.ISTF!1) 

BE TA=X/ (X+GAffA { 1  .  ,  R  ,  I  STEM)  ) 

RETURN 

END 

FUNCTION  GAfl  (AK,  ISTBM) 

K=  AK 
FK  =  K 
GAi=0 

IF (K)  30,30,10 
PROD= 1.0 
DC  20  1=1, K 
PEOD=PROD*DRAND (1ST  EM) 

G AM=-ALOG (PROD) 

DG=AK-FK 

IF  (DG-.015)  100,100,40 

IF  (DG-.985)  60,50,50 

»=1. 

GO  TO  90 
A=1./DG 
B=  A./  (1.-DG) 

X=DRAND (ISTRN)  **A 
Y  =  DR  AND  (ISTRH)  **B*X 
IF  (Y- 1 . )  80,80,70 
H  =  X/Y 

Y=-ALOG (DRAND (ISTRN) ) 

GAM=GAR4W*Y 

RETURN 

END 

FUNCTION  GAMA (B, A, ISTRN) 

IF(A.GE.I.O)  GO  TO  10 

GANA=GAN(A, ISTRN) *B 

RETURN 

K  A  =  A 

AK=KA 

C=  A- A  K 

IF  (A. LT. 5.  )  GO  TO  20 
IF  (DRAND(ISTRH) . LT.  C)  K  A  =  K  A ♦  1 
PR= 1 . 0 

DC  30  1=1, KA 
PR=PR*DRAND(ISTRN) 

G ANA=— A  LOG  (PR) 

IF  (A. GE.5.)  GO  TO  40 

TEST= (GANA/AK) **C*EXP(-C* (GANA/A- 1. ) ) 

IF (DRAND(ISTRH) .GT. TEST)  GO  TO  20 
GANA=GANA* (A/AK) 

GANA=GANA*B 

RETURN 

END 

FUNCTION  NPSSN (P, ISTRN) 

IF  (P.LT.O)  CALL  ERROR  ( 1 0  15 ,  1) 
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II  2  3  a  5  0  7 

1 23'* 567  89*>  12345137090  12  345*7890 12 345678901 2 3*156 71*10  123456  7090 12345673901 2 3 • 


536. 

NF5SN=0 

587. 

IF (  P-9)  10,10,40 

1 :  5€»  8  _ 

10 

Y=EXP  (-P) 

.599. 

X=1.0 

590. 

20 

X  =  X*DRA  ND  ( ISTEN) 

591  . 

IF  (X-Y)  50,30,30 

59  2. 

30 

NP53N=NPSSN+ 1 

15  93. 

3C  TO  20 

594. 

40 

Z  =  SQRT  (P) 

595. 

N  FSSN=  RNO  F  M ( P, Z , 1ST  FM)  +. 5 

'  596  . 

50 

RETURN 

1  597. 

END 

598. 

PO NOTION  DP ROB  (CPF OB, VALUE, NVAF , 1ST VI) 

1599. 

DIMENSION  CPROB(I)  ,  VALUE  (1) 

'  6C0  . 

RN=DEAND(ISTRH) 

601. 

DC  10  1=1 , N VAR 

:602  . 

IF  (RN.  LE.CPROB(I)  )  GO  TO  20 

!603. 

10 

CONTINUE 

604. 

CALL  ERROR  (10 16,  1) 

.605. 

20 

DPROB=V  ALUE (I) 

1606. 

RETURN 

607. 

END 

I 

I 

I 

I 


DISTRIBUTION  LIST 


18  Jun  1981 


"A  Generalized  Simulation  Language  for  Naval  Systems  Modeling" 
All  addressees  receive  one  copy  unless  otherwise  indicated. 


Dr.  B.  Schmeiser 

Department  of  Industrial  Engineering 

Purdue  University 

West  Lafayette,  IN  47907 

Dr.  T.  Schriber 

School  of  Business  Administration 
University  of  Michigan 
Ann  Arbor,  MI  48109 

Dr.  R.  E.  Shannon 

Department  of  Industrial  &  Systems 
Engineeri ng 

University  of  Alabama,  Huntsville 
Box  1247 

Huntsville,  AL  35807 

Mr.  Roger  Willis 
Mr.  Keith  Thorp 
U.S.  Army  TRASANA 

White  Sands  Missile  Range,  NM  88002 

Dr.  R.  Launer 
Army  Research  Office 
Durham,  N.C. 

Dr.  Seymour  Goldstein 
David  Taylor  Naval  Ship  Research 
&  Development  Center  (DTNSRDC) 
Carderock 

Bethesda,  MD  20084 

Mr.  Tony  Rouse 

U.S.  Army  Material  Systems 

Analysis  Activity 

Aberdeen  Proving  Grounds,  MD  21005 

Mr.  Robert  FI  urn 

ASW  Systems  Project  Office 

PM-4 

Department  of  the  Navy 
Washington,  D.C.  20360 


Mr.  Douglas  Kinney 
Code  122 

Naval  Weapons  Center 
China  Lake,  CA  93555 

Defense  Technical  Information  Center 
Cameron  Station 

Alexandria,  VA  22314  2  copies 

Center  for  Naval  Analyses 
2000  North  Beauregard  Street 
Alexandria,  VA  22311 

Office  of  Naval  Research 
Code  434 

800  N.  Quincy  Street 

Arlington,  VA  22217  2  copies 

Naval  Postgraduate  School 
Technical  Library 
Monterey,  CA  93940 

Dr.  Richard  E.  Nance 
Computer  Science,  Head 
Virginia  Polytechnic  Institute 
&  State  University 
Blacksburg,  VA  24061 

Dr.  A.  Alan  B.  Pritsker 
Department  of  Industrial  Engineering 
Purdue  University 
West  Lafayette,  IN  47907 

Dr.  J.  William  Schmidt 
Department  of  Industrial  Engineering 
&  Operations  Research 
VPT  &  SU 

Blacksburg,  VA  24061 

Dr.  Lee  Schruben 
College  of  Engineering 
Cornell  University 
Ithaca,  NY  14853 


Distribution  List  (Cont'd) 

AH  addressees  receive  one  copy  unless  otherwise  indicated. 


Dr.  S.  S.  Lavenberg 

IBM  Thomas  J.  Watson  Research  Center 

Yorktown  Heights,  NY  00598 

George  S.  Fishman 
Curriculum  in  Operations  Research 
&  Systems  Analysis 
Phillips  Annex 

University  of  North  Carolina,  Chapel  Hill 
Chapel  Hill,  NC  27514 

Toke  Jayachandran 
Naval  Postgraduate  School 
Monterey,  CA  93940 

Dr.  D.  L.  Iglehart 
Department  of  Operations  Research 
Stanford  University 
Stanford,  CA  94305 

Dr.  Averill  M.  Law 
1513  University  Avenue 
Department  of  Industrial  Engineering 
University  of  Wisconsin 
Madison,  WI  53706 

Dr.  Peter  A.  W.  Lewis 
Naval  Postgraduate  School 
Department  of  Administrative 
Sciences 

Monterey,  CA  93940 


Naval  Ocean  Systems  Center 
Code  16 

San  Diego,  CA  92132 

Newport  Laboratory 

Naval  Underwater  Systems  Center 

Newport,  RI  02840 

New  London  Laboratory 

Naval  Underwater  Systems  Center 

New  London,  CT  06320 

Naval  Weapons  Laboratory  (Code  K-30) 
Dahlgren,  VA  22448  2  copies 

Dr.  Donald  P.  Gaver 
Naval  Postgraduate  School 
Operations  Research  Department 
Monterey,  CA  93940 

Mr.  Carl  Bates 

MATH-STAT  Team,  MOCA-MRA 

U.S.  Army  Concepts  Analysis  Agency 

8120  Woodmont  Avenue 

Bethesda,  MD  20014 


Naval  Air  Development  Center 
Warminster,  PA  18974 

Naval  Research  Laboratory 
Washington,  D.C.  20375 
Code  2627  (1  copy) 

Code  5308  (1  copy) 

Code  7932  (1  copy) 


-2- 


END 

DATE 

FILMED 

7,-81 

DTIC 


